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The Experimental Approach to Aircraft 
Structural Research’ 


The Fifteenth Wright Brothers Lecture 


P. B. WALKER, M.A., PH.D., F.R.Ae.St 
Royal Aircraft Establishment 


SUMMARY 


Research in aircraft structures is constantly bringing forward 
character. Many of these 
paper 
work 


new problems of a fundamental 
respond to systematic experiment, and the 
the experimental approach in terms of recent 


in Great Britain and particularly at the Royal Aircraft Establish- 


presents 
done 


ment. 


One aspect of research considered is the testing of small speci 
mens made of Xylonite. Another is the strength testing of actual 
aircraft, which opens up several fields of research activity for dis 
Design studies of new testing equipment for large air- 
Consideration 


cussion. 
craft and for pressure cabins are also described. 
is given to the measurement of external forces in flight and the de- 
sign of measuring instruments, including a description of the new 
» counting accelerometer recently developed by the R.A.E._ Fi- 
© nally, structural fatigue is discussed as being one of the most diffi 
cult subjects with which the structural specialist has yet had to 


deal. 


INTRODUCTION 


nad STRUCTURAL RESEARCH now Covers a 
wide field.! For the purpose of this paper it 
is proposed to interpret it in the somewhat restricted 
Sense, which excludes the dynamic structural work 
such as transient loading, vibration, and aeroelasticity. 
| By selecting for discussion the experimental approach, 


Presented before the Institute of the Aeronautical Sciences in 
the U.S. Chamber of Commerce Auditorium, Washington, D. C., 
© December 17, 1951. 

* I wish to thank the Chief Scientist of the British Ministry of 
Supply and the Director of the Royal Aircraft Establishment for 
© permission to present this paper, though the opinions expressed 
are my own. I would also like to express my thanks to my 
colleagues at the R.A.E. who have assisted me in its preparation, 
and also those who have allowed me to make free use of their un- 
© published work. 

t Head, Structures Department. 


moreover, as distinct from what might be referred to as 
the theoretical approach, the field is still further re- 


stricted. 


The words ‘“‘experimental’’ and ‘‘theoretical,’’ how- 
ever, as used in this connection, may be somewhat mis- 
leading. The experimental approach is more scientific 
than some suppose, and it does not lead anywhere 
unless it is closely related to theory. Conversely, 
much of what is termed ‘“‘theory of structures”’ is really 
a particular branch of applied mathematics, with basic 
theoretical concepts already established. In some re- 
spects, therefore, it would be more appropriate to con- 
trast the physicist’s and the mathematician’s approaches 
rather than the experimental and the theoretical. 


The experimental approach to structural research, 
therefore, is here regarded as essentially the physicist’s 
approach. The experimental work described is aimed 
primarily at the fundamentals of aircraft structural 
engineering. Ad hoc tests and experiments that con- 
tribute to the development of particular aircraft types 
are not under consideration unless they have a general 
significance or need to be studied from the standpoint 


of method and technique generally. 


The particular subjects chosen for discussion are 
presented mainly as a study in experimental method. 
At the same time, the opportunity is taken to describe 
recent work done in Great Britain and particularly at 
the Royal Aircraft Establishment (RAE). It is hoped 
that some of the information will appear fresh to those 
who receive it; where this is not so, it is still hoped that 
there will be some interest in seeing familiar problems 


treated from a different point of view 


145 
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(1) THE EXPERIMENTAL STUDY OF STRUCTURAL 
PRINCIPLES BY MEANS OF SMALL SPECIMENS 


(1.1) Introductory 


To the aircraft structural engineer there is much 
satisfaction to be gained from small-scale experimental 
work done in a laboratory. This arises partly from 
economic considerations, but there is no doubt that the 
main benefits arise from the ease and speed with which 
small-scale work can be accomplished. New ideas can 
be quickly tried out experinientally, in contrast to the 
long delays that inevitably occur in full-scale structural 
work. 

In the technique now to be described, small and sim- 
ple structural models are made in ““Xylonite’’ and tested 
under applied load to determine the distribution of 
stress and strain. The work is aimed mainly at the 
study of structural principles and is not to be confused 


with attempts to produce exact scale-models. 


(1.2) Origin of the Method 


For some years, various research workers have been 
experimenting in the use of small-scale models for in- 
vestigating the properties of full-scale aircraft structures. 
Their objective was to produce in easily-workable 
material an exact reproduction of a structure that could 
conveniently be tested in a laboratory. The proprie- 
tary material ‘‘Xylonite’’ has been much used for the 
construction of such models, and there are now a number 
of skilled craftsmen in Great Britain who can work 
quickly in this material to exacting specifications. 

The value of direct scale-model work, however, is 
now much in doubt. It is difficult, if not impossible, 
to produce a structurally equivalent model on a small 
scale and in quite different material. This work on 
scale models, however, has had one important result: 
it has inspired what is virtually a new method of experi- 
mental structural research and provided for it a ready- 
made technique. 

In this new method, attention is concentrated on 
reproducing a particular structural phenomenon with- 
out detailed representation of the original structure 





| DIRECTION OF LOAD 


Fic. 1. Small-scale representation of engine-loading case. 
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on which it occurs. A small model is made in xylonite 
as before, but from the engineering standpoint it is ay 
extremely simplified version of the original. The 
dimensions and stiffness properties of the simplified 
members, however, are carefully chosen so that the 
behavior of the original structure, in respect of the 
particular phenomenon under review, can be deduced 
with reasonable accuracy from the model experiments, 
The technique for making and testing the simplified 
models is very similar to that used in true scale-model 
work. 

These tests of structural principle have only been in 
use at R.A.E. for a short time, but they have already 
proved their value. Furthermore, progress has been 
unhampered by minor difficulties of technique, owing 
to the experience already obtained in direct scale-model 
work. Xylonite has the advantage that deflections are 
comparatively large and cannot only be measured but 
can also be observed visually. The transparency of the 
Xylonite is a further advantage, since it enables the 
behavior of the inner parts to be observed. Facility 
for direct observation is useful when a phenomenon 
needs to be understood in a broad way before detailed 
investigation begins. There is, however, full scope for 
accurate quantitative work with dial gages and electri- 
cal-resistance strain-gages. 

The method has two distinct applications. Firstly, 
it can be used for investigation of a particular phenom- 
enon occurring on a specific aircraft structure. In 
this connection it is a valuable ancillary aid in the rou- 
tine testing of large structures. Secondly, it can be 
used for exploring general trends in design, and has thus 
a place in the structural research necessary to cope with 
new aerodynamic forms such as sweptback and delta 
wings. It is now proposed to consider examples in 
each of these two general categories. 


(1.3) Specific Application 


The best examples of direct application of the 
method to particular structures arise in the course of 
routine strength Most 


are tested before final acceptance, and from time to 


testing. aircraft nowadays 
time a structure behaves in a way that is entirely un- 
expected. Premature failure may sometimes be ex- 
plained by a mistake in calculation or by defective 
material, but frequently the explanation goes much 
deeper. The underlying theory or the basic assump- 
tions may be at fault, and it is in such cases that a 
small-scale test of structural principle is called for. 

A good example occurred only recently. Premature 
failure occurred on the wing of a wing-engined aircraft 
in the neighborhood of the attachments to one of the 
The weakness could not be explained in 
fundamental 


engines. 


any simple way, and some _ more 


explanation was sought. The most obviously unusual 
features were (a) the positions of the engines, which 


were much further forward than on any previous air- 
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craft tested and (b) the lightness of the ribs compared 
with those of previous aircraft. 

A Xylonite test-specimen was made reproducing a 
simplified version of the main structure of the original 
wing and engine mounting as shown in Fig. 1. The 
specimen was essentially a long rectangular box re- 
inforced by ribs and stringers and carrying a projection 
representing the engine. The resultant engine load 
js indicated by an arrow in the picture. The structure 
was made as simple as possible, but dimensions were 
chosen to ensure that the phenomenon could be in- 
vestigated quantitatively (to a reasonable degree of 
accuracy) as well as qualitatively. Deflection gages 
and strain gages were installed in the usual way (Fig.2), 
an identical second model being used for the dummy 
strain gages required for temperature compensation. 

In the early stages of the experiment, it became clear 
that something unexpected was happening, and the 
investigators were not left long in doubt as to the ex- 
planation of the premature failure of the actual wing. 
In the original calculation of the stresses, the customary 
allowance had been made both for the engine load 
and also for its concentrated application. The theory 
used, however, did not go far enough. 
assumed for the purpose of analysis that the ribs 
a legitimate working assumption 
In this particular case 


It had been 


were infinitely stiff 
for most structural problems. 
however, the actual flexibility, not only of the rib to 
which the engine was attached but also of the ribs 
outboard, was found to cause additional stress concen- 
trations in the wing which had clearly given rise to 
the premature failure. 

The error was shown by these experiments to be 
approximately 15 per cent which almost exactly ac- 
counted for the discrepancy found in the test of the 
original structure. Thereafter, remedial action was 
straightforward. The front spar of the wing was 
strengthened locally to the extent indicated by the 
experiments. 

It is of interest to remark that this phenomenon 
had been noticed for the first time partly because the 
engine was further forward than hitherto and partly 
because the designers had been extremely careful in 
reducing structure weight to the minimum—working 
to what had then been the accepted theory for wing 
design. Now that the explanation of the discrepancy 
has been obtained, a recurrence of the trouble on future 
aircraft is unlikely. 

In the light of this experience it is unnecessary to 
stress the value of this method of investigation, but 
there is one other point to emphasize—the Xylonite 
models were made and the essential preliminary test 


completed in about 3 weeks. 


1.4) The Study of General Trends 


The application of the method to the study of general 


trends in design presents a much less clearly defined 
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Fic. 2. Strain-gage system for engine-loading case 


field of activity. It is necessary to bear in mind that 
the method has not been in use very long and that its 
full possibilities cannot yet have been exploited. 
There are on record, however, a few general investi 
gations on sweptback and delta wings which are of 
interest. 

The first example treats the sweptback wing as a 
“conventional” straight wing merely inclined back 
wards (Fig. 3). There is thus introduced the additional 
wedge of wing structure indicated on the photograph. 
Simple xylonite models of both “‘straight’’ and swept 
back versions were tested and the results compared. 
The experiment showed that for the sweptback wing 
working inwards toward the root—there is a shedding 
of load (shear and bending moment) from the front 
spar, with a corresponding increase in the rear spar. 
In other words, the loads tend to take the shortest 
path to the fuselage. The extent of load transfer 
depends on relative sizes of the members and angle 
of sweepback, but it can be determined in any specific 





Fic. 3 


Specimen representing simple sweepback 
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Fic. 4. Sweptback wing—-comparison of 


structural’ ribs. 


“aerodynamic’’ and 


case with a fair degree of accuracy. For the average 
type of sweptback wing with, say, 40° sweepback, the 
increased load in the rear spar at the root is about 15 
per cent. 

The second example concerns the direction of the 
main ribs. As a structure analogous to a ‘‘straight’’ 
wing structure the ribs should lie approximately at 
right angles to the spars as shown on the left of Fig. 4. 
Considerations of aerodynamic efficiency, however, 
generally prescribe that the ribs should lie in the direc- 
The 
aerodynamic value of such ribs is still arguable in some 


tion of flight, as shown on the right in Fig. 4. 


cases, but it is clearly necessary that their associated 
structural penalties be properly assessed. To deal 
with this problem a pair of boxes of simple construction, 
reproducing the main characteristics of the structural 
part of a swept wing, was made as illustrated, with 
The 
showed that the “structural” ribs (perpendicular to the 


ribs in the two directions, respectively. tests 


spars) are superior to the ‘“‘aerodynamic”’ ribs (parallel 





Delta radial 


comparison of fore-and-aft and 
trailing ribs 


wing 
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to the line of flight), both for stiffness and strength, 
The difference, however, is not great, and, in general 
the increase of strength and stiffness to be obtained by 
inclining the ribs to the line of flight is of the order of 
5 per cent. 

The third example concerns the trailing portion of 
On the left 
of Fig. 5 the trailing ribs are shown parallel to the direc. 
tion of flight. 


the delta wing, aft of the main structure. 


The alternative is to support the trail. 
ing portion of the wing by members directly attached 
to the root as shown on the right of Fig. 5, where the 
ribs are spread out fanwise from the root of the rear 
spar. To deal with this problem, simple specimens 
with the two types of trailing ribs were made as before 
and tested. It was found that for the same structure 
weight the fanwise ribs are slightly superior structurally 
to the fore- and aft-ribs, both with regard to strength 
and stiffness. 

These three examples of preliminary work with Xylo 
nite models indicate that there is considerable scope 
for assisting the aircraft designer when he has to decide 
the general layout of his design. They also indicate 
that it may often be worth while making and testing 
a simple Xylonite model of a structure whenever it is 
proposed to depart appreciably from conventional prac 
tice. 


(2) STRENGTH TESTING OF COMPLETE AIRCRAFI 


(2.1) Introductory 


From the testing of small ‘‘models’’ we now pass on to 
consideration of the very different task of strength test 
ing real structures. Most new designs of aircraft nowa 
days are strength-tested in more-or-less complete form 
The 
purpose of each test is to check whether the aircraft 


before being released for service.* immediate 
meets the standard of strength previously laid down 
as being necessary for safety and operational useful- 
ness. Should any structure fail prematurely, however, 
there is then the further purpose of determining the 
nature and situation of the weakness so that it can be 
rectified. 

The work as a whole opens up several fields for re 
search activity, in spite of its superficially routine charac- 
ter, and it is proposed to consider three of the more 
general aspects of routine strength testing which have 
received particular attention recently in Great Britain. 


(2.2) Safety and Structural Efficiency 


From time to time it is necessary to consider whether 
the effort being expended on major strength testing 1s 
justified from the standpoint of safety and efficiency 
of the individual aircraft types concerned, and whether 


* Actually, it is common practice with large aircraft to test 
wing and fuselage systems separately. It is then usual to have 
one system complete and accompanied by sufficient of the other 
for accurate load representation, and there is no departure in 


principle from the testing of the complete aircraft 
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the work as a whole is being directed to best advantage. 
The best indication of the need for major strength 
tests is to be obtained from a study of past records. 
4 review of results obtained in the past 10 years or so 
was recently made in Great Britain.2 Among infor- 
ation of various kinds collected, the most comprehen- 


m 
concerned 24 wing systems, and it is to these that 


sive 
the following analysis refers. 

Initial results for the aircraft are indicated diagram- 
matically in Fig. 6. Each circle in the diagram rep- 
resents a particular aircraft, and the strength attained 
on test is plotted as a percentage of the design require- 
Of the 24 wing systems, ten are seen to be on 


ment. 
In other words 


the wrong side of the 100 per cent line. 
more than two-fifths of the aircraft as delivered failed 
to meet the standard of strength their designers had 
set themselves, while roughly one-fifth were completely 
unacceptable on safety grounds. This is regarded as 
clear proof that strength testing of all basically new 
designs is necessary, and that neither structural theory 
nor design skill can be relied upon entirely without 
testing as a final check. 

The outlook brightens considerably, however, when 
events following these initial failures are investigated. 
Of the ten understrength aircraft, nine were strength- 
ened locally in the light of information yielded by the 
initial test, and the test on each was then repeated. 
It is significant that in the great majority of cases the 
alterations were minor and the corresponding increase 
of structure weight almost negligible. This process 
of modification and retesting, moreover, was carried 
out once again in two cases where results attained in 
the second test were still unsatisfactory. It was also 
applied to a few designs that, though meeting the de- 
sign requirement, showed good prospects of useful 
extra strength being attained at small cost in structure 
weight. 

Fig. 7. shows the result when one alteration is allowed 
if required, or two alterations in rare cases. It is seen 
that there are now only four aircraft (or one-eighth 
of the total) on the wrong side of the 100 per cent line 
and that only one is less than 90 per cent. Even this 
result requires some qualification: it was not always 
considered necessary to retest aircraft achieving 90 
per cent or more, and two were retested only once 
while one was not retested at all. It is of interest to 
note, moreover, that even the one aircraft that did not 
respond to two minor alterations had fundamentally 
a sound structure; and, although something more than 
minor alteration was necessary, it was brought up to 
standard without much difficulty and at a surprisingly 
small increase in structure weight. 

These results attained after minor alterations are 
giving support toa philosophy of design which many have 
followed in the past intuitively. It is maintained that 
the most efficient designs—that is, designs with the 
lightest structure for a given strength—are most likely 
to be obtained when the designer deliberately takes a 
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Distribution of strength of wing systems as delivered. 


few chances and relies on strength tests to detect any 
weakness. It can be argued that, in order to have pre- 
vented any one of these premature failures, it would 
have been necessary to have disposed extra structure 
weight fairly liberally over much of the structure and not 
merely in one place, since the precise form of failure 
was unpredictable. The case for proceeding on these 
lines, moreover, is further strengthened by the fact that 
many of the structures that failed prematurely in the 
first test were ultimately among the most efficient of 
them all. 

This principle can, of course, be carried too far. 
There is always a possibility that too much confidence 
might produce a structure that is altogether too weak 
and beyond simple rectification. More likely, however, 
is that serious delay may be occasioned through numer- 
ous petty modifications and retests while an aircraft 
is waiting for the results before it can be allowed to fly. 
The solution to the problem thus presented may there- 
fore be some kind of compromise, and there is a good 
case for supposing that something close to the best 
compromise is already being attained in current prac- 
tice. 

Before discussion on this subject can be closed, it 1s 
necessary to consider to what extent such deduction 
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Distribution of strength on completion of modification 
program 


Fic. 7. 
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A typical test to destruction. 


Fic. &, 


from past experience can legitimately be applied to the 
future, and’ generally to temper statistical analysis 
with engineering judgment. The number of aircraft 
showing excessive strength gave rise to some concern 
until it was found that they were mostly early types 
created when the need for weight economy was not 
realized so acutely asitisnow. There are good grounds 
to believe, in fact, that future aircraft are most unlikely 
to be made excessively strong. On the contrary, there 
are indications that the proportion of premature failures 
on first test may well increase and even that the trend 
in this direction has already begun. This is consistent 
with the principle of modification and retest already 
mentioned and gives grounds for satisfaction rather 
than alarm, particularly since there is no evidence of the 
principle being carried to excess. 


(2.3) Contribution to Fundamental Knowledge 


The satisfactory clearance of specific aircraft types 
must always remain the first objective in major strength 
testing carried out as a routine. Nevertheless, the 
tests are not being used to full advantage unless they 
are made to contribute to the store of general structural 
knowledge and to improve the breed of aircraft gener- 
ally. Major strength tests are costly in time, money, 
and effort; and for this reason they can rarely be under- 
taken for research purposes alone. Every possibility 
must therefore be taken to obtain from the ad hoc 
tests as much fundamental and general information as 
possible. In the course of years a great variety of 
failures occur, since each test is normally taken to de- 
struction (Fig. 8) and many new and unsuspected modes 
of structural behavior are revealed from time to time. 

At the R.A.E. there now exists a small team of special- 
ists charged solely with the task of extracting funda- 
mental lessons from these major strength tests. Its 
members are not directly concerned with the ad hoc 
purposes of the test, though they naturally collabo- 
rate closely with those who are. In event of premature 
failure, moreover, they often help considerably in the 


diagnosis of the cause. Their main instrumental aid 
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is the electrical-resistance strain gage, and strain gages 
are distributed liberally over parts of the structyp 
deemed to be especially vulnerable or of special scientif 
interest. 

It is not possible, of course, to narrate here all the 
conclusions that have been drawn from such investi. 
gations. Often the work does not appear as an ae. 
count of a particular test but as a study of a specifi 
structural problem supported by further theoretica| 
conclusion, however, 


work. One 


worth mentioning. 


general appears 
This is that the most common 
single cause of premature failure is a ‘‘stress diffusion 
problem. In other words, while ample material may 
be present to carry the internal load in a particular 
region, some part of the structure carries a greater 


share than the designer intended or expected. 


(2.4) Technique of Force Control 


Another field of study related to major strength 
testing concerns the testing technique, which needs 
always to be kept up to date and in line with modem 
requirements. The most important branch of testing 
technique is the control of the applied forces. This 
presents a problem much more serious than is generally 
realized and for which it is probably true to say that 
no really satisfactory solution has yet been found. 

A few years ago the writer undertook an investiga 
tion into the principles of force control.’ The con 
clusion was reached that the main difficulties arise 
from the interaction of the various strainers, as the 
While the 


number of strainers is reduced, by means of the well- 


force-controlling units are generally termed. 


known lever system (Fig. 9), to a small fraction of the 
number of forces applied to the structure, there are 
still several in a test of any magnitude. 

The investigation also led to the further conclusion 
that there are two basic types of strainer according to 
whether force is controlled directly or indirectly through 
displacement. The strainer 
must be mechanically reversible, so that it is free to 
move to and fro without change in the force, as the 


direct force-controlling 





Use of lever systems with several strainers. 
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structure distorts under the action of the other strainers 
in the system. Such a strainer may be pictured as a 
frictionless hydraulic ram supplied with fluid at a con- 
stant predetermined pressure (Fig. 10, /eft). Con- 
versely, the displacement-controlling strainer must 
clearly be irreversible and is best represented by a screw 
jack (Fig. 10, right). 

Unfortunately, neither basic type of strainer is satis- 
factory in any simple form. The reversible strainer 
may give rise to dynamic effects and to instability of 
the whole system,‘ but its chief disadvantage is the extra 
damage it may do to the structure after initial primary 
failure. The irreversible (screw-jack) type of strainer 
has none of these disadvantages, but it does not give 
complete control of the force system, since the relation 
between forces and displacements depends on the charac- 
teristics of the structure and is never known precisely. 

All practical types of strainer are then seen as com- 
binations of these two basic types produced with a 
view to obtaining, so far as possible, the advantages 
of both and rejecting the disadvantages. The R.A.E. 
compound strainer (Fig. 11) is a simple illustration of 
this principle. An irreversible screw-jack is used but is 
suspended from a lever that is supported by a small hy- 


draulic ram. The travel of this lever, however, is re- 
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stricted to something small by stops. Thus the strainer 
is irreversible in respect of large movements and re- 
versible in respect of small. 

The leverages and ram diameters are actually ar- 
ranged so that all the rams in a complete system of 
strainers require the same pressure for the prescribed 
loading ratios, and they are interconnected by pipe 
lines. In operation, the main straining is done by 
screw-jacks and a coarse adjustment of forces obtained. 
Final adjustment is then obtained automatically on the 
hydraulic system. 

Fig. 12 shows the strainer in a practical form. 
are two coupled screw-jacks, handed to eliminate torque, 
As a result of 


There 


which are driven by electric motors. 
the investigation that has been made, it is proposed to 
develop this type rather than proceed on entirely new 
lines. In the existing system, contro] is automatic 
for the hydraulic part of the system, and all the screw 
jacks are controlled from a central control room, though 
each in some measure individually. Successful experi- 
ments have been made on a small scale to eliminate 
individual control entirely and to give a completely 
automatic system. The full application of this method 
to regular testing is an immediate aim for the future. 


(3) DESIGN oF TEST FRAMES AND PRESSURE TESTING 
EQUIPMENT 


(3.1) Introductory 


The next subject for consideration is the design of 
the testing gear used in the major strength tests just 
described and especially the test-frame structure. 
This is an important responsibility of the aircraft 
structural engineer and takes him quite outside his 
usual field. The design of a structure that is to be 
used for testing another structure is decidedly heavy 
engineering and as unlike aircraft engineering as it 
could possibly be. 

Interest in design of test frames at the present time 
is inspired by two testing requirements that are com- 





Electrically-driven compound strainer (R.A.E. type) 
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End view of R.A.E. ‘‘Cathedral’’ test frame. 


Fic, 13. 
paratively new. One is for the testing of larger air- 
craft generally, and the other for testing sweptback 
and delta types. The largest aircraft that can be tested 
at R.A.E. at present has a 120-ft. span and must be of 
the ‘“‘straight-wing’’ type. An aircraft not coming 
within this capacity has to be tested as a reduced-scale 
version or by improvised methods. 

In addition to the prospect of extending 
methods of testing, consideration has to be given to 
This is a new develop- 


‘ 


‘orthodox”’ 


the testing of pressure cabins. 
ment requiring the cabin to be filled with water and 
totally submerged in water. 


(3.2) The Universal Test Frame 


There are wide differences of opinion as to the best 
type of testing structure for major strength tests. 
Many systems call for a good deal of improvisation, 
the testing gear being designed to suit the particular 
aircraft that is to be tested. As is well known, the 
R.A.E. have kept consistently to the “‘universal test 
frame,’ which enables any aircraft within specified 
dimensional limits to be tested without any specia] 
constructional work apart from the lever systems. 
Two advantages normally conceded for the universal 
test frame are saving in time and labor during the proc- 
ess of rigging for a test and during the actual testing. 
These are regarded as over-riding considerations for 
a central establishment, though not 
necessarily for design firms concerned with one partic- 
ular aircraft during its design stages. A further 
advantage claimed is greater accuracy of measurement 


experimental 


and greater precision of load control through use of 
equipment that is permanently installed. On the other 
hand, first costs are high, and it has been maintained 
that the system is not so flexible as others involving 
a greater degree of improvisation. 

Before discussion of prospective developments in 
this field it is desirable to have on record a brief descrip 
tion of a typical universal test frame. Essentially, 
this is a structural box with open ends. Fig. 15 
shows an end view of the ‘“‘Cathedral,’’ the largest test 
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frame in use at R.A.E. The top of the frame consists 
of a series of bridges, each supported by a pair of two. 
wheeled carriages which can be moved along raijis 
These 
bridges also support rails on which rest the bogies 


traversing the whole length of the test frame. 
that carry the upward-loading strainers. The two 
degrees of mobility thus obtained are utilized for initia] 
positioning only, and, subject to minimum spacing limi. 
tations, enable the forces to be applied in any desired 
position. 

The comparatively few downward forces are applied 
in three possible ways: by strainers attached to the floor 
structure, by levers connected to the upper strainers 
by tie rods, or by dead anchorages. 


(3.3) Design Study of a Large Test Frame 


In discussing the design of larger test frames it js 
proposed to consider only the requirements for wing 
systems, since fuselages in general are then covered 
automatically. The choice of an upper limit of size 
depends upon three considerations. There is first the 
probable sizes of future aircraft generally, which can- 
not of course, be predicted with any degree of certainty, 
Next there is the possibility of not testing aircraft of 
extreme size, since it has been argued by some that 
testing of such aircraft is becoming uneconomic and 
must be waived. Finally, there is the possibility of 


testing specially constructed reduced-scale versions 


of very large aircraft. This last method has been 
employed in Great Britain with apparent success and 
is undoubtedly better than no test at all; but it is 
costly in time and effort, and there is always an element 
of doubt as to whether the reduced-scale structure and 
the test itself truly represent reality. 

Before a proper decision on the maximum size of 
test frame can be reached, however, it is desirable to 
know what limitations, if any, are imposed by considera 
tions of the practicability of design, construction, and 
operation of a universal test frame considerably larger 


than any in existence. Accordingly, a design study of 
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a large test frame was recently prepared without prej- 
udice to any later action that might be taken. 

The internal dimensions of the test frame in this de 
sign study are 200 by 100 ft. Any aircraft, irrespective 
of shape, with a wing system not exceeding these dimen- 
sions can be tested, though it is not expected that 
both extreme dimensions will ever occur simultaneously. 
Furthermore, a straight-wing aircraft of span es great 
as 230 ft. can be tested provided the extreme wing 
tips are removed, a practice that should not lead to 
appreciable error. In consequence, existing aircraft 
like the Brabazon can be accommodated. 

Fig. 14 shows a scale model of the test frame produced 
by this design study, with part of the building and some 
of the bridges removed for clarity. A model of a 
large sweptback wing system is also shown in its 
testing position. 

In the creation of this design the mobile bridges 
provided the greatest problem. As is well known, 
bridges increase in weight, depth, and width out of 
proportion to increase of span. The bridges shown 
weigh SO tons, which is not considered excessive. In 
order to obtain sufficiently close loading lengthwise 
along the frame, however, it is necessary to have two 
sets of rails with the bridge carriages arranged alter- 
nately. The most critical design factor is the bearing 
pressure on the rails, and it is considered that the 
present arrangement works close, though not too close, 
to the practical limit. 

The effect of length of the test frame (corresponding 
to wing span) is much less serious than the effect of 
breadth (corresponding to wing chord). To 
cope with the tip deflections of a 200-ft. span aircraft, 
however, requires a clearance height of 50 ft. measured 
from floor to the lowest point on the bridges. 

The building to house the test frame also requires 
serious consideration and, treated separately, costs 
In this scheme, 


gross 


about as much as the test frame itself. 
however, the building is constructed integrally with 
the frame (Fig. 15) partly for economy but also to sim- 
plify the design of the foundations. Though presenting 
a level surface, the floor has beneath it a stout struc- 
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ture incorporated in the foundations and forming an 
integral part of the test-frame structure. The height 


of the building at its highest point is about 110 ft. 


(3.4) Strength Testing of Aircraft Pressure Cabins 


From one point of view the strength testing of pres- 
sure cabins is simple because it is merely necessary to in- 
crease the pressure until failure occurs. Unfortunately, 
however, when the pressure is applied in the normal 
that is, to the air within the cabin—the effect 
In a calculation recently made by 


way 
can be disastrous. 
the writer, it was shown that the effect of bursting of 
a fair-sized pressure cabin filled with air is roughly 
equivalent to the bursting of a 100-Ib. bomb.° 

The testing of a pressure cabin by internal air pres 
sure thus involves grave risk to personnel carrying out 
the test, and certain damage to equipment and build- 
ings. Itis, of course, conceivable that these risks might 
be lessened by building protective walls, but there 
still remains a major difficulty. The cabin structure 
is usually blown into small pieces so that there is little 
hope of diagnosing the primary weakness and no chance 
at all of a repair and retest. 

The solution lies in using water instead of air. 
Water inside a cabin, however, introduces a static 
head that cannot be ignored, and the total weight 
may be so great that it cannot be supported without 
introducing extra loads of appreciable magnitude. 
The whole cabin is therefore submerged in water as 
well as filled with water. 

The calculations already mentioned show that the 
available energy stored in the compressed water is of 
the order of one ten-thousandth part of that stored 
by air under the same conditions. Furthermore, the 
presence of the water lessens damage still further by 
damping ‘down any tendency to violent motion. In 
consequence, failures in a water test are usually local 
ized even to a greater degree than in ordinary strength 
testing, and the diagnosis of the primary weakness is 
rendered easy. The difference between the effects of 
air and water is well illustrated by the results of some 
comparative experiments. Two identical cylinders 
were made (Fig. 16) about 9 ft. long and 3 ft. in diameter. 





Model pressure cabin before test (9 by 3 ft 
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Fic. 17. Model pressure cabin after air test. 





Fic. 18. Model pressure cabin (ends removed) after water test 
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R.A.E. strength-testing tank for pressure cabins 


In an air test the first specimen was blown to pieces, as 
shown in Fig. 17. A corresponding test in water, 
however, gave results as shown in Fig. IS. In this 
last picture, the undamaged ends have been removed 
for inspection, and it is seen that the only damage is 
failure of the rivets along the seam. 
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One of the main characteristics of water pressure 
testing is its extreme simplicity. Care has to be taken 
that little or no air is trapped in the cabin. Pressure 
is produced by a simple hand pump. Fig. 19 shows 
a testing tank, recently constructed at R.A.E., in the 
The tank is 
15 ft. deep, half of this being below ground level, ang 
pressure cabins of most aircraft can be accommodated. 


course of some preliminary experiments. 


(3.5) Design of a Fuselage Test Frame for Combined 

Pressure and Structural Loading 

The foregoing discussion on testing of pressure cabins 
relates to pure pressure tests. In the absence of any- 
thing better these serve a useful purpose, but in practice 
the cabin structure carries primary structure loads as 
well as pressure loads. A design study has therefore 
been made of a combined test frame on the lines already 


described for the normal test frame. 


Fig. 20 is a photograph of a transparent model of the 
test frame produced in this design study. Essentially, 
it is a normal test frame, with proportions and dimen- 
sions appropriate to the testing of complete fuselages, 
but it is almost completely submerged in a tank of 
water. Only the bridges and upper strainers, a few 
of which are shown in the picture, stand clear of the 
water. 

For the downward forces it is not essential to have 
the straining units with their electrical gear under water, 
though this is not necessarily considered impracticable. 
Where not produced by dead anchorages, the few forces 
normally required can be transmitted through levers 
and tie rods to the upper strainers. 

The length of the test frame is 200 ft., its internal 
width is 45 ft. and the clearance height is 35 ft. The 
ends of the tank can be lowered to give easy access 
during erection of the test specimen. The test frame 
can also be used efficiently for testing fuselages without 
pressure loads and without water being present. 

(4) FLIGHT ACCELERATIONS AND THE COUNTING 

ACCELEROMETER 
(4.1) Introductory 

Experimental work such as that so far described 
deals essentially with internal behavior of the structure, 
Safety 


and structural efficiency, however, depend as much 


with the external forces taken for granted. 


upon the proper assessment of these external forces 
as upon the correct determination of the stresses and 
strains they produce. In consequence, the study of 
external and internal loading have become closely 
integrated and both require to be investigated by the 
structural specialist. 

In a somewhat similar way, the responsibilities of the 
structural specialist have extended in another direction 

to include the design of instruments of measurement. 
The planning of any major experimental investigation 


is becoming increasingly a matter of deciding what 
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can be measured and how it can be measured. The 
result is that experimental structural research and 
instrument design are in many instances inseparable 
parts of the same problems. 

The new counting accelerometer about to be de- 
scribed illustrates both these trends. On the one hand, 
it is intended for measurement of the accelerations 
occurring in flight, upon which depend the main load- 
ing conditions for the aircraft. On the other hand, it 
has been invented by a structural specialist* as a stage 
in the study of gusts and their effects on an aircraft 
structure. 

This instrument is described in some detail as having 
prospects of becoming widely used and also as provid- 
ing a good example of the principles of instrument 


design for specific structural objectives. 


4.2) The Counting Principle 


One of the main problems in present-day design of 
recording instruments is to prevent them supplying 
too much information. There are already many in- 
struments that produce a quantity of data greatly ex- 
ceeding what can be properly analyzed and assimilated. 
Whenever this risk is encountered, attempts must 
be made to design an instrument that has in some degree 
powers of discrimination, so that it will either reject 
secondary or unwanted information, or else segregate it 
from that which is essential so that it can be ignored. 

The principle is illustrated in a remarkable degree 
by the lg recorder of United States origin. Ostensibly, 
this form of accelerometer records as a continuous curve 
the variation of speed and acceleration. Actually, 
its purpose is to record only extreme values of speed and 
acceleration, and it does so in a most effective way. 
The recording point, in passing over the middle regions 
of speed and acceleration a great many times, virtually 
obliterates all record of them. The peaks of speed 
and acceleration, which occur infrequently, are then 
left prominently displayed as projections from the cen- 
tral mass. 

The counting accelerometer bears a family resem- 
It also measures acceleration, 
Its purpose is to 


blance to the Vg recorder. 
though it does not measure speed. 
record accelerations over the full range, however, and 
especially the fluctuations at lower levels which are one 
of the primary causes of structural failure in fatigue. 
With this objective there is again the danger of pro- 
ducing too much information. A continuous record 
of acceleration would require too many man-hours 
to read and interpret when the instrument was in 
generaluse. In principle, the necessary simplification is 
achieved by simply counting the bumps as they occur in 
flight, though classifying them according to their sever- 
ity. In practice, this process is modified slightly for 
convenience in instrument design. A range of values 
of acceleration separated by convenient intervals is 


* James Taylor, of the R.A.E. 





Transparent model of submerged test frame for fuse 
lages under combined loading. 
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first decided upon, and the accelerometer is then made 
to count separately the number of times each of these 
values are attained. Thus the records for a flight of 
any length are reduced to a few figures, each correspond- 
ing to a given level of acceleration. 


Once the principle of counting is accepted, however 
a major problem is presented. Unless special precau 
tions are taken, the accelerometer will count not only 
the accelerations corresponding to overall motion of 
the aircraft but also the accelerations caused by air- 
frame vibration. The overall accelerations are wanted 
as being a measure of the forces acting on the aircraft. 
The vibration accelerations have no such significance 
and are unwanted in this connection. If the two are 
counted indiscriminately, there is no possibility of 
Hence, the first requirement for a 
shall completely 


later separation. 
counting accelerometer is that it 
reject accelerations caused by vibration while truly 
recording those produced by overall motion of the air 
craft. 


4.3) The Accelerometer Mechanism 


In the design of any accelerometer for recording ir 
regular fluctuations of acceleration such as occur in 
flight the main problem is to render the readings of 
acceleration effectively independent of rate of change 
of acceleration. In practiee, it is usual to express 
how nearly this independence is achieved in terms of 
response to a sinusoidally varying acceleration of con 
stant amplitude, applied at various frequencies. Thus 
the ideal is a uniform amplitude of response for all 
frequencies of the applied acceleration, though this, 
of course, is never completely attainable. In practice 
a working compromise is usually obtained over a limited 
range by an appropriate choice of the natural frequency 
of the accelerometer and by heavy velocity damping. 


Fig. 21 shows the simple accelerometer represented 


by a mass suspended on a spring. The response curve 


+ This is also at least a desirable characteristic for any other 
form of accelerometer, though vital to the counting accelerometer 


for the reasons stated. 
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Fic. 21. Simple accelerometer and its response curves 
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Fic. 23. Diagrammatic representation of actual accelerometer. 


for a small degree of damping then has the shape of 


the curve shown dotted in the diagram, the peak cor- 
With heavy 
damping the curve becomes much flatter, as shown by 
the full line. 
used for many purposes as a basis for accelerometer 


responding to the natural frequency. 
Such a fiattened response curve can be 


design, but, as will be shown, it cannot be used for 
the counting accelerometer or for any accelerometer 
that has to reject all accelerations due to air-frame 
vibration, 
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A study of the overall accelerations in flight shoy 
that they can be recorded with reasonable accuracy 
by any accelerometer that has a reasonably uniforn 
response curve for a frequency range from zero up to 
about one and one-half times the fundamental naturaj 
frequency of the wings—usually about 4 or 5 per see. 
On the other hand, the significant air-frame vibrations 
have frequencies exceeding about five times the wing 
frequency. The response requirement for the counting 
type of accelerometer is thus a level response curye 
from zero to one and one-half times the wing frequency 
with virtually no response at five times the wing fre 
quency and above. 

In order to meet this requirement, the simple 
accelerometer—a weight on a spring—is used with 
out any damping, but the weight is coupled to what js 
virtually a second accelerometer (Fig. 22) having a 
somewhat higher natural frequency. In order thet 
this second accelerometer respond only to movement of 
the mass in the first accelerometer, however, and be 
unaffected directly by acceleration of the aircraft, 
it has a rotational inertia instead of a translational one, 
The angle of rotation is then taken as the measure oj 
acceleration of the aircraft. 

The curves in Fig. 22 show the response curves for the 
compound system when the various parameters are 
suitably chosen. The dotted curve corresponds to a 
small degree of damping and has the two peaks appro- 
priate to two degrees of freedom. When heavy damping 
is applied to the rotary inertia, the response curve 
becomes the full line, and this is the one used for the 
counting accelerometer. It is seen that there is a 
fairly uniform response in the lower frequency range, 
with a very low response in the higher frequency 
range. 

The mechanism is shown diagrammatically in Fig. 
23. The translational inertia is suspended on a canti- 
lever spring and is connected through spiral springs 
to a continuous band that drives the rotary inertia. 
A recording pointer is attached to this rotary inertia 


to indicate the acceleration. Fig. 24 shows the actual 





Back view of accelerometer with cover plate removed. 
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mechanism. Damping is by eddy currents produced 
in an armature (not shown), which is driven by the 
rotary inertia through a suitable high-speed gearing. 


4.4) The Counting Mechanism 

The counting is done by the simple mechanical device 
shown in Fig. 25. The rotational inertia in the center 
carries a flexible arm or pointer (the second arm should 
be ignored at this stage), and as this swings in either 
direction from the position appropriate to steady flight 
it passes Over the series of ratchet wheels to be seen 
in the picture. These are spaced at intervals around 
the circle to correspond to intervals of, say, 0.2g 
acceleration and are each connected to simple counters. 
A ratchet and its counter are not moved when the pointer 
moves outward, but when it returns—as it inevitably 
must—the ratchet is gathered up one tooth and the 
counter registers one digit more. 

The second pointer is introduced purely as a matter 
of design convenience to give suitable spacing of the 
ratchet wheels. Both operate in the same way, but 
for overall recording intervals of O.lg one takes the 
odd multiples and the other takes the even. For weight 
economy, it will be noticed, unnecessary counter cylin- 
ders are avoided and the number of digits are reduced 
for high values of acceleration. 

This counting mechanism incorporates two features 
that, though apparently trivial, have fundamental 
significance. In the first place, if the counters were 
to register strictly every time a specified value of acceler- 
ation was passed in both directions, the results could 
be seriously misleading. The acceleration might hover 
about the specified value with the counter registering 
an indefinite number of times. It is essential, therefore, 
that only changes of acceleration exceeding some mini- 
mum amount be registered. The size of the ratchet 
teeth are chosen to give this latitude, and the needle 
can oscillate anywhere through an angle not exceeding 
that corresponding to one tooth without a count 
being registered. 

In the second place, there is a possibility of a mis- 
count through the inertia of the counter. On the one 
hand, if the flexible pointer is moving rapidly, it may 
be unable to move the counter. On the other hand, 
the counter may be moved, and its inertia then can 
carry it forward to register two counts instead of one. 
A flexible connection is therefore interposed between 
each ratchet wheel and its counters. This sets a defi- 
nite upper limit to the speed at which the counter can 
move, and ensures that it always registers a single 


count. 


4.5) The Counting Accelerometer in Use 


The counting accelerometer has been well tried out 
in flights on different aircraft over many thousands of 
It has given little trouble and so far shows 
Already valu- 


miles. 
every sign of being completely reliable. 
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Fic. 25. Front view of accelerometer showing counting mecha 
nism 


able information has been obtained concerning gust 
loadings which is being applied in the study of fatigue. 

In principle, the accelerometer has one limitation, 
such as is to be expected when information is deliber- 
ately suppressed in order to obtain usable data. For 
completely random variation of acceleration informa- 
tion is lacking concerning the order in which the various 
levels of acceleration are crossed. In practice, however, 
this characteristic does not amount to anything serious. 
The normal fluctuations of acceleration on an aero- 
plane follow a fairly simple pattern. 

For calibration and testing the accelerometer can be 
coupled to a continuous electrical recording system 
that gives the order in which the various thresholds 
are crossed. When a few hours preliminary flying 
have been done on a particular aircraft, it is possible 
to make a statistical analysis and to apply the results 
to the accelerometer counts obtained in the ordinary 
way. 

The counting accelerometer with electrical recording 
has incidentally many applications in research experi- 
ments, but the electrical equipment is too heavy for 
routine use or for statistical work generally. The 
mechanical version, on the other hand, weighs about 12 
lb. and is sufficiently small to be installed in aircraft 


engaged in routine flying. 
(5) AIRCRAFT STRUCTURAL FATIGUE 


(5.1) Introductory 

No discussion of experimental research in aircraft 
structures could be regarded as complete without some 
mention of fatigue. Little is known about the funda- 
mental causes of fatigue failure from the metallurgical 
standpoint. The structural specialist, therefore, has 
to rely mainly upon the experimental approach to meet 
the needs of the aircraft designer. Unfortunately, 
it is not possible to do much more than report progress 
at this stage in a task that is both long and difficult. 

Some 3 years ago the writer reviewed the fatigue 
situation from the aircraft structural standpoint in a 
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Fic. 26. Fatigue test on Meteor tailplane representing a wing. 


paper read before The Royal Aeronautical Society.° 
The conclusion then reached was that much research 
was required before the designer would have the scien- 
tific background to which he was entitled. Since then 
systematic experimental work has been pursued vigor- 
ously in Great Britain, concurrently with ad hoc testing 
of new designs of aircraft and components. 

Already much useful information has been obtained, 
but in the present state of knowledge the following 
remarks can only be regarded as a rather miscellaneous 
selection of tentative facts and preliminary ideas. 


(5.2) Fatigue Testing of Wings and Tail Planes 


One of the most profitable lines of fatigue investi- 
gation is the testing of wings and tail planes. These 
can be vibrated in their own natural modes in a way 
that reproduces a fair approximation to the conditions 
occurring in flight. At the present time, attention 
is being concentrated mainly on wings, since these 
present the more urgent design problem. 

A comprehensive investigation of the fatigue prop- 
erties of a wing is extremely difficult. Each test, 
which may take several days, gives information con- 
cerning only one particular loading cycle. For the 
simplest form of loading cycle, however, there are two 
degrees of variation corresponding to alternating and 
steady load, respectively. In consequence, the number 
of tests required is large. Furthermore, isolated tests 
cannot, in general, be accepted as conclusive, and more 
tests are required to allow for variations of nominally 
identical specimens. Still more tests may have to 
be contemplated moreover, to cover other variables 
such as temperature; or in order to follow up any partic- 
ular line of inquiry that the original results might 
suggest. 

The effort required for a complete survey is thus 
considerable and probably without parallel in any 
other field of structural research. The greatest diffi- 
culty of all, however, is to obtain the necessary number 
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of specimens, which are costly and need to be as nearly 
identical as possible. 


About 2 years ago, it was decided to test at least one 
type of structure exhaustively. To this end, some 50 
‘““Meteor”’ tail planes were diverted from the production 
line, with special arrangements to prevent modifi- 
cations being introduced. These tail planes were 
then treated as reduced-scale models of wings and were 
loaded accordingly as appropriate to civil and trans- 
port aircraft. In this way, not only was an adequate 
supply of specimens ensured, but the tests were 
rendered comparatively easy to conduct. The tests, 
of course, have no connection with the Meteor aircraft 
itself. 


Fig. 26 shows the set-up for the fatigue tests. The 
‘‘wing’’ system is vibrated close to its natural frequency 
by an out-of-balance rotating mass. Mean or steady 
load is applied by suspended weights carried on a low- 
frequency suspension system, the wings being, of course, 
inverted. Provision is also made for testing at low 
temperatures. In this case the specimen is enclosed 
in a heat-insulating box and cooled by air blown over 
solid carbon-dioxide. 


So far some 20 such tests have been made. The fa- 
tigue failures experienced usually take the form of frac- 
ture of the front or rear spars near the root, though 
such failure is usually preceded by cracking of the 
adjacent skin. 


Some interesting conclusions have already been drawn 
from these experiments. In the first place, when every 
care is taken over the tests, the degree of scatter thri ugh 
variation of specimens is much less than has generally 
been supposed. Secondly, the effect of mean load on 
this type of structure is highly significant. Application 
of a mean load of 25 per cent of the static ultimate re- 
duces the fatigue strength in the million-cycle region 
to less than half the value for zero mean load. Finally, 
with regard to temperature effect, the fatigue strength 
increases with lowering of temperature, and in the 
million-cycles region the fatigue strength at —30°C 
is about 10 per cent greater than at normal temperature 
(15°C.). Equipment is being designed, incidentally, 
to enable much lower temperatures to be reached, 
as appropriate to flight at high altitudes. 


Similar tests are also being made on actual wings of 
large size. The number that can be tested is necessarily 
limited, which is one of the reasons for having a small- 
scale version tested exhaustively to extablish general 
trends. At present, the technique for vibrating large 
wings with low natural frequencies is imperfect and 
still undergoing development. The alternative is a 
“repeated loading’’ system by which load indirectly 
applied and removed under automatic control. With 
this method, however, tests take much longer, and 
vibration methods are normally preferred. 
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(5.3) Aircraft Joints 


The major joints on an aircraft are usually the most 
critical parts of the structure from the fatigue stand- 
point. It is usual to test joints as single units before 
they are tested as part of a complete structure such as 
a wing. The tests on joints are made in one or 
other of the recognized tension-compression fatigue 
machines. 

These joints are usually bolted, steel bolts being used 
in preference to light alloy ones because of their greater 
strength and more reliable fatigue properties. Experi- 
ence shows that joints designed without special regard 
to fatigue can usually be improved upon considerably. 
The best designs have tapered cover plates or cleats, 
with bolt diameters graded to correspond. 


Bolted joints illustrate in a remarkable way most of 
the factors adversely affecting fatigue strength. Each 
hole itself produces a ‘“‘geometric stress concentration” 
in relation to the main load carried in the member. 
The high bearing stress on the side of the hole produces, 
in addition, an ‘‘applied’ stress concentration.” 
Then there is the effect of unequal distribution of load 
between bolts placed in a line longitudinally, by which 
the end bolt in particular carries an undue share of 
load. Finally, there are present, as on any other struc- 
ture component, surface imperfections that reduce 
fatigue strength and which can be particularly signifi- 
cant inside a bolt hole. 


These effects are mentioned because it is by attention 
to them in design that most improvements are to be 
effected. In this connection there is particular in- 
terest in methods that have been tested experimentally 
with the objective of improving a joint already made 
without redesign or major alteration. One method 
is to drill a small hole just outside the critical end hole. 
This diverts the flow of stress away from the loaded 
hole and reduces the stress concentration there. 
Another method is to replace the end bolt of a row of 
steel bolts by a light alloy one, thereby reducing the 
load by allowing greater deflection. In cases where 
the holes have been drilled in the ordinary way, 
moreover, it has been found beneficial to ream the holes 
to give greater freedom from scratches, oversize bolts 
then being fitted. Still better results are obtainable 
if the insides of the holes are given a “‘rolled’’ finish, 
though this may be considered too costly. 


Such palliative measures that can be applied retro- 
spectively may or may not be appropriate to an ab 
initio design. In the design of a joint there is an 
economic problem as well as a technical one. Some 
of the ways of improving the fatigue strength are costly 
in production manpower and may call for skilled work- 
manship. It may be worth while in many cases, 
therefore, to increase the dimensions of the components 
of a joint and accept the weight penalty rather than to 
call for special treatment or skilled workmanship. 


(5.4) Miscellaneous Fatigue Phenomena 


In addition to the systematic study of components 
and complete structures, there are several miscellaneous 
effects worth experimental investigation. Ideas may 
occur to those engazed on fatigue research that can 
be tried out experimentally, and, even if only a small 
proportion yield results, the time and effort has not 
been wasted. It is proposed to describe three such 
ideas that have been pursued to some effect. 

First, there is the effect of tightening steel bolts 
that are loaded in tension. Tests that have been made 
prove conclusively that a steel bolt tightened hard 
against a fairly robust fitting has a far greater fatigue 
resistance than a slack one. The tightening of the 
bolt introduces a steady tensile load that has to be 
accepted as an adverse factor, while there is a corre- 
sponding compression load in the usually much stiffer 
fitting. Any superimposed fluctuating load, however, 
is then carried by the bolt and fitting as a single member, 
the share of the fluctuating load carried by the bolt 
being small. From the fatigue standpoint the effect 
of reduced load fluctuation greatly outweighs the effect 
of the increased steady load in the bolt. Tests have 
shown that by proper pretensioning the fatigue life 
of a steel bolt may be increased to as much as 100 times 
the life of a lightly tightened bolt. 

This phenomenon raises two important issues in air- 
craft design and maintenance. Firstly, there is a 
conflict between positive tightening and positive lock- 
ing. Generally, the requirement for tight bolts should 
take precedence over locking, and the bolts should be 
unlocked and kept tight by regular inspection. Spring 
washers and soft packing should, of course, never be 
used if it is intended to take advantage of pretensioning. 
Secondly, there is the prospect of bolt slackening at 
low temperatures through the greater contraction of 
the light-alloy fitting. Where this is likely to have 
serious effects, steel bolts with a high coefficient of ex- 
pansion may be used. The alternative of having light- 
alloy bolts is not attractive because of the better and 
more reliable fatigue properties of steel bolts. 

The second subject of experimental investigation is 
the effect of bolts or pins, that are loaded in shear, 
being made a tight fit in their holes. The material 
around each hole is then in an initial state of radical 
compression and peripheral tension. The effect is, 
broadly speaking, that the hole ceases to exist as 
a discontinuity so far as concerns the main stresses 
in the member. This benefit is again offset by the 
initial preloading, but when the pretensioning is done 
properly, a considerable increase of fatigue life can be 
obtained. The negative tolerance in fit, however, is 
critical and it is not yet known whether advantage can 
be taken of the principle in practice. 

Finally, there is the effect of tightening bolts loaded 
in shear at a joint. The general effect of friction thus 
produced is, of course, well known but has not hitherto 
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been regarded as very significant in light-alloy struc- 
tures that have a low coefficient of friction. Experi- 
ments show, however, that tightening of shear bolts 
increases the fatigue life considerably. The practical 
significance of this has not been fully assessed. It 
is clearly necessary to establish what degree of reliance 
can be placed on the bolts being tightened and then 
kept tight. The phenomenon may explain the many 
discrepancies that have occurred in past results for 
fatigue life of identical joints obtained in different 
tests, especially where the tests have been made in 
different laboratories. Where somewhat involved ex- 
planations have sometimes been found, the true one 
may well be merely the different interpretations of the 
degree of tightening required. 


(5.5) Prevention of Fatigue Failure 


The study of fatigue failures in the past, as well as 
laboratory tests, emphasizes the importance of care 
in detail design. It is believed that the majority of 
serious fatigue failures that have occurred in the past 
could have been avoided if greater attention had been 
given to detail design, and particularly in the way of 
avoiding unnecessarily high stress concentrations. 
Thus, in this matter, the draughtsman takes precedence 
of the experimental scientist. With the great care 
in detail design which is now being shown, however, 
further improvement becomes increasingly difficult 
and calls for the best efforts of the structural special- 
ist. 

It is likely that for some time to come fatigue tests 
on major joints will be necessary for civil and transport 
aircraft. The sooner these can be carried out for 
any new design, in fact, the less is the delay likely to 
be caused in completing the design. It is desirable also 
that tests on complete wings, at least, should be made 
as a final check. 

The difficulty with all these tests, however, is to 
decide on what loading conditions to apply and what 
standard to accept. It is not possible to make many 
tests, and the few that are made must be reasonably 
representative of actual operational conditions. Un- 
fortunately, there is no satisfactory rule by which be- 
havior under one loading cycle can be deduced from 
another, while a satisfactory cumulative rule by which 
the effect of mixed cycles can be deduced has not been 
found and, in a general sense, may not exist. The 
study of gust loadings has been helpful in the 
choice of a standard testing condition for design accept- 


ance. 
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For aircraft in the civil or transport class an alternating 
load of +7'/2 per cent of the design ultimate, superim- 
posed on the steady level flight loading, is considered 
to provide a good working basis. The minimum number 
of cycles to failure that can be accepted, however, js 
still debatable and in any case must depend on the work. 
ing life expected from the aircraft. It is considered, 
however, that 2,000,000 cycles at the specified loading 
is not too high a figure at which to aim. 


CONCLUSION 


This completes an account of what is intended to be 
a representative selection from experimental work pro- 
ceeding in the aircraft structural field, excluding dy- 
namic work such as vibration and flutter. . 

It is proposed to conclude by making one comment 
in retrospect and to anticipate one possible criticism, 
Always the emphasis appears to be on the study of 
actual structures rather than hypothetical ones. This 
tendency has the advantage that effort is expended 
on structures that at least one professional designer 
regards as feasible and practical. Against this it can 
be argued that research should be ahead of design in all 
respects and applied to structures the designer has not 
yet thought of. 

The preference for designer's designs, however, is a 
matter of deliberate policy, though it is not carried to ex- 
tremes and tentative ventures into new types of struc- 
ture are by no means excluded. The policy is only ten- 
able with a reasonably large and vigorous aircraft indus- 
try, and such success as has been attained is largely due 
to the close cooperation that exists between scientists 
at places like the R.A.E. and designers in the British 
aircraft industry. 
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Discussions of the Lecture 


Dr. Karl Arnstein, Vice-President, Goodyear Air- 
craft Corporation: Before I go ahead with my re- 
marks I should like you to know that, when asked to 
discuss Dr. Walker’s paper, I had good reason to be- 
lieve that, as is customary, I would be one of a group of 
discussors. I, therefore, limited my prepared notes to 
a small section of Dr. Walker’s work leaving the bulk 
of the discussion to my more ambitious colleagues only 
to find that tonight I am alone on the platform. 

Dr. Walker has made an excellent presentation of the 
manner in which the many complex problems connected 
with the physicist’s approach to aircraft structural re- 
search can be solved. 

Out of the great wealth of information he has covered, 
I propose to discuss first the subject of testing of small- 
scale models, a field in which I have taken particular 
interest in the past. 

The author’s emphasis on the importance of careful 
selection of the elastic characteristics of simplified mem- 
bers in what he calls “Simplified Small-Scale Models” 
is well placed. Properly designed simplified small- 
scale models are of great value in the study of the 
underlying theory of a specific structural problem. One 
strikingly good example may be found in the model 
tests made in this country of a main frame of a large 
rigid airship of the ‘“‘Hindenburg’’ Class. This frame 
was composed of girders of approximately equal length 
located at the periphery of the circular hull section ex- 
cept at the bottom where the frame was built into a 
deep truss. The frame was braced by radial wires be- 
tween diametrically opposite corners of the polygon. 
It was found that certain loads, particularly inwardly 
acting aerodynamic forces, would cause instability of 
the frame unless a definite amount of initial wire ten- 
sion was maintained—an important piece of informa- 
tion that could not easily be gained by the mathema- 
tician’s approach. On the other hand, I know of efforts 
where misleading results have been obtained from the 
use of simplified small-scale models because some elastic 
characteristics were erroneously considered nonessential 
and were neglected. The method of analysis with 
simplified small-scale models is only as good as the 
simplifying assumptions are, an observation that ap- 
plies equally well to simplifying assumptions made 
in a mathematical approach. 

Both approaches—through models and 
matics—have been used extensively in the solution 
of structural problems in the field of rigid airships where 
exhaustive testing of a full-size airship is not practical. 

An example of tlie development of a mathematical 
approach can be obtained from a review of the history 
of the solution for the determination of the distribution 
of shear and bending stresses in an airship hull. In the 


mathe- 


early days, three different methods of analysis involving 


simplifying assumptions were tried: One used the beam 
theory as a basis; the second used the ‘‘method of 
transverse shears’’ as a starting point; and the third 
attempted to consider both, such as the one suggested 
by Hovgaard in 1922. There were good reasons to be- 
lieve that the last approximation was the closest. How- 
ever, it was not until 1938 that an exact solution was de- 
veloped in the United States by the use of synthesis of 
type solutions involving difference equations expressed 
in trigonometric series which proved the reasonable 
accuracy of the Hovgaard approximation. 

It is interesting to note that, prior to the availability 
of this exact theory, the degree of accuracy of the 
Hovgaard assumption was completely confirmed by 
tests on reduced-size true small-scale models, which 
were scaled down in correct elastic relationship with re- 
gard to axial stiffness, radial bending stiffness, and tor- 
sional stiffness. With an apology to Dr. Walker, I wish 
to state that my sympathies still rest very strongly with 
true small-scale model tests, which model tests I prefer 
even to measurements on prototype ships. This prefer- 
ence is based on the fact that in testing the prototype 
ship many strain gages are required at any section of a 
girder to obtain an accurate picture of the total load in 
the girder, while, by the use of models, the model mem- 
bers themselves can be designed so that they, in effect, 
become large strain gages for summing up the stresses 
imposed on the girders. Another point of interest is 
that the elastic range of such a model can be made to 
exceed the scaled-down prototype yield point, and, 
therefore, stresses may exceed the yield point in some 
members without damaging the model. 

While there are many applications for which models 
can be used, their use is not recommended in the field 
of design of joints. It is our practice to test full-scale 
specimens, for the greatest difficulty in joint design is 
that some portion of the joint is apt to pick up loads for 
which it is not designed. Consequently, any simplifica- 
tion of a joint may cause misleading results. 

Reviewing the available methods of analysis, we are 
reminded of the expression ‘“‘many roads lead to Rome.” 
This is also true for engineers seeking the answer to any 
one of their many problems. It is my opinion that 
every engineer should take the road with which he is 
most familiar and on which he has developed his greatest 
skill, whether it be the approach by simplified small- 
scale models, by true complete small-scale models, or by 
one of the many approaches of mathematics. We are 
greatly obligated to Dr. Walker for pointing out to us 
a road not very widely known or used which offers a new 
and speedy solution to certain problems. 

Another remark of Dr. Walker’s which interested me 
was his opinion that ‘‘the most efficient designs are most 
likely to be obtained when the designer deliberately 
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takes a few chances and relies on strength tests to de- 
tect any weakness’ and that ‘‘structures that failed 
prematurely in the first test were ultimately among the 
most efficient of them all.” 

Procedures of this kind were particularly indicated in 
the early days of aircraft construction, where great care 
was necessary in order to obtain a sufficiently light 
structure, since only inefficient materials were avail- 
able. 

When Count Zeppelin submitted his airship design to 
the Prussian Academy of Science about 60 years ago, he 
was told that an airship could not lift its own weight if 
it were to be designed strong enough to withstand any 
breeze of consequence. The Count succeeded in his 
work only by taking, nota few, but a great many chances. 
Consequently, he had a right to believe that the useful 
lift benefited from the amount of our scientific ig- 
norance and that the useful lift deteriorated with the 
increased application of scientifically not fully digested 
information. He considered an airship overweight 
when it returned from its first trial flight without struc- 
tural failure. 

It is still true today that the greatest advances are 
made by lessons learned from failure. 

In conclusion, I wish to emphasize that Dr. Walker's 
paper presents some of the important problems con- 
fronting the practical designer and gives the status of 
our knowledge regarding their solutions in an interesting 
way without the use of obscuring theoretical dis- 
cussions. The author’s beliefs and experience deserve 
careful study by every practical aircraft designer. 


A. B. Callender, Supervisor—Structures Test Labo- 
ratory, Glenn L. Martin Company: Dr. Walker's 
paper is most successful in pointing out basic aspects 
of aircraft structural testing, this in itself being no 
mean task with respect to a field of engineering which 
deals in the main, with problems involving one speci- 
ally tailored job after another. In addition, he is 
deserving of a vote of thanks from all of the structural 
test engineers in the airplane industry for presenting 
their story so expertly. 
craft work merits more publicity than it is usually 


This important branch of air- 


accorded, among other reasons, because of the com- 
petition it faces from newer fields of activity for interest 
and attention from young engineering graduates. 

The presentation of test results on a large number of 
wings is of great interest. We find that we have more 
or less subscribed in recent years to the design philos- 
ophy expressed by Dr. Walker. The same trends 
from old to new designs are evident as indicated by the 
lable 1) carried out at The Glenn L. 
Martin Company over a period of 12 years. The wing 


tests (shown in 


specimens which failed to meet the design requirement 
the first time were able to be modified and retested 
successfully without an excessive amount of trouble. 
However, it is apparent that a point of diminishing re- 


TABLE | 
Wing Specimens 
Percent of Design Re- Percent of Design Re. 
quirement Attained in quirement Attained jn 
Original Static Test Retest After Modifica. 
O07 ) 


Wing Year (% tion (%) 

y. 1939 138 

B 1942 123 

¥ 1943 107 

D 1945 94 108 (no failure) 
E 1947 100 

F 1947 94 108 (no failure) 
G 1951 90 Test not complete 


turns is being approached wherein design efficiencies 
for weight- or cost-savings have to be balanced against 
modification difficulties, retests, and possible produc- 
tion schedule delays. 

The section of the paper dealing with aircraft fatigue 
presents an excellent summary of the problem from the 
structural testing standpoint. We hope that Dr. 
Walker will be able to publish more information on the 
tests of the ‘“‘Meteor’’ tail planes when it becomes avail- 
able. It would be of interest to those working in this 
field if the lecturer would care to outline briefly the 
background of his selection of +7!/2 per cent design 
ultimate alternating load and 2,000,000 cycles as a 
reasonable compromise for a service-life test. An 
alternating load of 7'/2 per cent design ultimate on a 
4'/.¢ airplane is approximately '/;g. If, in approach- 
ing this admittedly tenuous proposition, the frequently 
mentioned Cumulative Damage theory is accepted as a 
criterion for estimating the load level producing the 
most damage, the magnitude generally seems to turn 
out somewhat higher than Dr. Walker’s figure, perhaps 
on the order of °/sg or */,g. 

It is becoming evident that in the not too distant fu- 
ture some structural tests on components may have 
to be made under elevated temperature conditions. 
The problems presented by this development will open 
up a broad new field of activity in experimental air- 
craft structural research wherein different technical 
knowledge and much ingenuity will have to be exer- 
cised by the test engineers. Recent investigations in- 
dicate that, for representation of transient tempera- 
ture effects at least, the best lines of approach may lie 
in the directions of induction heating or radiant heating. 
One aspect of this development which is quite ap- 
parent is that the designers and interested agencies will 
have to accept much more extensive test compromises 
with respect to load simulation in combination with 
temperature simulation than they have been accus- 
tomed to when the matter is one of load simulation 
alone. 


L. M. Hitchcock, Staff Engineer—Structures, Boeing 
Airplane Company: It is with interest that we learn 
of the small-scale model testing being. accomplished in 


Great Britain. A parallel development has recently 
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been in progress at the Boeing Airplane Company. 
It was begun for the same reasons that prompted the 
Royal Aeronautical Establishment to use these models 
as vehicles of structural research. It is believed that the 
most fertile fields for this method of research lie chiefly 
in determining load distributions and material effec- 
tiveness in complex structures not lending themselves 
to ready analysis and in making comparative efficiency 
checks of various types of structural design. 

With reference to strength testing of complete air- 
craft, Boeing policy in the past several years has been 
to static-test complete airplane structures of each new 
normally the first airplane from the 
production line. Static tests of prototype or experi- 
mental models are not, in general, required except for 
miscellaneous panel and fitting tests and proof tests 
of control systems. Dr. Walker has indicated that it 
is common practice in the case of large aircraft to test 
wing and fuselage separately. We found in planning 
static tests for the B-47 that it was difficult and im- 
practical to separate the wing and body tests for a 


production type 


swept-wing airplane because the critical condition for 
the inner wing was also critical for the fuselage center 
section. Therefore, the B-47 wing and body were 
tested concurrently for the most critical design condi- 
tion in order to ensure proper distribution of the stresses 
to each part of the structure including any interaction 
effects. 

In our experiences of static testing airplanes, we have 
made extensive use of electric strain gages for both con- 
trol of the test and check of the designer’s assumptions 
In the test of the B-47 airplane, 

These gages were read at each 


of load distribution. 
1,500 gages were used. 
load increment, and the gages at critical locations were 
watched continuously to help guard against major fail- 
ure. As many as 400 gages were read and plotted 
during each load increment for a major test. These 
gages proved invaluable in explaining the causes of 
any failure that occurred and in selecting the best way 
of increasing the strength of the airplane if desired. 

Our history of initial failures in static testing is in 
general agreement with Dr. Walker’s chart, although 
on wing primary structure our scatter has been in the 
order of 80 to 120 per cent of required ultimate strength. 
We have lost or destroyed two wings because of cata- 
strophic failure below 100 per cent of design load: One 
wing failed at 92 per cent and the other failed at 95 
We required a complete new wing for the 
However, 


per cent. 
test in which failure occurred at 92 per cent. 
the test of the other airplane wing to 95 per cent was 
considered proof of 100 per cent after determining the 
reason for, and the location of, failure and incorporat- 
ing the structural changes considered necessary for 100 
per cent strength. 

We have been concerned with the trend in recent years 
to design airplanes with an absolute minimum in struc- 
tural weight and also with the use of less ductile mate- 
rials, which have a tendency to produce more scatter 


on the side below 100 per cent. Also, if several 100 
per cent tests are to be performed on a particular com- 
ponent, the problem of local distortions occurring in 
one test to 100 per cent may influence the stress dis- 
tribution in a test for another condition to the extent 
of causing premature failure below 100 per cent. In 
addition, any local failures in the range of 95 per cent 
to 100 per cent are more likely to induce complete 
failure of the structure, since there is little reserve in 
immediately adjacent structure to absorb any resulting 
redistribution of load. We have therefore considered 
testing structures to 95 per cent instead of 100 per cent 
for all tests on each structural component except for 
the final test to destruction. This manner of testing 
relies upon the destruction test substantiating our 
methods of stress analysis. We then have reasonable 
assurance that all other conditions tested to 95 per 
cent can be evaluated as to their ultimate strength. 
This approach, of course, is primarily directed toward 
saving of calendar time in the test program and reduc- 
ing costs of replacement of structures by eliminating 
those few failures that might occur between 95 and 100 
per cent. Structural changes required for airplanes in 
production must be determined as early as possible so 
that the least number of production airplanes will be 
affected. 

With regard to the philosophy of deliberately design- 
ing structures by taking calculated risks, it has not 
been Boeing’s policy to follow this course. We have 
enough trouble with premature failures due to im- 
proper detail design and errors without adding the 
risk of further delays and structural changes that 
might be required because of such a policy. Depending 
upon the quantity of production airplanes being pro- 
duced, this can have a far-reaching effect upon rework 
of airplanes in the factory or possibly in the field. I 
might add that not too long after static tests are com- 
pleted, the problem of airplane growth usually faces us. 
To my mind, it is preferable to design to scratch on the 
basis of minimum guaranteed properties and hope to 
demonstrate by static tests that the airplane has excess 
strength that can be utilized in permitting airplane 
growth rather than end up with a structure that might 
reduce permissible gross weights below the design 
value. In this regard we attempt to learn as much as 
possible about the strength of various sections or com- 
ponents of each major structure by progressive destruc- 
tion tests. For example, the ultimate strength of the 
B-47 wing was obtained at three locations along the 
span. 

Dr. Walker’s discussion of the universal test frame 
is convincing for the case of a central test facility that 
is used to conduct tests on many types of airplanes. 
The initial high cost of installation is justified in view 
of the saving in setup time and man-hours for a num- 
ber of static test airplanes and the probable high utiliza- 
tion of the test frame. However, if the tests are con- 
ducted at the manufacturer’s plant, the initial expendi- 
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Fic. A. Pilot’s compartment after pressure test XB-29. 


ture for a static test facility may be of a considerably 
lower magnitude, since substantially less flexibility is 
required. 

In the considerations of testing of large aircraft at a 
central testing agency in contrast to testing at a manu- 
facturer’s location, we must consider the cost involved 
in preparation of the static test airplane for flight to 
the central testing agency, the delays in the test pro- 
gram caused by preparation of the airplane for flight, 
and the additional delays and costs in the test program 
if extensive repairs or replacement of major components 
It is 
aiso possible that premature failures in tests of large 


are required away from the contractor’s plant. 


aircraft will require replacement of major components 
that cannot be shipped from distant points for con- 
tinuation of the test program. 

In addition to the above items that affect the timing 
of changes required on production airplanes, we also 
have to consider the advantages gained by the manu- 
witnessing static 


facturer’s structural engineers in 


tests conducted at the contractor’s plant. A consider- 
able amount of this experience is lost if the tests are 
performed at a central test facility, because a large 
number of design engineers cannot be spared to wit- 
ness these tests. 

Testing of pressure cabins can indeed be disastrous. 
In the pressure test of the B-29 airplane, the complete 
nose section failed, as shown in Fig. A. This has been 
the exception rather than the rule because of our basic 
design policy for pressure structures. Because of 
fatigue considerations, it has been our practice to 
provide strength in excess of that required by the 
Thus, in the 
test program this fatigue reserve is reasonable insur- 
Protection 


usual factors of safety on the static load. 


ance against primary structural failure. 
from failure of doors and windows is provided, since 
these are the sources of most of our troubles in pres- 
sure sections. 

The difficulty attached to a program submerging a 
full-scale pressure cabin of a large airplane into water 
and pressure testing is so great without special facili- 
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ties that at present no such test procedure is contem. 
plated at Boeing. 

The counting accelerometer discussed by Dr. Walker 
is an interesting development with several ingenioys 
features incorporated to reduce (using Dr. Walker's 
words) “‘the danger of producing too much informa. 
This raises a question fundamental to experi. 
Should we develop instru. 


tion.”’ 
mental work of all kinds. 
mentation purposely to exclude information that ap. 
pears at the moment to be superfluous? The N.A.C.A, 
V-g recorder is a case in point. When this instrument 
was originally designed, the need for statistical study 
of accelerations significantly less than the extreme 
values was not generally realized. Accordingly, a de- 
vice was conceived which, intentionally or not, ob- 
scured all data in régimes not then of immediate inter- 
est. We now look back on the data recorded in those 
thousands of flight-hours and wish we had been smart 
enough to anticipate our present intense interest in the 
statistical characteristics of the lesser accelerations, 
In fact, it is the limiting characteristics of the V-g re- 
corder which have made necessary the reinstrumenta- 
tion of further aircraft with Dr. Walker’s acceleration 
counter. Even this instrument will give little if any 
direct information on airplane gross weight, altitude, 
and air speed at the instant of encountering the ac- 
celeration, gust gradient distance, and other factors of 
major importance to various aspects of the pressing 
structural design problems. I should like to suggest a 
diametrically opposite approach to problems requiring 
experimental study. 

It would seem more desirable in experimental study 
to design instrumentation to record as completely as 
practical all factors of likely interest with regard to the 
particular subject. Having thus obtained a permanent 
record that may 
some present pressing problem, it becomes a matter of 


include information extraneous to 
applying ingenuity to the development of data reduc- 
tion techniques. Later, when other aspects of the sub- 
ject are of interest, the same experimental records can 
be reanalyzed to provide additional useful information 
without the delays incident to creating new specialized 
instrumentation, obtaining necessary program coordi- 
nation, and then waiting the time necessary to accumu- 
late statistically significant data samples. 

In addition, I should like to venture the opinion that 
any high-frequency accelerations of magnitudes near 
the same order as the lower frequencies for which the 
instrument was designed would be of extreme signifi- 
cance to the fatigue problem if they were transmitted 
through the structure. While the amplitude of mo- 
tion required to produce a given acceleration is inversely 
proportional to the square of the frequency, the shear 
and bending stresses for a given amplitude are higher 
for the higher modes because of the associated more 
complex mode shapes. If these higher frequency ac- 
celerations are vibrations in the usual sense rather than 
the structural response to some transient external ex- 
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citations, they may be of even greater significance be- 
cause of their tendency to be ever present. 

The tests of the Meteor tail planes are of consider- 
able interest, and at such time as the information is 
made available it should contribute considerably to 
basic data on fatigue characteristics of structures. It 
is noted that the fatigue tests have been conducted at 
415° and —30°C. We wonder if the Royal Aeronauti- 
cal Establishment plans to conduct tests on major 
components at elevated temperatures. 

We agree with the comment that detail design of 
structural joints is a most important factor from a 
fatigue standpoint and that improvement in their 
design is the most expeditious method of obtaining in- 
creased fatigue life. Dr. Walker does not mention the 
comparative fatigue properties of various high-strength 
aluminum alloys. Our own tests, as well as other 
data, show that increased fatigue life can also be ob- 
tained by choice of material, such as usage of 24ST 
instead of 75ST. This is considered particularly im- 
portant when tension material of a structural com- 
ponent is concentrated in one or two members. 

The preloading of tension loaded bolts by the usage 
of nut-tightening torques is universally accepted as a 
method of reducing the alternating stress in the bolt. 
The proposed usage of high coefficient of expansion 
steel bolts to reduce the effect of low temperature 
slackening of the bolt tension is novel. We would be 
interested in knowing if special steel bolts have been 
evaluated for this application. I would presume that 
Dr. Walker is familiar with our practice of using bolts 
with threads rolled after heat-treatment with the at- 
tendant increase in fatigue life. 

The amount of interference used in the interference fit 
bolted joint tests would be of interest. The tests con- 
ducted by Battelle Memorial Institute for the N.A.C.A. 
as reported in N.A.C.A. Technical Note No. 1030, using 
bolt clearances from —0.0005 to +0.050 in., showed no 
clearly defined effect on the fatigue strength of bolted 
lap joints from bolt clearance. It would be presumed 
that the lecturer’s data were based on tests of consider- 
ably higher interferences. We would be concerned as 
to the possibility of stress corrosion with the use of 
high interference fits. However, interference fits pos- 
sibly tend to alleviate the problem of fretting corrosion 
which has occurred on bolted connections with clear- 
ance holes. With regard to bolted connections, has any 
work been done on relative fatigue life of bolted versus 
redux connections? 

Dr. Walker has indicated that an alternating load of 
7'/» per cent of the design ultimate provides a good 
basis for the study of wing fatigue life. Our studies also 
show that this general region of alternating load for a 
four-factor airplane causes the major portion of fatigue 
damage in service. However, this is based on cumula- 
tive damage theory and published data on gust fre- 
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N. J. Hoff, Head, Department of Aeronautical Engi- 
neering and Applied Mechanics, Polytechnic Institute of 
Brovklyn: Perhaps the most interesting idea brought 
up by Dr. Walker is the distinction he makes between 
the theoretical and the experimental approaches to re- 
search. I should like to quote his statement that 
“the experimental approach is more scientific than 
some might suppose, and it does not lead anywhere 
unless it is closely related to theory.’’ This may sound 
strange to those experimentalists who still believe that 
the only purpose of testing is the measurement of the 
failing load of a structure. Dr. Walker’s lecture illus- 
trates clearly the great amount of theoretical knowl- 
edge that is required nowadays in order to plan tests, 
design test rigs, develop measuring instruments, and 
interpret the results obtained. It is safe to say that in 
reality there is little difference between theory and 
experiment. The theoretician who disregards experi- 
ment engages in mathematical deductions from first 
principles established by experimentalists decades or 
centuries ago and does not admit the possibility that 
present-day experimenters may add to or modify these 
principles. The experimental man who disregards the- 
ory simply puts his own personal experience above the 
fund of knowledge collected and ordered into a con- 
sistent system by great men who lived before him. 
Sound progress is possible only if theory and experi- 
ment merge as they do in this Wright Brothers Lecture, 
even though the emphasis in it is on the experimental 
approach as signified by the title. 

One of Dr. Walker's great contributions to the de- 
velopment of our knowledge of structural behavior is 
his analysis of force control in testing. The design of 
improved strainers followed from this work, which in 
turn permitted the construction of a test frame of the 
efficiency and size of the ‘‘Cathedral."’ Not far from 
Farnborough in the south of England, tasks of a simi- 
lar magnitude were accomplished by churchmen 
hundreds of years ago when they built cathedrals of 
unrivalled beauty dedicated to purposes different from 
those of Dr. Walker’s Cathedral. 

Another favorite subject of the lecturer is fatigue 
His words on this topic inevitably recall that un- 
paralleled source book of information on fatigue written 
by Nevil Shute entitled No Highway. At the time the 
book was written, there was only one chief of struc- 
tures at Farnborough; his name was not Scott but 
Walker. Dr. Scott's efforts were toward 
merging the purposes of scientific research with the 
practical requirements of airplane construction and 
operation. He had to fight the group of research men 
who were characterized in the book as follows: 

“There comes a time when the research worker, dis- 
appointed in promotion and secure in his old age if he 
avoids blotting his copybook, becomes detached from 
all reality. He tends to lose interest in the practical 
application of his work to the design of aeroplanes, and 
turns more and more to the ethereal realms of mathe- 


directed 
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matical theory; as bodily weakness gradually puts an 
end to physical adventure he turns readily to the adven- 
ture of the mind, to the purest realms of thought where 
in the nature of things no unpleasant consequences can 
follow if he makes a mistake.”’ 


Judging from this Fifteenth Wright Brothers Lec- 
ture, Dr. Walker has been singularly successful in 
minimizing this kind of thinking in his group and has 
obtained results that are theoretically interesting, as 
well as useful in practice. 


Paul Kuhn, Langley Aeronautical Laboratory, 
N.A.C.A.: Dr. Walker has presented a survey of 
structural research methods and results which is bound 
to be of great interest to anybody concerned with such 
work. I should like to make some comments on those 
phases of the lecture that fall within the sphere of 
N.A.C.A. activities. 

First of all, I want to congratulate Dr. Walker and 
his staff on the interesting and useful work with plastic 
models. Aircraft structural designers are confronted 
nowadays with many difficult problems of which delta 
wings are one example. There is little theory for such 
structures, even for the sample case of a pure shell, 
let alone the much more complex case where the shell 
structure is heavily modified by adding spars and heavy 
ribs to take care of a landing gear, for instance. For 
such structures, plastic models can quickly provide 
answers that are much more reliable than those ob- 
tained by means of arbitrary assumptions. There is at 
least one case on record where a bad assumption led 
to a 30 per cent error, with serious consequences. 
There are some research men and some designers who 
look askance at plastic models. I hope sincerely that 
Dr. Walker's lecture will be instrumental in securing 
for such model tests the attention they deserve. 

Dr. Walker mentioned a comparison between what 
he called aerodynamic ribs and structural ribs in swept 
wings. Similar comparative tests have also been made 
by the N.A.C.A. The differences were found to be as 
small as those mentioned in the lecture—of the order 
of 5 to 7 per cent—but in favor of the aerodynamic ribs. 
As far as bending deflections are concerned, this result 
agrees well with that estimated by elementary theory if 
one takes into account what civil engineers call par- 
ticipation stresses—that is to say, the contribution 
made by the inclined ribs to the moment of inertia of 
the beam. In the light of this experience, the results 
quoted by Dr. Walker are a little surprising, and it 
may be prudent to interpret differences of less than 10 
per cent with caution until the test techniques em- 
ployed on plastic models are more fully developed. 

Somewhat surprising, also, but for a different reason, 
is the first test mentioned—a wing with overhung 
engines. It would seem that a simplified model test 
was hardly necessary, particularly after a full-scale 
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test had been made, because a theoy of bulkhead flex. 
ibility was established some time ago, and one or two 
reports published in England dealt exclusively with 
this problem. Is this, perhaps, another case of g 
prophet being without honor in his own country, or was 
the specific source of trouble one not covered by the 
existing theory? 

Dr. Walker mentioned as a general conclusion from 
full-scale tests that the most common cause of pre- 
mature failure is a stress-diffusion problem. I should 
like to underscore this remark heavily and hope that it 
receives full attention. In the past, there was some 
justification for treating such problems rather lightly 
in design. However, hidden factors of safety are being 
eliminated more and more, and the materials are be- 
coming less ductile. As a result, stress diffusion prob- 
lems have become extremely important indeed in 
static strength; they are, of course, also a prime source 
of trouble in fatigue. 

The counting accelerometer described in the lecture 
embodies an ingenious approach to the difficult prob- 
lem of eliminating extraneous vibration. I wonder if 
Dr. Walker would care to say a few words on the in- 
tended use of this instrument. It will evidently give a 
record of the load history of the airplane on which it is 
installed; this record might be used to estimate the 
remaining fatigue life of the structure. However, the 
application of the records to the design of new air- 
planes would seem to be rather uncertain, because the 
loads experienced are directly proportional to the speed 
of the airplane, and information on the speed is not ob- 
tained. 

The fatigue tests on tail planes are also interesting, 
but they prompt another question. Our predictions 
of fatigue life are highly uncertain at present. It ap- 
pears that some fundamental work has been done in 
England which promises to eliminate this vexatious 
problem. I wonder when Dr. Walker expects to be 
able to persuade the originator of this theory to pub- 
lish it. I am referring to the theory of fatigue de- 
veloped by that eminent staff member of the R.A.E., 
well known on both sides of the Atlantic, Mr. Honey. 


J. F. McBrearty, Structures Division Engineer, 
Lockheed Aircraft Corporation: 


indeed when an engineering profession is favored with 


It is a rare occasion 


a thorough and lucid inventory of the major accom- 
plishments and problems in that field such as Dr. 
Walker has assembled in this Wright Brothers Lecture 
for 1951. The aircraft structural engineering profes- 
sion is indebted to him for the insight and practical 
touch with which he has seasoned this paper. In addi- 
tion to important concepts of engineering philosophy 
and design policy, he has interlaced this lecture with 
most of the important down-to-earth facts discovered 
in the course of structural research and testing in 
recent years. 
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One of the most important, and perhaps controver- 
sial, points discussed by Dr. Walker is the philosophy 
of designing a structure ‘‘on the light side’ and then 
developing suitable reinforcement in static tests, and 
| would like to comment on this one aspect of Dr. 
Walker’s subject. It is true that this approach has be- 
come more widely accepted in recent years due largely 
to increased demand for weight economy and to more 
and more complicated structures, and yet I feel that 
this approach can be carried too far for our own good. 

In the first place it must be remembered that the 
design condition and the corresponding test does not 
constitute an end in itself but is merely the means to 
an end—that end is to produce a structure that will be 
adequate in service. It is not always true that the de- 
sign conditions are just the necessary and sufficient 
criteria to result in a satisfactory service structure, 
and all too often the structural elements that show dis- 
tress in static tests are not at all the elements that 
prove to be the Achilles’ Heel in service. Frequently, 
100 per cent of some design condition in a static test 
is viewed as a ‘‘magic number” and exact achievement 
of it assumed synonymous with the structure living 
happily ever after. I wish this was true—certainly the 
life of the structural designer would be a happier one 
if it were—but it seldom is. Design conditions are 
rarely so rational and realistic and often are completely 
arbitrary, and in any case an arbitrary safety factor is 
superimposed on any condition, however realistic. 
Consequently, I feel that a good-stress analysis fre- 
quently is the better index of structural adequacy and 
would be equally fruitful in minimizing weight if the 
stress analyst were permitted to assume the very things 
that are taken for granted in a static destruction test, 
such as above average material properties, favorable 
load distributions, and so forth. 

It has been our experience in many cases that the 
stress analysis, moreover, is a better indicator of 
troublesome details than is a static test, which often 
glosses over stress raisers and the like by plastic relief. 
Thus, it should be made clear that the static test of the 
structure is a supplement to, and not a substitute for, 
a thorough and careful stress analysis and does not 
relieve the structural designer of the necessity for ap- 
plying judgment to preclude service 
troubles. 

Perhaps the most practical consideration in this ap- 
proach to structural development is the timing. If 
the structural test isn’t early enough in the schedule, 
then it produces none of the advantages yet all of the 
disadvantages of this procedure. The attractiveness of 
the development by test scheme fades if the test is 
going on after all the production tooling is built and 
long lines of production articles are in the jigs or (as 
sometimes the case) even in the skies. Thus, while this 
approach may be economical and practical for the de- 
sign of a prototype or small quantity contract, it is 
doubtful that it is the wise course to pursue when a 


experienced 


RESEARCH 167 


design is scheduled into early production. In any event 
great emphasis should be given to providing the struc- 
tural! test laboratory with facilities for the rapid con- 
struction of test structuresand the speedy conduct of the 
tests. Dr. Walker has described some of the advances 
made in this latter field, and they are certainly to be 
encouraged. No less important, however, is the prob- 
lem of getting the test data and report into the hands 
of the designer who uses them, and, consequently, there 
should be more and more mechanical aids to this phase 
of the work. Automatic strain recording devices, me- 
chanical and electronic computing machinery, and 
dictaphone equipment should be employed to mini- 
mize the interval that the research results spend be- 
tween production and delivery to the consumer—the 
designer. 

It is believed that structural development by test is 
here to stay and that it enables the designer to exploit 
at the outset much of the airplane “‘stretch’’ that here- 
tofore was a slow process spread over the years as the 
airplane demonstrated in service its further capabilities. 
The designer, however, must be discreet in his employ- 
ment of this powerful method and not forget the pitfalls 
mentioned here. 

Once again, I wish to commend Dr. Walker for pre- 
senting us with a most informative lecture for which 
all aircraft structural engineers should indeed be 


grateful. 


Walter Ramberg, Chief, Mechanical Division, Na- 
tional Bureau of Standards: Many of us in the field of 
structural testing will envy Doctor Walker his impres- 
sive ‘“‘cathedral’’ test frame for applying vertical loads 
to full size wings and fuselages of aircraft and will wish 
him all success in obtaining for his laboratory the still 
more ambitious test frame described by the models in 
Figs. 14 and 15. 

Facilities such as these give the Royal Aircraft 
Establishment an exceptional opportunity to collect 
data on the modes of failure under static load of aircraft 
of many different types. As a result, Doctor Walker 
and his associates should be in an excellent position to 
answer a question that has bothered all those who have 
had to rely on analyses and tests of relatively small 
components in estimating the strength of the full size 
structure. The question is simply to what extent can 
one rely on careful design and on tests of small com- 
ponents? Or, to put it quantitatively, in what per- 
centage of the tests of complete structures was the 
failure due to a “general instability’’ in which the 
entire structure was involved rather than a local weak- 
ness that could and should have been caught before the 
test of the complete structure? Figs. 6 and 7 of the 
paper indicate graphically the improvement that re- 
sults from tests of complete wings and fuselages, but 
the failures are not separated into those that might 
have been prevented by more careful design and con- 
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struction and those that could not have been predicted. 
Doctor Walker may recall that the aircraft structures 
people became seriously worried about ‘“‘general insta- 
bility’’ following the tests of reinforced cylindrical shells 
conducted by von Karman, Dunn, and Sechler at the 
California Institute of Technology about 10 years ago 
(see for instance N.A.C.A. Technical Notes 905-911). 
These showed convincingly that it was sometimes 
unsafe to neglect the complicated interactions between 
different portions of a stressed-skin aircraft structure in 
estimating its margin of safety against instability. 

This discusser would like to compliment Doctor 
Walker upon the development by his colleague, James 
Taylor, of a counting accelerometer that appears to 
function without a single electronic component in it. 
This should be a challenge to some of our American 
test engineers who have become prone to capitulate to 
the electronic engineers before fully exercising their 
ingenuity in devising a purely mechanical instru:nent 
and making it work. 


E. E. Sechler, Professor of Aeronautics, California 
Institute of Technology: In these days of new and 
elaborate computational methods and machines, it is 
refreshing to hear so eminent a scientist as Dr. Walker 
advocate not only the continued use, but the expansion 
of, experimental techniques for the solution of some of 
the complex problems of aircraft design. It is par- 
ticularly fitting that the experimental method should 
be re-emphasized in a Wright Brothers lecture, since it 
was essentially through experimentation that the 
Wright Brothers were finally able to build and fly their 
history making airplane. 

So many of our young engineers grow up with the 
belief that if, by some magical means, a problem can be 
fed to an analog or digital computer, the answer that 
comes forth in a matter of seconds will be the solution 
to end all of their difficulties. Unfortunately, this is 
not true, primarily because the complex computing 
machine, with all of its capabilities, has no power to 
modify assumptions based upon actual experience as to 
the manner in which a structure will react under vari- 
ous loading and restraint conditions. 

Experimentation such as that discussed by Dr. 
Walker is, in many instances, the only way in which 
proper assumptions as to the behavior of complex 
structures can be made. Even when full-scale assem- 
blies of metal structures are tested, the designer is fre- 
quently limited to visual observations on the outer 
portions of the assembly and must rely upon remote 
reading instruments to give limited information as to 
the behavior of internal elements. The use of optically 
transparent material, such as Xylonite, for models often 
reveals interaction effects that would not otherwise 
have been taken into consideration. 

I do not wish to give the impression that the use of 
computers should be eliminated or even decreased. 
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What is required is the use of all available methods ty 
obtain the best possible design for a given aircraft 
structure. This will entail the integrated use of the 
highest level theoretical knowledge, the best computa- 
tional techniques, and the fullest possible application oj 
experimental methods. Dr. Walker's paper shows that 
the experimental method can well take its proper place 
in this trinity of tools for design analysis. 


C. R. Strang, Strength Engineer, Douglas Aircraft 
Company, Inc.: Someone has said that nobody believes 
a theoretical analysis except the analyst who prepared 
it, while everybody believes a test except the test en- 
gineer who prepared it. This too true statement aptly 
conveys the idea that tests are convincing regardless, 
The admirable thing in Dr. Walker's philosophy is pre- 
cisely that it never occurs to him to test for the test’s own 
sake but that there isalwaysa purpose. Dr. Walker has 
effectively shown that the experimental method is nowas 
important as it ever was and that the situation will 
probably continue that way. At least in the foreseeable 
future, there are certain to be radical changes in struc- 
tural configuration from one model to the next, as well 
as a constant pressure on the designer to hold structural 
weight to the minimum. This means no letup in the 
need for tests to prove theoretical strength analyses. 

The test methods and testing philosophy of the 
manufacturer of transport air frames are considerably 
different from those of an agency such as the Royal 
Aeronautical Establishment. For example, where the 
latter is able to concentrate on destruction tests of 
complete airplanes, the manufacturer is forced by con- 
siderations of timing and economy to do most of the de- 
struction testing on components of the structure (see 
Figs. B-F) during the design stage of the airplane, fol- 
lowed by a limit load test on the complete airplane. This 
procedure is not all that might be desired, but it has re- 
sulted in some highly successful models. In the follow- 
ing discussion of our test methods it will be seen that 
there are other differences in technique, which should 
not necessarily be interpreted as adverse criticism of 
either one. 

Our custom is to use hydraulic tension cylinders as 
much as possible to apply and measure the test loads. 





Section of the DC-4 wing set up for repeated load testing. 
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AIRCRAFT STRUCTUR 


The jacks are manifolded together in the minimum 
number of groups, each at a common pressure, which 
will result in the test specimen being stable. A con- 
ventional system of levers is used to make this possible. 
The system is such that, if one part of the structure fails, 
the remaining loads will be relieved without, it is hoped, 
any excessive dynamic forces. Incidentally, the Austra- 
lian Council for Scientific and Industrial Research has 
published a simple and efficient method of arranging link- 
age systems to obtain such common pressure groups of 
jacks. * The all-hydraulic system of loading has the 
advantages of simplicity of installation and of operation, 
accuracy, and sensitive control. It stores practically 
no energy within itself during loading (providing no 
accumulators are incorporated), with the result 
that damage to the structure after failure is mini 
mized. 

The problem of handling large deflections is becoming 
more and more acute as wing spans increase, while rela- 
tive wing thicknesses simultaneously decrease; In this 
case it seems that a combination of a screw jack with a 
moderate-size hydraulic jack would be more practical 
than a long-stroke jack by itself. For this purpose the 
mechanical and hydraulic jacks would probably be con- 
nected in series, rather than arranged as in the R.A.E. 
lever system. 

Air-pressure tests of windows, astrodomes, and other 
glass or Plexiglas structures almost always result in 
violent failures with all evidence of the locality of the 
initial failure destroyed. This is due (in addition to 
the causes discussed by Mr. Walker) to the brittle 
nature of the structural material. Such structures, 
therefore, are ordinarily tested with water. Small pres- 
sure vessels of metal are most conveniently and safely 
tested in this manner also. 

Recently, we tried out a variation of this technique 
in which we ran a vacuum (instead of positive pressure ) 
test on a vessel. The specimen was a barrel-like en- 
closure for an auxiliary power plant. The water was 
drawn from the test specimen by means of a small 
rotary pump and was discharged into the tank in which 
the test specimen was submerged. A small volume of 
air, unexpectedly trapped in the specimen, expanded to 
a large volume when the pressure was reduced. 
Also, some of the air in solution in the water ‘‘boiled”’ 
out to further increase the volume. The net result was 
amore general type of failure than would have occurred 
if the air had been completely removed at the start of 
the test (Fig. C). 

The failure of a fuselage under negative pressure is 
more violent than would be expected, as we learned at 
the end of a test program on a 15-ft. section of DC-6 fuse- 
lage from which the air was evacuated until the struc- 
ture failed. The fesults of the “implosion” that oc- 
curred at —3.4 Ibs. per sq.in. are shown in Figs. D 


and E. 


* Report SM52 by W. W. Johnstone, “Strength Testing of 


Aircraft Wings—A Review of Current Methods.” 
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Auxiliary power-plant enclosure after failure due to 
external pressure. 





Fic. D. 


Test failure of a fuselage section due to external pressure. 





Fic. E. 


Test failure of a fuselage section due to external pressure. 
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Fic. F. Pressure test setup for the forward section of the DC-6 
fuselage. 


For pressure testing of major fuselage sections we 
have favored air, rather than water, as the pressure 
medium for the following reasons: 

(1) Using air, the structure is easily and quickly 
accessible for observation and for inspection and 
repair of damage during test. 

(2) The structure is not so subject to corrosion. This 
is an important factor, since such test programs usually 
extend over a period of several months, and, in some 
cases, parts of the specimen will be used for structural 
development a year or more after completion of the 
original test. This phase of the problem is further 
aggravated by the use of magnesium alloys in the struc- 
ture. 

(3) The time, expense, and complication of the sub- 
merging tank are avoided. 

Although we have found it more practical to use the 
air-pressure type of test, for the reasons already men- 
tioned, we still have a great deal of healthy respect for 
the potential energy stored in a pressurized cabin and 
therefore take precautions to prevent injury to per- 
sonnel or plant facilities. However, it has been our 
experience that it is possible to obtain a great deal of 





Fic. G. Details of the damage to the fuselage forward section 


caused by the pressure test. 
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design information in a reasonable length of time by 
repeatedly pressurizing the structure, with air, to sye. 
cessively higher pressures, reworking the structure each 
time a weak point develops. A considerable portion of 
this part of the program is usually accounted for in de. 
veloping the latches and seals for the doors. Usually, 
we try also to work in enough cycles of limit pressure 
to ensure that there are no details of the structure which 
are obviously deficient from the fatigue standpoint, 
When sufficient static and fatigue strength information 
has been obtained by this cut-and-try procedure, we are 
ready to determine the ultimate strength of the strue- 
ture. As we all realize, the failure that is then obtained 
can be rather violent, although this is not invariably 
true. 


per 
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Fic. H. Details of the damage to the fuselage forward section 


caused by the pressure test 


Fig. F shows a pressure test on the forward section of 
a DC-6 fuselage. 
of this specimen. 


Figs. G and H show the final failure 

At 10 Ibs. per sq.in., the vertical 
bulkhead immediately forward of the pilot’s position 
blew out. The pieces were thrown into the nose-wheel 
well and were stopped after barely puncturing the nose 
cap. There was recently a more dramatic demonstra- 
tion of the ability of a fuselage to hold together after 
explosive decompression. A DC-6 while in pressurized 
flight was cut almost in half by a propeller blade thrown 
from one of the engines. Except for one fatal heart 
attack brought about by the excitement, there were no 
injuries to passengers or crew, and the airplane was able 
to fly to the nearest airfield and land without further 
damage. 
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AIRCRAFT STRUCTURAL 


It is generally agreed that more accurate gust data, 
as well as a more suitable method of interpreting the at- 
mospheric gust phenomenon and airplane response, are 
urgently needed. In the U.S. we are presently engaged 
in this study with V-g-h-instruments of N.A.C.A. 
origin. These instruments record time, altitude, speed, 
and c.g. acceleration of the airplane. Most of the ex- 
isting data have been difficult to use because all of these 
variables have not been measured simultaneously. 

The N.A.C.A. and the aircraft companies of this 
country have asked the commercial air lines, as well as 
the military, to cooperate in obtaining flight data with 
presently available and, later, improved recording in- 
struments. In order to minimize the labor and delay in 
the reduction of the vast amount of data that would 
necessarily be collected, they have recommended that, 
as far as possible, the data reduction be mechanized. 

While some repeated load testing of aircraft wing 
structures has been, and is being, done in this country, 
we feel that data from fatigue tests of small specimens 
and components can be used to good advantage. For 
the past IS years we have been fatigue testing various 
types of specimens from simple riveted joints to full 
size wings and have been fairly successful in correlating 
the results of small specimen tests with those of tests of 
large complicated specimens. One of the test programs 
involving the 15-ft. fuselage specimens showed especi- 
ally good correlation between the two types of tests. It 
had been suspected that the attachment of the fuselage 
frames to some of the stringers would be critical in 
fatigue because of repeated pressure cycles in the fuse- 
lage. Small specimens simulating the attachment con- 
ditions were tested in a fatigue machine with results 
that were believed to confirm this feeling. The full- 
scale test that followed duplicated the results of the 
smaller test. It was then a relatively simple matter to 
devise several improved attachment designs, to try 
them out with simple specimen fatigue tests, and, 
finally, to incorporate the best redesign on the produc- 
tion airplane. (See Fig. K.) 

Dr. Walker has referred to the practical necessity of 
balancing economic considerations against design re- 
finement for maximum fatigue strength in structural 
joints. We also feel that such costly steps as reaming 
holes to extremely close tolerances, shot peening, using 
relieving slots, etc., should be used only after all of the 
simpler design refinements have been exhausted. For- 
tunately, there are steps that can be taken to increase 
the fatigue strength of parts without requiring too 
skilled workmanship. For instance, it may be worth 
while to use bolts having threads rolled after heat- 
treatment, which gives much greater fatigue strength 
than ground or lathe-cut threaded bolts. This is illus- 
trated in Fig. N. The effect of pretensioning of bolts 
that Dr. Walker referred to can also be seen on this 
figure. The test loads to which these bolt specimens 
were subjected were carefully controlled, and the load 


ranges were constant during the testing. It should not 
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be implied, however, that this increase in strength is 
indicative of what can be attained in an aircraft joint 


in service. Lab tests have shown that, after a few 
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thousand cycles of loading, the pretensioning that had 
originally existed in the bolt decreases and may de- 
crease to as little as one-half its value. Overhaul and 
maintenance checking every few thousand hours are the 
only preventative measures for the fatigue of major 
structural bolts in aircraft. 

In summarizing the fatigue problem, we agree that 
we need more knowledge of the flight loads imposed on 
aircraft. We do have a reasonable knowledge of the 
effect of cumulative loading. We also have some suc- 
cess in correlating small laboratory coupons with full- 
scale tests. We can also, with a fair degree of accuracy, 
calculate or predict fatigue life and strength of test 
specimens from past experience and a wealth of data; 
but, as to the calculation of the life expectancy of an 
aircraft, little progress has been made. Fig. O, the form 
of which you have seen many times, shows in some 
degree the dispersion in laboratory test mortality data 
and service loading mortality. Mortality dispersion of 
the joint between the inner and outer wing panels of the 
DC-8 are plotted, using laboratory test data and actual 
service failure data. Life insurance company data on 
human mortality are also plotted for comparison. 
Mortality of the laboratory test specimens and of hu- 
man life appears to follow a similar law, while that of 
airplane wing joints follows a different one. It is inter- 
esting that a calculation of the life of the part from the 
stress analysis and accepted gust frequency statistics 
revealed the part would never fail. Yet failure did 
occur in the form of skin cracking which is the most 
common type of fatigue cracking of light members. 
Serious? No, not when it takes 24,000 service-hours 
per airplane to develop cracks in 6 per cent of 920 wing 
joints. 

Fatigue cracks characteristically develop slowly, and 
finding them is a matter of careful inspection. Here, 
then, is the essence of concern for the eventual struc- 
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Determination of Elastic Wing Aerodynamic 
Characteristics 


WILLIAM J. GAUGH* anp JOSEPH K. SLAP?+ 
Northrop Aurcraft, Ine. 


SUMMARY 


The basic aeroelastic integral equation is formulated by treat- 
ing the wing as a tapered cantilever beam mounted normal to its 
straight elastic axis. By means of an iterative process, a series 
solution is obtained for the aeroelastic change in angle-of-attack 
distribution. For sweptback wings, this Neumann solution fails 
to converge at dynamic pressures that exceed the absolute value 
of the divergence dynamic pressure, even though the latter is 
negative and the wing physically stable. This eigenvalue prob- 
lem is handled by expanding the series about a new point and 
thus effectively increasing the radius of convergence to agree 
with the physical case. It is shown that the new point, expressed 
by a convergence factor, can be chosen so as to improve the rate 
of convergence of the series. 

The aeroelastic effects on the stability derivatives, on the hinge 
moments, and on the effectiveness of the control surfaces are in- 
vestigated and shown to be expressible as simple functions of the 
dynamic pressure. Wind-tunnel data for a full-scale sweptback 
wing are presented and agree well with results calculated using 


the present method. 


NOTATION 


a, dg = constants—Eq. (8.5) 

b = wing span 

bh, bs = constants—Edq. (8.6) 

c = chord, parallel to centerline 

Cy = root chord 

d = per cent semispan center of pressure location 

EI = bending stiffness for sections perpendicular to elastic 
axis 

e = distance of elastic axis aft of local aerodynamic 
center, measured as fraction of local chord 

GJ = torsional stiffness for sections perpendicular to elastic 
axis 

h = bending deflection of elastic axis, positive for bending 
up 

' = identity operator—Egs. (6.1) 

om = indices for series summation 

k = constant, used with various subscripts—-Egqs. (8.7), 
(8.8) 

l = length of elastic axis of one wing panel 

Mo = section lift curve slope for zero sweep 

P = linear operator, defined to operate on an angle of 
attack—Eq. (3.3b) 

P, = operator—Edq. (3.3d) 

p = rolling velocity 

p = rolling acceleration 

q = free-stream dynamic pressure 

ep = eigenvalue of least absolute value of Eq. (3.3a)—Eq. 
(5.1) 

S = total wing area 


* 


Presented at the Aeroelasticity Session, Ninteenth Annual 
Meeting, I.A.S., New York, January 29-February 1, 1951. 

* Aerodynamics Engineer. 

+ Aerodynamicist. 


s = dimensional coordinate measured from root positive 
along elastic axis of starboard wing panel. 


coordinates—Section (2) 


3% 
a = angle of attack in free-stream direction 
a), @ = successive approximations to final angle of attack 


distribution 
a = denotes loading independent of root angle of attack 
Eq. (3.3d) 
= angle of attack for zero lift 


ay 

a = see Eq. (3.3c) 

aR = angle of attack distribution before elastic deforma- 
tion 

i = dihedral angle 

A = used before other symbol to denote change 

6 = angular deflection of control surface 

0 = convergence factor—Eq. (6.2) 

A = sweep angle of elastic axis, positive for sweepback 

S = nondimensional spanwise coordinate—Section (2) 

¢ = angle of twist for sections perpendicular to elastic 
axis 

identically equal to 
- = approximately equal to 
(? = value of () for wing before deformation 
¢? = mean aerodynamic chord (M.A.C.) 


INTRODUCTION 


i pammegee A LIGHTWEIGHT WING with a low drag 
polar and a high critical Mach number requires 
developing a highly swept, thin wing having a large 
aspect ratio and designed for a low ultimate load factor. 
The large deformations encountered with such a wing 
produce large changes in the magnitude and distribu- 
tion of the air loads acting on the wing. These, in turn, 
cause large variations in the stability and flight attitude 
of the aircraft. Thus, the elastic wing structure 
created aerodynamic and structural problems without 
precedent and therefore without available methods of 
solution. To eliminate the necessity of building full- 
size wings for static and flight testing, several methods 
of solving the aeroelastic problem were investigated. 


One approach, the direct structural, which consists of 
applying the rigid wing air load, using the resulting 
angle-of-attack distribution to calculate a new air load 
and rm ating this procedure until there is no further 
chans .n the air load, requires an increasing number of 
adj iuuents as the dynamic pressure is increased. A 
method that not only speeds up and ensures the con- 
vergence of this process but also conveniently expresses 
the sum of the series expression required for mathe- 
matical treatment is developed herein for use with a 


given wing design. 
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Fic. 1. 


Three views of a swept wing with nomenclature. 


The elastic deformations of the wing caused by verti- 
cal acceleration (wing inertia or load factors) and various 
air-load conditions (elevon deflections) produce large 
changes in the spanwise distribution of the local lift co- 
efficients and angles of attack of the wing. In addition 
to these longitudinal displacement derivatives being 
altered by aeroelastic effects, the damping moment 
derivatives due to forward velocity and to rate of pitch 
must be corrected for wing deflection before the longi- 
tudinal equations of motion can be determined for any 
instantaneous flight condition. The rolling moment 
derivatives due to slideslip, due to rate of roll, and due to 
differential elevon position must also be corrected for 
aeroelastic effects before the lateral equations of motion 


can be obtained. 
‘ 


Two types of deformation occur in the plane normal 
to the elastic axis of a swept wing. One type is caused 
primarily by the deflections of the wing near the root. 
In this region, there is coupling between the bending 
and torsion so that bending induces torsion and vice 
versa. Also, for nonsymmetrical loading, there is 
carry-through from one panel (through the root) to the 
other because of the difference in the root bending 
moments of the two panels. Since these types of root 
deformation are not affected by the geometry of the 
wing away from the root, they may be assumed to be 
constant along the outer portion of the wing span. 
Methods of analyzing these effects are still being in- 
vestigated. However, it can be stated that, as the 
aspect ratio of the wing is increased, the root effects de- 
crease in relative importance. Static and wind-tunnel 
tests of a full-scale, high-aspect-ratio wing have verified 
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this conclusion. Moreover, for sweptback wings, neg. 
lecting these effects is conservative. 

The second, and aerodynamically more important, 
type of deformation is analyzed in the present report, 
that is, 
the assumed root section is mounted skewed so that it js 


The wing is treated as an unswept cantilever 
normal to the elastic axis. It intersects the actual root 
section at an angle equal to the sweep angle of the wing, 
This small change near the root, however, is not im- 
portant. The elastic axis of the wing is assumed 
straight. For the swept-wing case, an elastic axis curved 
toward the wing trailing edge at the root in the plane of 
zero dihedral would have been a more general assump- 
tion, but the mathematical complexities that are intro- 
duced are not warranted by the small increase in ae- 
curacy which might be achieved, especially in view of the 
other approximations involved. Experimental data 
show excellent correlation with theoretical results from 
the present method. 


(1) ASSUMPTIONS 


(A) Deformation of sections near the root due to 
geometric coupling between bending and _ torsion is 
neglected. 

(B) The cantilever wing panel root section is mounted 
(e.a.). Stiffness 
distributions for both bending and torsion can be de- 


normal to a straight “elastic axis’’ 


termined for sections normal to this e.a. 

(C) The spanwise loading acts along a line of aerody- 
namic centers (1.a.c.) 

(D) There is no deformation in the chordwise direc- 
tion. Also, the control surfaces are connected to the 
wing with a rigid hinge fitting across their span. 

(E) The aeroelastic change in the geometrical angle 
of attack and dihedral distributions at a given span 
station, &, for chord sections parallel to the plane of 
symmetry, are given, respectively, by 


Aa(é) = g(&) cos A — (Oh/0s)(€) sin A (1.1) 
AT(é) = ¢g(&) sin A + (O//0s)(~) cos A (1.2) 


(F) The effects of drag are neglected. 


(2) COORDINATE SYSTEM 


Standard N.A.C.A. notation and conventions are 
used where possible (Fig. 1). The stability axis system 
has its origin at the center of gravity of the aircraft, with 
x measured positive along the direction of flight, y 
measured normal to the plane of symmetry and positive 
to the right, and z measured positive down. The elastic 
axis 1s a body axis with its origin at the plane of sym- 
metry and with s measured positive out along the 
elastic axis of the right wing panel. The dimensionless 


span coordinate, £ is defined by 
Ss y s cos A 
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(3) ELastic EQUATIONS 


The total bending moment due to aerodynamic and 
other external loads applied at a given wing station may 
be equated to the internal elastic moment at that sta- 
tion. The total torsional moment applied about the 
elastic axis at this wing station may be expressed 


similarly. 


2] 
EI = -f ds» fia q(c, + €,,)e cos A ds; 
ds? Ss 


Og é . ° 
GJ a = q(Cmy + C0Ci, + €C,)C" COS* A ds, 
Ss s 


The moment and lift coefficient distributions, C,,, 
and ¢,,, are independent of root angle of attack. For 
example, Cm, might be the section moment due to cam- 
ber, while c,, might be the section lift due to elevon de- 
flection or to inertia loads. On the other hand, c, does 
vary with root angle of attack, being due to some aero- 
dynamic condition such as straight flight, rate of pitch, 
or rate of roll. The above moment and lifts are the 
equilibrium values for the elastically deformed wing. 

Introducing the dimensionless span station* and 
integrating the above equations give 


h © dé, . 
0 = gl* cos A [FZ & def c(G,, + &) dé 
Os 0 te 


(3.1) 


& dk 1 
g = gi* cos? A i oS f C?(Cmy + CC + €€r) dé 
0 GJ £2 


(3.2) 


An expression for the geometrical angle of attack can 
be obtained by substituting Eqs. (3.1) and (3.2) into 
Eq. (1.1). It is convenient to collect the loads which are 
not dependent upon the final equilibrium angle of at- 


tack distribution. 


a= arg + Aa = a + GPa (3.3a) 


where 


“§ (]? soe A 1 . sin A cos A 
Pa = | a f cec, dks EI x 
, dé. fo CC, det dés (3.3b) 


The term Pa is to be construed as the application of an 
operator P to an angle of attack a. First P acts on a 
to give c, (by some span loading method) and then on 
this c, to give Pa [by the operations indicated in Eq. 
(3.3b)]. “hat is, P consists of two distinct operations— 
namely, finding the c, due to a and integrating this ¢, 
after substitution in Eq. (3.3b). 


* The functional dependence of the variables on per cent span 
station will be implied, although the notation will be omitted for 


brevity—e.g., g(t) = ¢. 
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The term a» includes the initial angle of attack dis- 
tribution, az, plus the effect of any other condition, 
such as control surface deflection. Explicitly, 


a = are + gPia, (3.3¢) 


where Pia, is given by 


— |! 7? cos® A : 
Pia, = G (Cm, + CoC) dé — 
0 r 
[3 sin A cos A 
dt; (3.30) 


The term Pia, is similar in appearance to Pa but is 
somewhat different in concept. The multiple operator 
P, is defined to yield first c,, and ¢,, which do not de- 
pend on root angle of attack, and then to yield Pia, by 
means of the operations indicated in Eq. (3.3d). The 
term a,, then, is not an angle of attack per se but signi- 
fies a condition whose loading is independent of angle 
the loading of a deflected control surface 


. 


of attack—e.g., 
as given by c,,, and C),. 

The problem is to find that value of a which will 
satisfy Eq. (3.3a) once ap has been given—that is, to 
find the value of a for which the applied loads will be in 
equilibrium with the elastic forces of the wing. 


) ITERATIVE SOLUTION 


The solution to Eq. (3.3a) may be generated by 
means of an iteration process. Before using this ap- 
proach, it is important to note that P is a linear opera- 
tor. Then the first step in the solution is to assume 
that the equilibrium angle of attack at a given station 
is equal to ap. Substitution in the right-hand side of 


Eq. (3.3a) yields 
Qa, = a +- gP cw 


as the first approximation to the final answer. Next, 
assume that a is equal to a; and obtain 

a = a + gPa = a + gPa + q?P?a 
as the second approximation. Continuation of this 
process yields the solution in the form of a Neumann 
series—namely, 


a= = g'P! a 


i=0 


(4.1) 
where Pa) denotes i successive applications of the 
operator P, and P®ay equals ay. 

The discussion to follow in this and the next few 
sections is based on the premise that the initial loads 
If they are negative (e.g., for up elevon at 
then the conclusions must be 


are positive. 
zero angle of attack), 
modified. However, this will be obvious and will pre- 
sent no difficulties. 

For a positive initial load, the sign of both P‘a, and 
P'P,a, and, hence, of P‘ay will always be positive for 
sweptforward wings because both the bending and 
twisting of the wing tend to wash in the tip a. For un- 


swept wings, bending has no effect, while twisting 
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washes in a, resulting in all positive P’a . For swept- 
back wings in general,* the washout due to bending will 
predominate over the washin due to twisting, and 
therefore the sign of Pag and Pia, will be negative. 
Hence, the loading due to each of Pag and Pia, will 
tend to wash in the tip a, so that P*az and PP, will 
be negative. As more adjustments are made, the 
P'a, and P'P,a, will alternate in sign. 


(5) WING DIVERGENCE 


For each flight condition, the solution to Eq. (3.3a) 
is an angle-of-attack distribution for which equilibrium 
exists between the applied loads and the internal forces 
of the wing. 
dynamic pressure, whereas the wing strength, of course, 
does not. It is thus apparent that there might be a g, 
say qi, for which no equilibrium position exists—that is, 
A picture of the conditions 


However, the applied loads increase with 


at g, there is divergence. 
that exist at divergence can be obtained by considering 
an extremely flexible wing, such as one made of rubber. 
As g approaches g;, the deformations become extremely 
large, and at qg, the wing actually ‘“‘corkscrews.’’ How- 
ever, for any practical wing, as q; is approached, the 
large deformations result in stresses that exceed the 
values allowed by the strength of the material, and 
structural failure occurs. The concept of divergence 
here is analogous to that of the critical loads for an 
Euler column. 

For consistency, g; must be defined mathematically as 
that value of g for which a finite nonzero ap leads to an 
indeterminate value of a. Such values of g are the 
eigenvalues of the integral equation. That is, if ap is not 
equal to zero, there is, in general, no unique finite value 
of a which will satisfy Eq. (3.3a) when gq is an eigen- 
value.' Hence, if g; is the smallest (in absolute value) 
of the eigenvalues of Eq. (3.3a), then Eq. (4.1) will not 
converge for g = g:. Again using the Euler column as 

* For slightly sweptback wings, it is possible to have the washout 
due to bending be offset approximately by the washin due to tor- 
sion so that only the dihedral has an appreciable aeroelastic effect. 
The condition “in general’’ will be implied, but omitted, in the 
following discussion. 
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an analogy, the critical loads are found as the eigen. 
values of the equation that describes that problem. 

In reference 2 it is shown that the eigenvalue of lowest 
absolute value can be found from 
qi = lim (P"a/P" 


NAA 


+1 » 
a) (2.1) 


As was pointed out in the previous section, the P” ay are 
all positive for unswept and sweptforward wings, so that 
for such wings there is a real physical divergence speed,} 
corresponding to g, > 0. On the other hand, the P"qy 
alternate in sign for sweptback wings, so that q, < 0) 
and no real divergence speed exists. 

In order to calculate g,; from Eq. (5.1), the ratio 
P"a/P"*'a) must be calculated at several span sta- 
tions. The calculations should be done for successive 
values of m until the ratio becomes approximately con- 
stant across the span (Fig. 2). It usually will not be 
necessary to consider m greater than three or four. 

The question naturally arises as to which value of a 
to use in the calculations for q;. As has been stated 
previously, g; is an eigenvalue of the operational equa- 
tion, Eq. (3.3a). The eigenvalues, however, belong to 
the operator itself, in this case P, rather than to the 
function upon which it is operating. Hence, from a 
mathematical point of view, once P is defined, the eigen- 
values are determined and are independent of a. The 
same conclusions can be reached by physical reasoning. 
For the purposes of computation, then, any nonzero 
distribution can be used for ap to ascertain the succes- 
sive ratios P"a/P"*'ao for Eq. (5.1); but, with an eye 
toward future calculations, it is convenient to use the 
constant distribution ay) = 1 rad. Further discussion 
of the actual calculation procedure is to be found in 
Section (7). 

It already has been shown that, for unswept and 
sweptforward wings, mathematical divergence of Eq 
(4.1) and physical divergence of the wing both occur at 
g = 4: > 0. For sweptback wings, no physical diver- 
gence is possible because g; < 0, but mathematical 
divergence of Eq. (4.2) does take place when g > /q.. 
It is apparent, then, that Eq. (4.1) must be modified in 
some way in order to obtain a true picture of the 
physical situation for sweptback wings. Such a modifi- 
cation is developed in the next section. 


(6) EXTENDED SOLUTION 


Consider Eq. (3.3a) 


a= a+ Pa (3.3a) 
UI — gP)a = a (6.1a) 
a = (I — gP)—'a (6.1b) 


where operation on a function by the identity operator / 
leaves the function unchanged. The inverse in the 
right-hand side of Eq. (6.1b) can be rewritten 

{ The indicated air speed corresponding to q; is called the 
divergence speed of the wing. 
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(—gP)'=(1 a{l+(1/a)[( — @)J — gP)} -! (6.2) 


where @ is some fixed nondimensional number.  Ex- 
panding the inverse in a series and inserting in Eq. 


(6.1b) yield 


fl, 1S fO- ‘| a |! 
a= frt+5d ( A aa ar tal ee 


(6.3) 


where the term [J + (¢/@ — 1)P]' is to be expanded by 
the binomial theorem. 

Since Eq. (6.2) expresses an identity, it can be stated 
that, if Eq. (6.3) converges, it does so to the correct 
solution of Eq. (3.3a). It is also possible to show that, 
whenever the wing is physically stable, a @ can be found 
such that Eq. (6.3) will converge, whereas if the wing 
diverges physically, no such @ exists. It is interesting to 
note that, if @ = 1, Eq. (6.3) reduces to Eq. (4.1), as can 
be seen by an inspection of Eq. (6.2). It is thus empha- 
sized that the convergence of Eq. (6.3) is strongly 
influenced by the choice of 6. 

The problem of finding the exact effect of @ on the 
convergence of Eq. (6.3) is greatly complicated by the 
complexity of the operator in Eq. (3.3b), especially if 
an air-load method is used which takes induction effects 
into account. Convergence criteria can be set up using 
function theory, and these will provide appropriate 
However, the problem of finding the 
optimum rate of con- 
A numerical 


choices for 6. 
value of @ which vields the 
vergence permits of no facile solttion. 
approach was deemed feasible, and it indeed provided 
values for 6 which yielded convergence where physically 
possible; but it, too, failed to give generalizations con- 
cerning comparative rates of convergence. 

It is convenient to rewrite the solution in terms of 
Aa, the aeroelastic change in the angle-of-attack dis- 
tribution. Then, the actual angle-of-attack distribu- 
tion can be found by adding Aa to the initial angle-of- 
attack The equation for Aa may be 


written 


> Ma P \' 
aa = 247+ > ( ) (r- W\l x 
al a 4 1—o6/ 


(Pra. + Par ) (6.4) 


distribution. 


By means of a bilinear mapping in which g was con- 
sidered as a complex variable, O. K. Smith, of Northrop 
Aircraft, originally and independently extended the 
radius of convergence of the Neumann series. Smith 
suggested a particular value of his transformation 
parameter, and it was found to yield solutions with 
rates of convergence satisfactory for engineering pur- 
poses. Because of the isomorphism between Smith’s 
general result and that of the present report, the par- 
ticular value of the parameter mentioned above can be 


expressed as a special value of @—namely, 


6 = (1 — Q)/h (6.5) 
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This value of @ has been used in calculations that show 
good agreement with static and wind-tunnel experi- 
mental data. Substitution of Eq. (6.5) in Eg. (6.4) 


yields 


x i+] 
¢ 
Aa = -—n>, ( d ) CU —_ aP)'(Pia + Pap) 
=O \g —- qi 
(6.6) 


Eq. (6.6) is to be used when g; < 0. If gq > 0, it is sug- 
gested that Eq. (7.1) be used. Then, in both cases, 
mathematical convergence will occur over the same g 
range as for physical convergence. 


(7) SUPERPOSITION OF EFFECTS 


Successive approximations to the final solution for 
Aa, the aeroelastic change in angle of attack, are ob- 
tained by taking increasingly large numbers of terms of 
the series 


x 


q'P'(qPia + qPar) (7.1) 
7 


i=0 


Aa = 
if g, > O, and 


Aa = —n” > ( = 


i+1 
) (J = “P)'(Pice 4- Par) 
i=0 \d—- hh 


(6.6) 


if g, <0. Because of the linearity of the operators P and 
I, the principle of superposition can be employed. That 
is, the unit load for any condition (e.g., load for ag = | 
rad., or 6 = 1 rad., or 1g condition, etc.) can be applied 
and the results plotted versus span station for various q 
values. Superposition of different multiples of these 
loadings can be used to calculate the aeroelastic effects 
for any flight condition. By separating the effects in 
this way, considerable simplification is obtained in the 
numerical work. 


The final lift distribution can be obtained by sum- 
ming the distributions corresponding to the incremental 
angle-of-attack adjustments. The eigenvalue can be 
determined from the limiting value of the ratio of the 
lift distributions corresponding to two successive angle- 
of-attack adjustments (Fig. 2). 

(8) VARIATION OF WING LIFT AND MOMENTS WITH 
DYNAMIC PRESSURE 


Eq. (5.1), for finding the divergence dynamic pres- 
sure, g;, suggests an approximation that simplifies the 
expression for the Neumann series, Eq. (4.1). Let 


P”“ a p™+ | ms a q1 (8.1) 


where is sufficiently large so that the ratio is approxi- 
mately constant across the span. M, say, will be 
three or four. Then Eq. (8.1) implies 


P!@t"q, = Pa qi” -.) (8.2) 


ae 


(n = 


Hence, Eq. (4.1) becomes 
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k,>I| SWEPTBACK | a = c + [big/(b: — q)] (8.6 
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Fic. 3. Aeroelastic correction factor versus dynamic pressure. 
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a~a+ E: q'P' a + qu P™ a pe (q ‘q)? (8.3) 
Fray j= 
which, for |q| < |q:|, reduces to 
M-1 M pM 
ore ng PP % 
a~ at z q' P* a + i (8.4) 
i=1 a ¢ 


The above simplification can be extended by assum- 
ing that the final spanwise shape of a is approximated 
by 
a + [qgPao/(m — 9)] 


ant 
which, because of the simplifying assumptions made, 
can be written in the more general form 





a = a + [aig/(a2 — g)] (8.5) 

The resulting section lift distribution can be expressed 
as 
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Fic. 4. Comparison of predicted and experimental spanwise 
bending, angle of attack and twist. 


where ¢;/’ is the lift due to a. This lift distribution* cap 
be integrated across the span to yield the total lift 
moments, and static stability variations with g. 


Cr/Ci® = 1 + [Risg/(Ri2 — @)) (8.7a) 
Cu/Cu® = 1 + [Ranig/(Rue — gq) | (8.7b 
C/C& = 1+ [kuyg/(Re — @)!} (S.7¢ 


Cu L 


Cu ng 


_ (bs - 


- (&8) (@) 
‘i Cu Cr, - (Rus — 


Full-scale wind-tunnel tests of a swept wing have 
yielded results in good agreement with those predicted 
using Eqs. (8.7) and (7.1). 


q) [Rare + q(ku, = 1) ] 
q)(Ri2 + g(Rku —1)] 
(8.7d) 


Consider the aeroelastic correction factor 


k=] + [Rig (Re a q) | (S.S) 


which appears in Eqs. (8.7). On Fig. 3, R is plotted 
against g, with the assumption that kz = q:. In prac- 
tice, k. generally differs from q:, but the assumption of 
equality is sufficiently accurate for such trends as are 
shown on Fig. 3. The indicated air speed, correspond- 
ing to the g for which R = 0, is called the reversal speed. 


This g can be estimated from 
Grn; = ke (1 — ky) 


In order to obtain the expression for the correction 
factor R, the constants k, and ki, which appear in Eq. 
(8.8), must be determined. After the equilibrium angle- 
of-attack distribution, a, has been found for several 
extreme values of the flight g, the Cz, Cy, etc., corre- 
sponding to each of these a distributions can be divided 
by C,°, Cy", etc., the values for the undeflected wing. 
R is then known for each parameter at the two given 
values of g. These two values permit k; and k» to be 
found for each of Eqs. (8.7) in turn. The values of k; 
thus obtained will agree only approximately with the q 
found from Eq. (5.1), since the derivation of Eqs. (8.7) 
involves several simplifying assumptions. For this 
reason, Eqs. (8.7) and (8.8) have been written in the 
more general form, using k2 instead of q. 


It should be pointed out that Eqs. (8.7) and the con- 
cept of the correction factor R can be used validly even 
if a method other than that of the present report is 
employed for calculating the equilibrium angle of attack 
distribution. 

* The various c; distributions needed can be obtained from 
strip theory, Weissinger lifting-line induction theory,** or any 
other process for calculating span loading. Where a reasonably 
refined estimate of the static stability variation is required, lifting- 
line induction theory should be used. However, strip theory 
gives force variations that are usually consistent with preliminary 
estimates of the wing strength properties, such as the estimate 
that the wing torsional and bending stiffnesses vary as the fourth 
power of the chord. 
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ELASTIC WING 


(9) THEORY AND TEST CORRELATION 


A full-scale wind-tunnel test was made on an elastic 
wing having a reasonably large panel aspect ratio and 
angle of sweepback. The wing was mounted shoulder 
high on a long cylindrical fuselage. One-foot high, 
streamlined posts were located perpendicular to the 
upper wing surface along the leading- and trailing-edge 
chord lines. A camera was located in the aft portion of 
the fuselage to photograph the relative changes in the 
post deflections as the dynamic pressure, root angle of 
attack, or control surface angle was changed. The ex- 
perimentally determined changes in the bending deflec- 
tion, angle of twist, and angle of attack relative to the 
assumed elastic axis location (40 per cent chord line) 
were compared with the theoretical results. The com- 
parison is shown on Fig. 4. The theoretical results are 
based on strength distributions obtained from full-scale 
static tests for an induction air-load distribution. 

The variation with dynamic pressure of the total 
determined experi- 
As on 


force and moment coefficients 
mentally and theoretically is shown on Fig. 5. 
Fig. 4, the results indicate excellent agreement. 


(10) SwePT-WING STATIC TRIM REQUIREMENTS 


The variation with g of the trim requirements for an 
elastic swept wing can be investigated by using the con- 
venient R factor developed in Section (8). Rather than 
complicate the analysis by attempting to account for the 
aeroelastic effect on the rate of change of downwash on 
a horizontal tail elevator, the trimming moment can be 
obtained from a control surface located along the 
trailing edge of the outboard region of the swept wing. 
In order to simplify further the analysis, the lift and 
moment produced by the tail and body will be neglected. 

The change of the moment coefficient of the wing with 
the root angle of attack, Cy,, can be expressed relative 
to the quarter-chord of the M.A.C. as a function of the 
center of gravity position, a/c.. The form of the 
aeroelastic correction factors appears as Eqs. (8.7a) 
and (8.7b). The subscript a on the bracket denotes 
the dependence of the factor on the root angle of at- 


tack. 


: : *) ( a ) ( <t) 
Cue = Cue. (= —~ (—)c,.°(— (10.1) 
M M 6 (' u! - Ce L ( o o 


For reasonably large amounts of sweep (35° or 
greater), the lift generated by the deflection of the con- 
trol surfaces on the wing will produce a moment about 
the center of gravity which, by comparison, makes it 
possible to neglect the relatively small effect of the aero- 
dynamic couple due to the effective camber produced by 
their deflection. The elevator effectiveness can then 
also be expressed as a function of the center of gravity 
location relative to the quarter chord of the M.A.C. 


’ Cu ) (= ) ( Ce ) 
c.. =c,, (&") — (%\c0(&) aoe 
Mé Mb, 25 (= 5 Ce ” C,° 6 : | 
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slope and static stability. 


A convenient method of expressing the effect of the 
longitudinal location of the wing center of pressure rela- 
tive to the center of gravity of the aircraft can be ob- 
tained by dividing Eq. (10.1) by the lift curve slope of 
the wing. The form of the aeroelastic correction factor 
has been given as Eq. (8.7d). 


Cu, = Ki — (a/c) (10.3) 
where 


°(Cur,/ Curr) a 


% 


KA; = Coy, = 


Cc M Lo 


The lift and moment equations relative to the center 
of gravity can be written as follows: 


Cral(a — on) + Cro = Cr (10.4a) 


Cuala@ — ap) + Cu,6 = —Cu, (10.4b) 
where az» is the angle of attack for zero lift. 

These equations can be solved simultaneously to de- 
termine the amount of control surface deflection, 6, and 
root angle of attack, a, required for trimmed constant- 
speed flight for a given lift coefficient. 


6 = K2(Cy,Cr + Cu,) (10.5) 
where 
Ke = 1/(Cuz,,Crs — Cms,.,,) 
a= a, + K3C, — K4 
where 
K; = 1/C,,; Ky = K;C_, 


If the undeformed swept-wing values of C,,", C;,°, 
op Cu, m and Cy, are known for a given Mach 
Number and if the aeroelastic correction factors have 
been evaluated, then Cy,,, 6, and a can easily be calcu- 
lated for a given value of C; and dynamic pressure by 
reading the values of K;, Ke, K3, and A, from plots 
similar to Fig. 6, which was constructed for a typical 
sweptback wing design. 
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It is interesting to consider the case for which Cy,) and 
az. can be neglected. Then Eqs. (10.5) and (10.6) for 6 


and a become 
6/C, = KeCu, = K2[Ki — (a/c) ] (10.7) 
and 


a/Cy = K3(1 — KeCrsCu,) 
= K3 1] ie K2Cz5[Ky —_ (a Ce 


(10.8) 


Fig. 7 indicates the trend of the trim requirements 
with indicated air speed and dynamic pressure for a 
range of values of Cy, and c.g. locations. 

Of course, to trim an infinitely rigid sweptback wing, 
the elevators must be deflected further up and the root 
angle of attack must be increased as the value of the 
static stability Cy, is increased negatively or the center 
of gravity is moved forward. However, Fig. 7 indicates 
that for flight at a constant value of C, and Cy,, the re- 
quired up elevator deflection increases and the root 
angle of attack increases as the indicated air speed or 
dynamic pressure is increased. 


(11) SrasiLtiry DERIVATIVES 


If the lowest natural frequency of the elastic wing is 
greater than four or five times the natural frequency of 
the longitudinal short period oscillation of the aircraft, 
then a reasonable first approximation to the dynamic 
response of the aircraft can probably be obtained by 
neglecting the coupling effects that would be introduced 
by considering the additional degrees of freedom due to 
the dynamics of the elastic wing. In the present in- 
vestigation the aeroelastic corrections to the undeformed 
wing have been analyzed on this quasi-stationary basis 
That is, the spanwise mass distribution of the wing is 
assumed to experience the same dynamic response as 


the center of gravity of the aircraft. 

The aeroelastic effects on the longitudinal X axis 
derivatives are neglected. 

If compressibility effects are neglected and the air- 
craft is assumed to be in a steady-state trimmed flight 
condition such that the angle of attack, a, and the 
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elevator deflection, 5°, are constant, then a perturbatioy 
change, u, in the steady-state forward velocity, U, wil 
produce a vertical force, 2,upUS, and a pitching mo. 
ment, m,upUSc., where 2, and m, are nondimensional. 


1 dZ “A nW a2) dq 


™ pUSdu pU\ gS * dg/ du 


nW ( 9 dCi Pa 5 acu) 
gS i dq dq ie 


I dM q dCydq _ 


%, = —— . 
pUSc, du pU dq du 


1C. IC. 
q ( — 5 go ut) (11.2 
dq dq 


The variation of the displacement derivatives, (C,,, 
Crs, Cua, and Cys, with dynamic pressure can be ob- 
tained by differentiating the suggested algebraic form 
for the aeroelastic correction factor given in Eq. (8.8). 


dC . d kiq kyk2C 
= ( (: > = - (ina 
dq dq ko —q (ke — gq)? 


Substituting Eq. (11.3) into Eqs. (11.1) and (11.2) gives 
the following expressions, where 7 is the load factor. 


nW Ba | kisku 
Zu = — —_— ‘ YH 5 2 
qs — (Riz — 9)" Je e 
. Risks | 
6°C 0 2 
a — gq)’ 5! 
k k 9 
mM, =  q erCue!| aael 4 
M2 a 
" Rukue | e) f 
OC, 0 « 1 7” ; 0C ‘ 
ss Fe = q)” ' q (* a La x 


Rrikre | © | Riikie it 
&°C 0 > . 
r —~ 2" de . - (Ria — q)* Je) 


Vertical acceleration of the spanwise mass distribu- 
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(11.5) 


tion of the wing has an effect upon the lift and moment 
coefficient expressions. This effect can be calculated by 
determining the angle-of-attack distribution due to 


wing weight (corresponding tom = land g = 0). The 
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resulting lift and moment distributions can be integrated 
across the span in order to determine C,,° and Cy,,’. 
Their variation with g can be determined by using the 
iteration procedure previously discussed. The results 
can be expressed as follows: 


k 
Cr, = Cz,° (: + - — ) (11.6) 
kre — Qn 


k, 
Cu, = Cutng s (1 4 mq ) 
Rue = 


The change in the vertical force and pitching moment 
due to a change in the vertical velocity includes an 
aeroelastic correction term. 


_— Cra _ = Cie ) 
a 2 2 \nW/qS 


= -(#)'[e +<¢ ( Cra )| (11.9) 
a sh vedi . 


For a conventional aircraft with an unswept wing, the 
horizontal tail is the main contributor to the damping 
in pitch. Consequently, the effect of the damping in 
pitch for a swept wing alone has received little atten- 


a 
: Cry, (11.7) 


n Ce 


(11.8) 


tion. 


my = (1/4)(¢-/Re)* Cure (11.10) 


Before an aeroelastic correction factor was deter- 
mined for mj, an assumed lift distribution along the 
quarter-chord line of the undeformed wing was derived 
for an induced angle-of-attack distribution along the 
three-quarter chord line due to the angular pitching 
velocity of the wing relative to the center of gravity 
location of the aircraft. When this lift distribution was 
integrated across the undeformed wing span, it gave 
somewhat better agreement with some experimental 
data” than that predicted by strip theory. 


: ae C ae C, 
Cua = Afi May Py (S) = fol Mita, hy (=) sig 
— Se) — (=) . (=) | 
he E- foCre | = — ate x 
| as (< ore eat Cr) t Ha Cu’ a 
a ° C, a\? 
~ Star i 
Ce C1" a \Ce 


fi = 1/e. (co — b& tan Ao.25) 
f. 


) 


(1E557) 


where 


II 


b/c, tan Ao.75 


The small letter ¢ subscript denotes the effect of a 
linear twist distribution. Notice that Eq. (11.11) ex- 
presses Cy aS a quadratic function of the distance be- 
tween the quarter-chord of the M.A.C. and the center 
of gravity of the aircraft. 

The change in pitching moment due to vertical ac- 
celeration, m,, was also investigated. 


my = (1/4) (c/Re)*® Cura (11.12) 


The virtual air mass moving with the wing was as- 
sumed to be proportional to the area of a circle described 
by tlie local wing chord. The inertia reaction that was 
assumed to act along the mid-chord line was integrated 
across the span to obtain the following expression for 
the moment coefficient relative to the center of gravity. 


, | : a Cr : 
Cuz = Cun. ~ (*) Cr2 (=) (11.13) 


where 
‘ -0 Cr 7 0 Cur 
Crt arg 2 = faCra (&) r Cua (=). 
fs = (1/4 ¢) (2b tan Ao.25 — Co) 


The subscript 0.5 denotes mid-chord loading, and /; is a 
function of the wing geometry. 

It might well be argued that analyzing the damping of 
the pitching moment is somewhat facetious, since it is 
known that the terms are small and generally will not 
affect the short-period response. However, since this is 
not always true and since the purpose of the investiga- 
tion is to indicate an approximate method for correcting 
those stability derivatives that are affected by aero- 
elasticity, these terms were investigated. For a given 
wing, it might be desirable to determine the magnitude 
of the effect on the longitudinal response of the aircraft 
for a range of estimated values of Cyj and Cy, before 
attempting to evaluate the above expressions. 

The aeroelastic effects on the side forces, yawing 
moments, and the rolling moment due to rate of yaw 
lateral stability derivatives for the elastic swept wing 
are neglected. 

The elastic swept wing damping in roll due to rolling 
can be determined by integrating the final equilibrium 
value of the local lift coefficient across the semispan 
relative to the X axis. If the method discussed pre- 
viously is used, the undeformed wing local lift coefficient 
distribution, c,°, should correspond to the antisym- 
metrical case where each wing panel has an antisym- 
metrical linear twist distribution expressed as a func- 
tion of the nondimensional parameter, (p)/2U).’ It is 
also convenient to evaluate the aeroelastic correction 
factor [Eq. (8.7c)]. 


(C,/C?)» - {1 + [kng/Ri2 gh» 


Because of the mass distribution of the wing, rolling 
acceleration will produce an antisymmetrical loading 
across the wing span. The corresponding antisym- 
metrical local lift coefficient distribution can then be de- 
The aeroelastic correction factor can be ex- 


(11.14) 


termined. 
pressed as follows: 


(C/CP)5 = {1 + [kng/(Rn — 1} 5 
The method can also be used to determine the aero- 
elastic effect on the aileron effectiveness. It is con- 


venient to apply the assumption of linearity and 
superposition of results as discussed in Section (7) 


(11.15) 
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The undeformed swept-wing initial angle of attack and 


corresponding lift distribution are taken to be zero, so 


that the initial deformation is due only to the local 
lift distribution corresponding to the aileron deflection, 
c,°. The final equilibrium antisymmetrical lift dis- 
tribution can be integrated across the semispan in order 
to evaluate the rolling moment due to aileron deflection. 
The aeroelastic effect can then be expressed as a form 
similar to Eq. (11.15). 


(C,/C)s = {1 + [kng/(Re — Q)]}s (11.15) 


For a given Mach Number, an expression for the 
aeroelastic effect on the steady-state rolling effective- 
ness, pb/2l’, can be obtained by dividing the rolling 
moment coefficient produced by a given amount of 
differential aileron deflection, Aé, by the damping in roll 


( 4 kung ) 
pb AéC 18° Ris — OF 5 ~ 
i. Bee - (11.17) 
? rm (1 4 ng ) 


Rie ae q 


of the elastic wing. 


An equation for the integrated results for the effects 
of angle of sweep and dihedral on the rolling moment 
coefficient due to sideslip, Cig, is given in reference 8. 
This equation can be used to illustrate the form of the 
aeroelastic correction. The effect of the skewness of the 
trailing vortices is assumed negligible. For a given 
elastic wing, the dihedral can be expressed as a simple 
function of the dynamic pressure.* A method of es- 
timating the effect of dihedral on the air-load distribu- 
tion on unswept wings is presented in reference 9 and on 


sweptback wings in reference 7. 


i k, 
Cig =a - (1 oa k nf ) ' Cr, tan Ao.25 + 


2 d2 — | 


E + My COS a 
mw ALR. dC,° 


kid : 
: 1 =f. lm x 
2m cos Ao. | da Ris — g 
l + 


w A.R. 
kriq )! 
1 
( ai kro — q f 


The method of determining the k constants has been 
discussed previously. 
semispan center of pressure location d should be deter- 


The lift curve slope and per cent 


mined for asymmetric loading. However, the values 
obtained from symmetric loading should provide rea- 
sonable accuracy. 


(12) EFFECT ON THE HINGE MOMENTS 


The hinge moment coefficient is dependent upon the 
magnitude of the control surface deflection, 5, and the 


* The form is the same as that for angle of attack [Eq. (8.5).] 
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effective angle of attack, a,, in the region of the contro] 
surface. For an elastic wing at a given root angle of 
attack, ag, a change in g produces an incrementa| 
change, Aa,, which can be approximated as follows: 


Aa/ ar ~ [kig/(k2 — 9) }a (12.1 


Similarly, for zero root angle of attack, the a, pro- 
duced by a control surface deflection, 6, can be approxi- 
mated as follows: 

Aa, ~ [kig/(k2 — q)]s6 (12.2 
The constants k, and ke for each condition (denoted by 
subscript a or 6) are evaluated by the usual procedure 
[Section (8) ]. 

The conventional hinge moment coefficient equation 
-an be expressed in terms of the root angle of attack, 
the control surface deflection, and an approximate 
correction for the aeroelastic effects. 


al j kiq 
C, = Go + \ l + are + 
Re — q/a 


( kiq ) sb C. + aCe (128 
kp ‘a g ’ ( ha hé <0) 


Although these indicated corrections are only ap- 
proximate, they do indicate trends and relative orders 


of magnitude. They should be useful in flight-test 


analysis of hinge moment data. 
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A Theory of the Direct and Inverse Problems 
of Compressible Flow Past Cascade of 


Arbitrary Airfoils 
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ABSTRACT 


A simple unified approach to the inverse and direct problems 
of compressible flow past cascade of arbitrary airfoils is pre- 
sented. In this approach, the close relation between the shapes 
of the mean line of the blade and mean stream line in the channel 
and that between the variations in channel width and specific 
mass flow on the mean stream line are utilized. The customary 
assumption of a linear relation between pressure and specific 
volume is not required in this approach, and the camputations in- 
volved are also relatively simpler and quicker. The method is 
particularly useful for high-solidity blades. 

In the inverse problem, with the inlet and exit angles, a desir- 
able blade thickness distribution, and either a desirable mean line 
of blade or a certain mean stream-line shape chosen, the flow 
along the mean stream line can be properly determined and easily 
extended out along the pitch direction by the equations of con- 
tinuity and motion. The blade boundaries are then determined 
from the inlet mass flow and best velocity distribution over the 
blade. For a quick approximate solution of the direct problem, 
the same process is used with one or more adjustments on the 
mean stream line to get the correct blade shape. The solution 
obtained in this manner can also be used to give a good starting 
value for a more accurate solution by either relaxation method for 
hand computation or matrix method for automatic machine 
computation; these are also briefly described. 

The theory is illustrated with solutions of both inverse and di- 
rect problems of compressible flow past typical turbine cascades 
of highly cambered thick blades in which accurate results are ob- 
tained by using three terms in the series. The results compared 
well with experimental data. The method is directly applicable 
to two-dimensional flows on cylindrical surfaces in axial turbo- 
machines and can be extended to two-dimensional flows on arbi- 
trary surfaces of revolution in radial and mixed-flow turboma- 
chines. Families of blade sections for these machines can be 
built up by this method in a relatively simple manner. 


INTRODUCTION 


A BASIC AERODYNAMIC PROBLEM of turbojet and 
turbopropeller engines is the flow of compressible 
fluid past a series of blades in circular arrangement. In 
axial-flow type turbomachines, if the blades are rela- 
tively short in the radial dimension and are bounded 
by cylindrical walls, the theoretical flow passing 
through the blades is usually computed on the basis of 
two-dimensional flow on cylindrical surfaces, which are 
developed into planes for convenience of calculation. 


Presented at the Propulsion Session, Nineteenth Annual Meet- 
ing, I.A.S., New York, January 29-February 1, 1951. 

* Aeronautical Research Scientist. Now, Professor of Me- 
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neer, Minneapolis-Honeywell Regulator Company. 


A number of methods have been proposed to obtain 
the theoretical flow through such a given cascade of 
airfoils or to obtain a desirable airfoil shape for the blade 
section, including analytical methods using conformal 
mapping, interference technique, and series expan- 
sion;'~" graphical procedure ;''~'* and other mechani- 
cal and electrical devices.'* 


In the present paper a simple unified approach to the 
direct and inverse problems of compressible flow past 
cascade of arbitrary airfoils is presented. In both prob- 
lems, calculation is first made for the flow along a par- 
ticular stream line in the channel formed by two neigh- 
boring blades, preferably the mean stream line that 
divides the mass flow in the channel into two equal 
parts. In this calculation, the close relation between 
the shapes of the blade mean line and the mean stream 
line and that between the variations in channel width 
and specific mass flow on the mean stream line are 
employed. The flow is then extended in the pitch di- 
rection by the use of Taylor series, whose successive 
terms are obtained by the use of the equations of con- 
tinuity and motion. In the inverse problem, the blade 
boundaries are determined by the given mass flow at 
the inlet, and, by interpreting the starting mean stream 
line as dividing the mass flow in the channel into two 
slightly different portions, a number of blade profiles 
with different velocity distributions are obtained and 
the best one is chosen. In the direct problem, suc- 
cessive correction of the shape of, and the flow on, the 
mean stream line is necessary to match the given blade 
shape. Because the computation involved is rela- 
tively short (in the turbine example given herein, less 
than 16 hours were needed for one complete calcula- 
tion), iteration in the inverse problem for most desir- 
able blade thickness distribution or velocity distribu- 
tion and iteration in the direct problem to fit the given 
blade shape are practical. This iteration, however, 
will be reduced to a minimum if a number of families 
of cascades are built up by this method. A check of 
the accuracy of the solution in either problem and its 
improvement, if found desirable, can be readily made 
by finite difference methods using the available flow 
variations obtained in the solution. 


Several numerical examples are given at the end of 
the paper to illustrate these methods. Results ob- 
tained are compared with available experimental data. 
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SYMBOLS 
B = differentiation coefficients 
[C] = coefficient matrix 
[D] = diagonal matrix 
E] = [L] + [uw] -— [7] 
H = total enthalpy, h + '/2W? 
h = static enthalpy 
[7] = unit matrix 
J = VW ,(0p/ds) — W-(dp/dy) 
K = constant 
L = blade length in the ¢ direction 
[L| = lower triangular matrix 
{/] = lower triangular matrix with unit diagonal elements 
M == mass flow 
N = number of interior grid points 
P = pitch or spacing 
p = static pressure of gas 
S = stream line 
t = blade thickness 
{u] = upper triangular matrix with unit diagonal elements 
W = velocity of gas relative to blade 
y = distance in the direction of pitch 
Zz = distance in the direction perpendicular to pitch 
ia}, {8B}, fv}, {6}, fe} = column matrices 
B = flow angle, tan~' (Wy/Ws) 
y = ratio of specific heats 
A = correction in pw, 
6 = grid spacing 
p = Static density of gas 
yy = stream function 


Superscript 


* = dimensionless value 
Subscripts 

7 = at inlet 

e = at exit 

m = mean stream line 

p = pressure surface of blade 

s = suction surface of blade 
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FLow ALONG THE MEAN STREAM LINE 


The steady two-dimensional isentropic flow of a non. 
viscous compressible fluid is governed by the following 
equations of continuity, irrotationality, and the isep. 
tropic pressure-density relation: 


[O(pW,)/dz|] + [O(pW,)/dy] = 0 (| 
(OW./Ov) — (OW,/d0z) = 0 (9 
p = Kp’ (3 


Consider now the gas flow along a stream line some- 
where in the mid-part of the channel, such as ab jn 
Fig. |. The coordinates of the stream line and their 
differentials are related, respectively, by the following 
two equations: 


S(z, y) = 0 (4 
(OS/0z)dz + (O0S/dy)dy = 0 (5 


Following the motion on this stream line, it is con- 
venient to consider any quantity g on the stream line 
as a function of z only—that is 

q = gz, y(z)] (6 


The total derivative of g with respect to zis 


dq gq , Og dy 
dz Os oy dz 
But 
dy 0S/dz W, 
-—=— = — = tan 8 (9 
dz Os/Oy W. 
Hence, 
dq Og . Og 
dz Oz ° Oy si 


Using Eq. (9), the continuity and irrotationality con- 
ditions can be written as 


l(pW. O(pW,) O(pW, 
= ev) _ Kl) tanB=0_ (10 
dz Oy Oy 
and 
dW, OW. OW, 
— _ - tan B = 0 (11) 


dz Oy oy 


The variation of density throughout the flow domain 
can be most conveniently expressed in terms of 1ts inlet 
value, through the use of Eq. (3), as follows: 


p _ pr r — 
x _ _ va p os er ( ) 12) 
p= = = 1 — (12 
Pi Pi PT,i Pi 2H 
where // is a constant throughout. A tabulated func- 
tion of p/pr,; in equal intervals of W?/H can first be 
calculated, from which either a table for p/p; in equal 
intervals of W? or a graph between the two can be easily 
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FLOW PAST CASCADE 
constructed for each individual case and used for the 


evaluation of density from the velocity. 


In a solution of the isentropic flow through a typical 
cascade of turbine blades with two inlet Mach Numbers 
bv a numerical method to be described later, the fol- 
lowing results are obtained: (1) The variation in spe- 
cific mass flow pW,, on the mean stream line follows 
closely the variation in channel width in the pitch di- 
rection inside the channel, the increase in specific mass 
flow is about 4+ per cent higher than that given by a 
one-dimensional calculation based on the reduction in 
channel width (Fig. 2). (2) The shape of the mean 
stream line follows that of the mean channel line (or 
that of the mean blade line obtained by taking the 
mean of the y coordinates of the two blade surfaces) 
but with lower curvature, and the position of the mean 
stream line is closer to the suction surface than to the 
pressure surface (see Fig. 3). Utilizing these rela- 
tions, the flow on the mean stream line can be properly 
estimated as follows: From the blade thickness dis- 
tribution desirable (from 
stress and Mach Number in general and the consider- 
ation of additional requirement of coolant passage in 
the case of cooled turbine blades) in an inverse problem, 
and from that given in the direct problem, the variation 
in specific mass flow on the mean stream line is esti- 
mated from the relation given in Fig. 2 or, even better, 
from a known solution of similar blading. In the in- 
verse problem, the designer may either specify a cer- 
tain desirable shape of the mean blade line and esti- 
mate the shape of the mean stream line according to 
Fig. 2 (or other similar relations) or specify the shape 
of the mean stream line itself or the rate of change of 
W, on the mean stream line (to lead to a certain de- 
sirable loading of the blade). With these data, either 
the coordinates of the mean stream line and the fluid 
state on it may be determined at a finite number of 
, in Fig. 1) or equations may 


the consideration of blade 


stations (such as 2), 2... 
be obtained to represent smooth variations passing 
through the values at these stations. In the direct 
problem, the shape of the mean stream line is estimated 
from the mean blade line according to Fig. 2 or other 


similar relations. 

With this variation of specific mass flow pW’, and the 
mean stream-line shape or its slope tan 8, the density 
and then the velocity components are computed in the 
following manner: Eq. (12) is rewritten as 


> 2 poner os os 1/(y-1) 
W.. :270**W.*? sec? B 


] — 
pg? = — 2H p* (12a) 
pi 1 — (W,2/2H) 
or ; 
2p*? ( _ 
(9*W,* sec 8)? = 1 — ptr (y — 
oe: W.. 2 . oH 
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Comparison between variation of specific mass flow along 
mean streamline and channel-width ratio. 


Fic. 2. 


where the asterisk indicates values made dimensionless 
Either Eq. (12a) is used to pre- 
*W,* sec B 


by the inlet values. 
pare a table of p* for equal intervals of p 
by an iterative process, or Eq. (12b) is used to com- 
pute p*IV,* sec 8 for a number of values p* which are 
then plotted as a graph. After densities have been 
obtained from table or graph, the velocity components 
on the mean stream line are readily computed. When 
the fluid state is completely determined on the mean 
stream line, it can be extended out in the pitch direc- 
tion by the formulas given in the next section (com- 
pare with references 5, 15, and 16). 


VARIATION OF FLUID STATE IN THE PITCH DIRECTION 

Eqs. (8), (10), (11), and (12) directly give the first 
derivatives of W,, W,, and p in the y direction in terms 
of the known variations on the mean stream line as 


follows: 
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Fic. 3. Stream lines obtained in incompressible solution. 
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1 Op l 


=— xX 
p Oy (y — 1) [(H# (W?/2)] 


, OW, _ OW, : 
(1. + W, (15) 
oy oy 


The second-order partial derivatives of I’,, 1, and 
p with respect to y can be obtained as follows: Dif- 
ferentiating the continuity Eq. (1) with respect to y re- 
sults in 


O*(pW.)/Ozd0y| + [0°7(pIV,)/Ov? QO (16) 


Eq. (16) can be written as, through the use of relation 
(9) 


0*(pW,) 0*(pIV,) d Opi.) 
_ _ —— tan 6 + = 0 (17) 
Oy’ Oy dz Ov 
which may be expanded to obtain 
Transposing and combining terms 
OW, | Op d(plV’ ld ( OW 
oy" p’ oy dz p dz Ov 
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OV, (= 1 d(pW,) ) ; P 
- + tan 8} cos* B (13) 
Ov dz p dz 


OV, dW, 1 d(plV,) 
= tan 6B — cos 2B (14) 
Oy dz p dz 


MARCH, 1952 
OW, | , Op OW, 
Oy" ~ Oy Oy 


( ow, 49 Op a) ; ; d O(pW.) . 
“ all {| | 
4 Oy" Oy OY dz Oy 8 


p 


Differentiating the irrotationality Eq. (2) with re. 


spect to y and using the relation (9), 


ow. d OW, ow, 


-_ — tan 6p (19) 
oy’ dz Oy Oy" 
Eq. (10) can be written as 
OVW, OVW, 1 d(pW.) 
= tan 6 - (20 
Oy Ov p dz ; 


in which the last term is only known on the mean 


stream line. Substituting Eqs. (19) and (20) into 


Eq. (18) results in 


o°wW, Op | Me) d O(pW,) 
_9 = 
' Oy" tle: Oy ( p dz t dz oy 


\ OW, 0*W, ‘ | , 
— an an 6 = 
F dz Oy oy? " ’ 


~ 


oO 1 f/oVw, 
W, ") + tan 6 ( )| cos? B 2] 
Oy dz \ ov 


After this is evaluated, the second partial derivative of W, with respect to y is obtained by using Eq. (19). 


O°’, /Oy? = (d/dz) (OW,/dy) — (0°?W,/dy?) tan B (22 


The second-order partial derivative of p with respect to y is obtained by using relations (12) and (15). 


1 0°p 2— vy (/Op\? l OW, OW, OW, \? ow, \? a 
oa" ° wae? i WW 2 , ss WW y 4 “-" ( i, (23 
p Oy" p° Oy (y — 1) [A — (W?/2)] Oy" Ov" Ov Oy 


Third- and higher-order y derivatives, if desired, can be obtained in a similar manner. 


The complete variation 


of any fluid property g across the channel can then be expressed in a Taylor series in (y — y,) using the various 


derivatives evaluated at the mean stream line: 


Og 
gy) = dm + (Y — Ym) “a zo 


(vy — ¥m)? (d°q Y — Vm)? (0% : 
(52) +! ~ (22) #2, (24 
2 Oy] m 3! Ov?) m 


Examples so far obtained in typical turbine cascades indicate that the first three terms in the series are sufficient 


for engineering accuracy in these cases. 


DETERMINATION OF THE BLADE 
INVERSE PROBLEM 


PROFILE IN THE 


The blade profile can be obtained by a consideration 
of mass flow. At a number of z stations, mass flow 
across the y line is computed as a function of y accord- 


ing to the following formula: 


y 
un f pW, dy (25) 


Because the condition on the suction surface is more 


critical, the blade shape on the suction side will be de- 
termined first. From the plot of mass flow against y 
(Fig. 4) at different station z, a number of y,’s are chosen 
for a number of mass flows in the neighborhood of half 
the inlet-mass flow, thus obtaining a number of suc- 
tion surfaces. The corresponding velocities on the suc- 
tion surfaces are read from the velocity plots. A suc- 
tion surface is then chosen for the best velocity dis- 
tribution. After the suction surface is determined, the 
pressure surface is obtained by the total mass flow re- 


quirement, 
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FLOW PAST 

If the shape or thickness of the blade or the velocity 
distribution on the blade obtained is not the one de- 
sired, the shape of, and flow on, the mean stream line 
can be modified accordingly and the process repeated. 
Because each case takes only a small amount of com- 
putation, it is not difficult to build up systematically 
4 series of families of blades for various applica- 
tions 

These data will also be helpful in solving direct prob- 
lems of similar blades by the methods to be discussed 


in the following sections. 


APPLICATION TO THE DIRECT PROBLEM 


The procedure just described can be applied to ob 
tain a quick approximate solution for the direct prob- 
lem. There are only two main differences: One is 
that the exit angle is not definitely known, in contrast 
to a specified value in the inverse problem, and the 
other is that the blade surface obtained must approxi- 
mate closely the given blade and the velocity compo- 
nents be consistent with the slope of the blade con- 
tour. 

For the first one, if the blade has a sharp trailing edge, 
the Kutta hypothesis may be used; but in the case of 
practical blading with a trailing edge of finite thickness, 
the available empirical rules are more appropriate. 
For the second difference, it is necessary to correct the 
solution to fit the given blade shape. <A simple method 
to adjust the value of plV, according to the mass flow is 
suggested as follows: In Fig. 5a is shown a number 
of stations a to g used in obtaining the mass flow across 
the y line which are equally spaced but have different 
values in the two segments. The variations of plV,, 
and mass flow along y (integrating from y,,) obtained 
in the solution are shown by the solid lines in Figs. 5b 
and 5c, respectively. Assume that the blade coordi- 
nates obtained in the solution are 7 and k. Whether 
(yj) — ¥,) 1s longer or shorter than (y, — y,) indicates 
whether pIV’. value is, in general, lower or higher than 
the correct value. Assuming that the variation of 
plV, with respect to y is correct, a constant A may be 
added to the pl’. values obtained in the solution. The 
magnitude of A is then determined by the correct 


amount of mass flow 


Vp 
inlet mass flow — [ plV, dy 


A= 
Vo — Ve 


(26) 


A corresponding increase in mass flow equal to A(y — 4y,) 
is then added to the curve in Fig. 5c, resulting in the 
dashed line. A new position of the mean stream line 
Ym is obtained by reading off this new curve at half the 
given mass flow. 

This process of correction is to be applied at each 
station along the mean stream line. Further modifi- 


cation can be made to get the right flow direction at the 
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STATIONS 





oO 2 4 6 8 10 i2 14 16 |8 


Fic. 4. Variation of mass flow across channel, 


blade surfaces. The mean stream line may be cor- 
rected a short distance outside the channel by using the 
finite-difference equations given in the next section, or 
the end-point differentiation formula may be used at 
the first and last station inside the channel. If accu- 
rate determination of the velocity distribution near the 
leading edge is desired, either shorter spacings can be 
chosen there and differentiation coefficients for unequal 
spacings used, or a portion of the flow region near the 
leading edge be improved by the finite-difference meth- 


ods to be discussed below. 








(a) STATIONS ALONG y AND BLADE 
COORDINATES OBTAINED 
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Fic. 5. Corrections for the direct problem. 
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grid points. 


CHECKING AND IMPROVEMENT OF SOLUTION 
BY FINITE- DIFFERENCE METHOD 


In a new type of blading where the relationship be- 
tween the characteristics of the mean stream line and 
the blade and the number of terms necessary in the 
series to represent the flow variation across the channel 
are not known, it is desirable to check the accuracy of 
the solutions obtained for both direct and inverse prob- 
lems. Because the variation of fluid state throughout 
the domain is available in the solutions, this checking 
and further improvement, if found necessary, of the 
solution can be readily made by some numerical meth- 
ods employing finite-difference approach."~” This 
process also offers a quick way of establishing some 
standard solutions for a number of typical bladings. 

In the numerical methods, it is convenient to solve 
the problem through the use of a stream function, which 
is defined according to the continuity equation [Eq. 
(1)] with 


pW, 
pW, 


= dy/oy \ 
—dy/dz4 


‘ (27) 


SCIENCES—MARCH, 1952 


Substituting Eq. (27) into Eq. (2) results in 
0” 0” oO fe) 
v vy (1, ~ oe ") =0 (2% 
Oz Oy 


Oy? 02? 
Denoting the third term by / and the differentiatioy 
coefficient that multiplies the y value at point j to give 
the second derivative at point ¢ using n” degree poly. 
nomial as 7Bj, the finite difference form of Eq. (28 
at the grid point whose y value is y’ can be written as 


+ 


>~ 3BW¥ +> BY + J' =0 (29 
j=0 k=0 


where y’ and y‘ denote y values along the z and y axis, 
respectively (see Fig. 6). For just checking, the equally 
spaced differentiation formulas*! can be used through. 
out by employing fictitious grid points lying inside the 
blade but equally spaced to the grid points in the chan- 
nel. For further improvement calculations, it is desir- 
able to use the unequally spaced differentiation for. 
mula” to avoid the interpolation of the boundary values 
at these fictitious points and to include the exact value 
at the blade boundary in the calculation in one step; 
this is most desirable in automatic machine computa- 
tions. If the second degree polynomial representation 
is chosen (nm = 2), three points are involved in the nu- 
merical differentiation in each direction, and the center- 
point formula can always be used for all interior grid 
points. If fourth degree polynomial representation is 
chosen (x = 4), five points are involved in each direc- 
tion, and an off-center point formula has to be used at 
the point next to the boundary (Fig. 6). It is there- 
fore desirable to have the spacing near the boundary 
somewhat shorter than the regular equal spacing so that 
accuracy at the point next to boundary is comparable 
to that at other points.” By using a five-point system, 
many less grid points are required in comparison with 
the three-point system for the same accuracy.” In the 
present problem, the variation of y is only rapid near 
the leading and trailing edges, and the grid pattern 
should be determined for these regions first. 

After the degree of polynomial representation and 
the size of grid are decided upon, a finite difference 
equation [Eq. (29)| can easily be written for each in- 
terior grid point using the differentiation coefficients 
for equal and unequal grid spacings given in references 
21 and 20, respectively. The constant y value on the 
suction surface can be arbitrarily chosen, and the con- 
stant y value on the pressure surface is determined by 
the mass flow passing through the channel or at the in- 
let as follows: 

Up oy Up 


dy = 


W.dy = M, (30 
Us oy Us siti 


Vy iain Vv; = 
Outside of the channel, two reference lines, a pitch dis- 
tance apart, can be drawn either parallel to the 2 axis 
or to the inlet and exit angles. (The latter, as shown i 
Fig. 7, is more convenient for the purpose of drawing 
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FLOW PAST CASCADE 


stream lines after the solution is obtained.) Because 
the fluid state repeats itself after each pitch distance 
along the y direction, only the grid points lying between 
the reference lines need to be included in the calcula- 
n, and the central five-point formula can always be 


tio 
For ex- 


ysed at every grid point between these lines. 
ample, at grid point c (Fig. 7) the y values at a and 6 are 
obtained in terms of the y values at 7 and 7 by the fol- 


lowing relation: 


y=v t+ (vy — We), (31) 


Y=W7+(% — vs) 


Inlet and Exit Boundary Condition 


The solution of the cascade flow problem in terms of 
y has an interesting feature in the boundary conditions. 
As we have seen in the previous section, the y value is 
constant along either the suction or pressure surfaces, 
which are fixed boundaries. Outside the channel, there 
is no longer a fixed boundary, but there is the condi- 
tion that the flow repeats itself for every pitch distance. 
Sufficiently far upstream of the blade, an inlet station 
can be chosen so that the boundary conditions there 
require the fluid state to be uniform in the ¥y direction 
and prescribe the flow angle 8;. Two methods have 
been devised to account for this boundary condition as 
follows: 

(1) Fixed-Angle Method. 
flow region near the inlet station 7-7. 
points to the left of 7-7 are not included in the calcu- 
lation.) To write the finite-difference equation at 
point a, the value of y at point c to the left of the first 
station can be obtained by using the given conditions 


Fig. 8 shows a part of the 
(The grid 


: , 6, — 6, tan B; ., 6, tan B; ., 
y=y'= : iin teak Cais die 


6. a 6, b 99 
1 ——tan B,)y + tan B;) py’ (32) 
b, 5, 


These coefficients of y“ and y’ in Eq. (32) are to be 
added to the regular coefficients at points a and ); 
otherwise the points a and } are treated in the same 
manner as other interior points. Whether the inlet 
station i-7 is chosen far enough out will be indicated by 
the linearity of the y variation at that station as ob- 
tained in the solution. 

(2) Stream-Line-Adjustment Method.—lIf the 
station 7-7 is chosen sufficiently far from the blades, the 
variation of stream function to the left of the station 
The value of the stream 
If solu- 


first 


-1 is linear in the y direction. 
function, however, depends on the inlet angle. 
tions for a range of inlet angles are desired, they can be 
obtained by specifying a number of linearly varying 
stream functions to the left of station 7-7 as fixed bound- 
The stream line obtained in the solution 
If, however, 


ary values. 
then gives the value of the inlet angle. 
the solution for a certain inlet angle is desired, the 
Stream line obtained in the solution must be adjusted 
according to that inlet angle (for example, as ad in 
Fig. 9b is adjusted to position ac), thereby obtaining 
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Fic. 8. Treatment of inlet boundary condition 
C [L u 
Z 
& 
= x 
Fic. 9. Zero elements in [C], [ZL], and [U’] 


an improved set of boundary values of the stream func- 
tions to be used in the next calculation. This method is, 
of course, not so accurate and convenient as the fixed- 
angle method for obtaining a solution for a given inlet 
angle but is desirable in the matrix solution (to be dis- 
cussed in the following section) because the inlet angle 
is then not involved in the matrix factorization, making 
the same matrix factors usable for a range of inlet 
angles or Mach Numbers. 


The same methods can be applied at the exit station. 
For a sharp trailing edge, the Kutta condition may be 
used, and the correct exit angle is the one that gives 
the flow at the trailing edge satisfying that condition. 
For round trailing edges, it is more appropriate to use 
some available empirical rules. If the calculation is 
made to compare with certain experimental results, 
either the measured exit angle or the exit angle com- 
puted from the exit pressure measurement may be used. 


Solution of the Finite-Difference Equations 


With the grid system and degree of polynomial rep- 
resentation chosen and the boundary conditions taken 
into account, the problem remaining is the solution of 
the set of N linear algebraic equations [Eq. (29)] writ- 
ten for N interior grid points. For a small number of 
solutions with a given blade, the best method is South- 
well’s Relaxation Method.’ '* A modification of this 
method involving the use of higher-order differences is 


Formulas and the table of co- 


suggested by Fox.” 
efficients obtained in reference 20 enable the direct use 
of higher-order polynomials for problems with curved 
boundaries. For the present flow problems, it is neces- 
sary to include a large domain to get to the boundary 
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conditions that are given at places far away from the 


blades. 
reduces the numerical work. 


If a number of cases are to be solved for a given cas 
cade, it is advantageous to solve the problem on a 
digital computing machine. If a high-speed machine 
is available, the simultaneous equations may be best 
solved by the Liebmann iterative process, which is most 
simple to set up. If only a relatively slow-speed ma 
chine is available, the following matrix process is found 
to be very satisfactory. The set of N equations are 
first put into a matrix form as 


Civ a or /C 


iat _ 29 
Wy yay (+55) 


where |C] is a square matrix and {yw} and {a} are col 


umn matrices, all of order N. Superscript 7 refers 
to the number of the row, and subscript 7 refers to the 
number of the column. In order to keep the special 
feature of the matrix C in having a large number of zero 
elements during the calculation (Fig. 10), it is found 
best to solve Eqs. (82) by factorizing into a lower and 
an upper triangular matrix according to the theorem on 


25 


triangular resolution of matrices.”? 


C] [7] [D] |] (34) 


where |/] and {w~]| are lower and upper triangular ma 
trices with unit elements along the diagonal and where 
|D] is a diagonal matrix. It is convenient to combine 


[7] [D] into a lower triangular matrix |/], resulting in 


[C] [L] [u] (35) 


The elements of |] and [wv] are given by the following 
multiplication rules: 


Operation along a column 7, 















































j—1 
Lf = C} — > L's} (¢ => 9) (36) 
k=1 
Operation along a row /, 
¢—1 
Ci — > L,iu} 
ut = k=t (i < j) (37) 
bg 
~ ¥ civ gizi+ See €'=ity’' 
‘ N we jz J jzitl a 
il 
i<j 
N C ay — 2 Bs] —— Wie 
i=j 
4 ae 
: wi: e' _ c i a, e* is Sei j 
it is i oe ¥ °F Zev 
if, 
C- SEE | 
i | st J ee 
i<) E;- i €*8-TE€ 
C; jpn j 
Fic. 10. Operating diagram and formulas. 
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The use of higher-order polynomials greatly 
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L| and |u|, the 
(33) involves the following forward and backward sy) 


Having obtained solution of E 


stitution process: Defining a column matrix 3 by 


so that 
paa 
a 2. / 3 
1 
J on 
Le 
then, 
ul IW} 1 Bi () 
Hence, 
i+1 
y 3 Ds i y () 
j \ 


These operations are summarized in Fig. 10, in whieh 


the matrices [1] and [uw] are put together into a com 


posite matrix 
| Je] L| + [ul] va 


Also shown are the column matrices | vy}, }6!, and }« 
for checking purposes with 


\ 
y » 3 ( j 1 a } 
j=1 
V+1 
5 i+ DO 1 
j i+ 1 
€ 1+ y' (43 


This method of solution has been successfully carried 
Richard 
the Lewis laboratory, on an IBM Card Programmed 
UNIVAC. On_ the 
CPEC, about 60 machine-hours are needed in the fac 


out for a similar problem by L. Purner, of 


Electronic Calculator and =a 
torization process and 2 hours for the solution of | ¥| 
for a given {a}. On the UNIVAC, the two operations 
take about 11 
calculations, 9 to 11 digits are used, and substitution of 


and 2.5 min., respectively. In_ these 


tw} into Eq. (29) is found to check up to within one in 
the next to last digit. Because | a} is to be determined 
on {p} or {W}, it should be successively improved using 
successively improved {|W} values. The use of these 
automatic machine computations greatly reduces this 
The calculation ol 


ta} from each improved value of |W} is also greatly 


tedious repetitious manual work. 


reduced by setting the machine to compute partial de 
rivatives of y and p with respect to z and y. 


SPECIAL CASE OF INCOMPRESSIBLE FLOW 


For incompressible flow, all the previous calculations 
are greatly simplified, especially in the finite-difference 
solution of the direct problem, which now does not re- 
quire the successive correction in the density term, and 
the solution can be obtained in a single cycle of com- 
putation. 
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Number 0.42 


The first and second derivatives at the mean stream 


line are now reduced to simply 


fa! (— dw 
) tan B ) cos” Bp, 
Ov / » dz SS Ja 
(5) (= dw 
| tan 3) cos” B,, 
Oy J» dz dz = 
® W, d (+) d () ; 
) tan 0,, cos” PB», 
Ov? /», dz Ov dz Ov ” 


(46) 


tan a] (47) 


(44) 


( 1) 


(<) P | (~*) ow, 
Oy? J» dz \ Ov Ov" 


ILLUSTRATIVE EXAMPLES 


Direct Solutions 


1) Finite-Difference Solutions. —A highly combined 
thick blade cascade as shown in Fig. | is chosen for this 
analysis. The distribution ob 
tained will be compared later on with the experimental 
data obtained by H. W. Plohr, of the Lewis laboratory. 
Solutions for two inlet Mach Numbers, 0 and 0.42, are 
obtained by the finite-difference method. The inlet 
18’ and —52° 57’, respectively, 


theoretical velocity 


and exit angles are 41 
conforming to the experimental value. The ratio of 
pitch to axial chord is 1.017 to 1.5. The ratio of [7 to 
W,, ° is 25.78 in the test. y is taken as 1.4 in order to 
compare with experimental data obtained with cold 
ui. 

The grid chosen is shown in Fig. 7. The finite- 
difference equations are obtained by using five-point 
coefficients given in reference 20. The system of equa- 
tions is solved by the relaxation technique. Because 
of the density variation involved, the accuracy ob- 
tained in the compressible flow is lower than that in the 


incompressible case. 
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The stream lines obtained for the two cases are shown 
in Figs. 3 and 11, respectively. The stream lines are 
seen to be pushed further away from the suction sur 
face in the compressible case, and the compressible 
mean stream line also has a curvature closer to that of 
the mean channel line than the incompressible mean 
stream line does. In this comparison, the shape of the 
mean channel line or the mean line of the blade ob 
tained by taking the mean of the blade y coordinates is 
found to be better than the usual mean camber line of 


the blade. 


The variations of specific mass flow pW, on the mean 
stream line are compared with the variation in the 
channel-width ratio in Fig. 2. It is interesting to see 
that the former follows the trend of the latter. Inside 
the channel, the former is about 4 per cent higher. In 
the neighborhood of the leading and trailing edges, 


higher difference is noticed. 


The variations of the velocities across the channel are 


shown in Figs. 12 to 17 for the two cases. For clarity, 
only about half of the stations used in the computation 
are shown. In general, it is seen that the velocity, 
inside the channel, varies faster than a linear rate from 
pressure surface to suction surface and that the vari 
ation is larger in the compressible case. Analysis of 
these variations shows that they can be closely approxi 
mated by second-degree polynomials 


(2) Vean-Stream-Line’’ Solution Solutions are 
also obtained for the same cascade and the same inlet 
and exit conditions by the approximate method utiliz 
ing the The 


stream line obtained in the incompressible solution by 


mean-stream-line’’ technique mean 


the relaxation method, as shown in Fig. | or 11, is used 


for both incompressible and compressible flows. For 
the incompressible flow, the blade coordinates obtained 


in the first solution is sufficiently close to the given 
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Number = 0.42. position 


in the approximate solution of direct problem 





blad 
velo 
solu 
one 
ties 
tion 
seri 
tion 
of s 
indi 
puts 
T 
inle 
com 
Fig. 
pres 
disc 
mer 
and 
tion 
mat 
trib 
blac 


r 
in 

Nui 
vari 
The 


exct 


Inve 


ver! 
line 
thic 
of t 


tail 

wit! 
for 

atin 
The 
alot 
den 
Fig 
Fig 
to « 
den 


cal 


san 
file 





——_ 


LUTION 
AXATION 
ROXIMATE 


’ ACTOss 


—. 


JTION 
CATION 
XI MATE 


CTOSS 


ES 
Ss 6 


ine 





FLOW PAST 
blade, and no further adjustment is necessary. The 
yelocities obtained are compared with the relaxation 
solution in Figs. 12 to 14. For the compressible case, 
one adjustment is necessary (Fig. 18), and the veloci- 
ties obtained are compared with the relaxation solu- 
tion in Figs. 15 to 17. Only three terms in the Taylor 
series are used in these calculations, and only six sta- 
tions are used inside the channel. This small number 
of stations was chosen for the purpose of giving some 
indication of the accuracy obtainable with a short com- 
putation. 

The velocities on the blade surface obtained for the 
inlet Mach Number of 0.42 by these two methods are 
compared with the available experimental value in 
Fig. 19. It is seen that the agreement is better on the 
pressure surface than on the suction surface, and the 
discrepancy involved is within the accuracy of the nu- 
merical calculations and experimental measurement 
and the neglect of viscosity in the theoretical calcula- 
tion. The results also indicate that the quick approxi- 
mate method is adequate to determine the velocity dis- 
work in cooled turbine 


tribution for heat-transfer 


blades. 


Fig. 20 shows contours of Mach Numbers obtained 
in the compressible solution. With an inlet Mach 
Number of 0.42 and exit Mach Number of 0.55, the 
variation of Mach Number in the channel is large. 
The maximum Mach Number on the suction surface 


exceeds 0.8. 


Inverse Solution 

In order to test the accuracy obtainable by the in- 
verse method suggested in this paper, the mean stream 
line obtained in the incompressible solution and a blade 
thickness distribution (Fig. 20) corresponding to that 
of the same blade used in the direct solution are taken 
as the prescribed value. Blade shapes are then ob- 
tained for incompressible flow and compressible flow 
with inlet Mach Number of 0.42. The computation 
for the mean stream line and first and second y deriv- 
atives for the latter are completely shown in Table 1. 
The computation of velocity, density, and mass flow 
along y are shown in Table 2 for one station. The 
densities on the mean stream line are obtained by using 
Fig. 22, which is constructed according to Eq. (12b). 
Fig. 23, which is plotted according to Eq. (12), is used 
to obtain p elsewhere. For convenience, the velocity, 
density, and mass flow are nondimensionalized in the 


calculation as follows: 


p* = p/p, M* = M/p,W,, iP 


Because of the same mean stream-line shape and the 
same blade thickness distribution used, the blade pro- 
file obtained for incompressible and compressible flow, 
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Fic. 19. 
in the direct solutions with experimental data. 





Inlet Mach Number = 0.42. 


Mach Number contour. 


Fic. 20. 
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Fic. 21. Blade thickness distribution 

Fig. 24, look the same, except that the suction surface 
of the compressible case is situated farther from the 
given mean stream line, which is consistent with the 
results previously obtained in the direct solutions. 
Both blades have thickness distributions close to the 
velocities obtained 


prescribed values. However, the 


in the compressible solution are, in general, higher than 
those in the incompressible case (Fig. 25). This is 
mainly due to the higher velocity on the mean stream 
line resulting from the use of the same pV, but with a 


decreasing p in the compressible case. 
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CONCLUDING REMARKS 


These examples serve to illustrate the simplicity in 
the present approach to the solution of the direct and 
inverse problems of the compressible flow past cascade 
of airfoils. For blades with high solidity, such as the 
one investigated, three terms in the series are sufficient 
for engineering accuracy. More accurate results for 
blade profiles and velocity distributions than those 
given can be obtained by using more stations than the 
minimum six used in the examples. This approach is 
particularly suitable for the design of turbine blades 
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TABLE 2 
Computation of 1/* Along y at Station 25 (Ay 
0.1316) 


= 0.0937, Ayp = 
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Fic. 25. Velocities on blade surface and mean stream line ob- 
tained in the inverse solutions 


when a certain thickness distribution is essential from 
the requirements of blade stress and of coolant passage 


in the case of a cooled turbine. 
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Besides the direct application to two-dimensional 
flow on cylindrical surfaces in turbomachines, the 
method can be extended to two-dimensional flow on an 
arbitrary surface of revolution in radial or mixed-flow 
turbomachines. Families of cascades of blades for 
these machines can be systematically built up by this 
method in a relatively simple way. The method can 
also be used in the design of channels in a plane or on 
an arbitrary surface of revolution. 
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Compressibility Corrections for Bodies 
of Revolution 


BORIS A. RABINEAU* 
Douglas Aircraft Company, Inc. 


SUMMARY 


An approximate formula for a compressibility correction on 
bodies of revolution at 0° angle of attack is derived. The maxi- 
mum perturbation velocity calculated by this formula on a 12 
per cent thick body at a Mach Number of 0.9 agrees within 0.5 
per cent with the results of more elaborate calculations based on 
potential theory and Goethert’s three-dimensional rule. The 
perturbation-velocity peaks are located in this case at 10 per cent 
of the body length from the body ends and increase 76 per cent 
between 171 = 0 and 0.9. The same formula applied to the ex- 
perimental data at a Mach Number of 0.9 on a 4.5 per cent thick 
body with elliptical ends (each end having a thickness ratio of 15 
per cent) gives the maximum perturbation velocity (at 7.5 per 
cent of the body length from the nose) with an error of 0.08 to 
0.15 per cent of the free-stream velocity and about 2 to 3 per 


cent of the perturbation velocity. 


(1) INTRODUCTION 


; es INFLUENCE OF COMPRESSIBILITY in axially sym- 
metric flow around slender bodies has been the 
subject of a number of papers in recent years. At 
present, the most reliable method appears to be the 
Prandtl-Goethert rule,! which requires (a) the reduc- 
tion of the body transverse dimensions in the ratio 
v1 — M?*, M being the free-stream Mach Number; 
(b) the calculation of the perturbation velocities Au,/u 
in incompressible flow around the reduced body; and 
(c) the multiplication of these velocities by 1/(1 — M?). 


It is essential to note that velocities around affine 
bodies of revolution in incompressible flow are not 
similar. In particular, the points of zero and maxti- 
mum perturbation velocity are at different locations on 
the body axis and the maximum velocities are not pro- 
portional to any power of the thickness ratios. Hence, 
the correct application of Goethert’s rule requires a new 
calculation for each Mach Number. Accordingly, a 
formula was sought, giving the influence of compres- 
sibility for any Mach Number, which would take into 
account the general body shape and would, in particu- 
lar, be valid for a general and practical class of bodies 
of revolution having cylindrical central portions. The 
application of most of the published compressibility 
correction formulas to these bodies is generally not 
satisfactory. The peaks of velocity, located near the 
body ends, increase much more rapidly than the ve- 
locities near the center of the body, and the compres- 


Received June 27, 1951. 
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sibility influence is also apparent on the cylindrical por- 
tion where formulas based on the local slope and curva- 
ture give identically zero. 


(2) PROPOSED FORMULA 


The proposed formula is as follows :t 


Au, Au, 7 - 
u u 2(1 — p?) 


a ~ ?(1— M? 
log (1 — M)(DA,Ps'cn) 1 — 4 (1) 
1 2(1 — yp?) 


where Au./u is the ratio of the perturbation velocity at 
any point of the body surface to the free-stream velocity 
at the free-stream Mach Number, /; Au,/u is the same 
ratio for the original body in incompressible flow; M is 
the free-stream Mach Number; ¢ is the body thickness 
ratio = maximum diameter/total length, 2/; uw is the 
distance from the point where the velocity is deter- 
mined to the body midstation, expressed in body half 
—1 at the nose, 0 at body midpoint, +1 
and 


lengths |u = 
at body tail (see Fig. 1 for notation and axes) |; 
n 

>> is the summation of the infinite series A,P,'Cp, 
I 

where 


+ See derivation in the Appendix and Theory of reference 2. 
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2n l saath 1 
An = .. : 2” na P,(u) du), 
M 


THE 


4 =-1 l dz 


r being the body radius, z the distance to the r axis, the 
origin being at the body midpoint, and n an integer, 
taking all the values from 1 to ~. P, is a Legendre 
polynomial of the first kind of order n in u;* P,,’ 


derivative of P,, with respect to uw; c, is defined as the 


is the 


average value of the expression 

Cn = dQ,/d(log t?) 
in a certain interval /,, tf, of thickness ratios to be de- 
fined below. 


kind of order nin = V1 + 22.3 
Eq. (1) was derived by applying Goethert'’s rule, as ex- 


Q, is a Legendre polynomial of the second 


plained above, to the expression of the perturbation 
velocity in reference 2 


Au 1 — p? i 
uo ( +? — p»? m 
(vi +P-x>Dd AnPx'Qu ) —~1 (2) 
l 


All the terms of order higher than /? were neglected in 
this formula both in the compressible and the incom- 


pressible case. (The corresponding body thicknesses in 
incompressible flow were taken as fy =/, and fy V1 — M? 
= fy.) It wasalso assumed that 

QO, ~ Cn log t? + d, (3) 


c, and d, being constant in the interval 4, fz (see Appen- 
dix). 

The values of c, for n varying from 1 to 7 have been 
tabulated in Table 1 for the intervals ¢; = 0.04 to ft = 
0.10 and 4; = 0.10 to 4, = 0.20 corresponding to the 
most usual thickness ratios. 


(3) APPLICATION OF THE FORMULA 
Some remarks regarding the determination of the 
n n 
sums )}> A,P,’Q, and >> A,P,'c, fora given body 
I I 


might be in order. 

(1) If no experimental velocity distributions are 
available, the method of reference 2 can be applied 
using the tabulated values of P,,, P,,’, Q,, and u and the 
values of c, of Table | of this paper. 

(2) If experimental velocity distributions at two 
Mach Numbers /, and W, (one of the two can be zero), 
or for two thickness ratios ¢; and f: at J = O are avail- 
able, the calculation can be considerably simplified. 


Values of c, for ¢ = 


AERONAUTICAL 


TABLE | 


SCIENCES—MARCH, 1952 


Eq. (2), applied to the equivalent incompressible ve. 


locity (Au/u) (1 — M?) and thickness ratio tf) V1] — M 
gives the values of z. Aegt 'O. corresponding to the 
1 


thicknesses /; = ty V1 — Mand tf, = tp V1 — yp 
Since it is evident from the definition of A, that they 
vary as the square of the body thickness ratio /, whereas 
P, and P,,’ are independent of it, we can write for the 
original body (subscript 0) 


, 
A ee _ 


since A,, corresponds to the affine body of thickness 


ratiot; = f)V 1 — M,?. Also, 


Qn, * Cy log to? + d, = ©, log (to?/t?) th? + 
d, = Cy, log th? + d, + Cy log (t/t; 


Qn = Qu, + Cn log [1/(1 — M,2)] 5 


Hence, 


n n 


Awl n'On = = A, F'Gu 
u 0 ¢ ( (1 Th & l Qn, 


log (1 — My?) A,,P,'c, (6 
1 


Similarly, 


n l n 
4 . m — 9 44h ny . Pas pen 


log (1 — AM,?) >. Aa FP by. He 
i 


These which 2YA,,P,’Q,, and 
YA,,P,’Q, are determined by application of Eq. (2 
with known velocity values Au(1))/u and Au(M,)/u at 


two equations, in 


corresponding points (same values of u!), give the value 
of the = term appearing in Eq. (1) without further com- 
putations. Eqs. (6) and (7) assume that the same 
value of c, applies in the JJ; — 0 and J/, — O inter- 
vals; this is sufficiently exact if the original body thick- 
ness ratio does not exceed 20 per cent and the ends are 
not too blunt. 
(4) COMPARISON OF RESULTS ON A 12 PER CENT 
THICK Bopy 

In order to have an idea of the error involved in the 
application of Eq. (1) instead of the complete calcula- 
tion outlined in reference 2, both were applied to the 
case of a 12 per cent thick body of revolution at a Mach 


0.04 tot = 0.20 


n l 2 3 } o 6 ( 
Group I: ¢ = 0.20 to 0.10 

tn = —0.47243 —0.443806 —(0.41076 —0.37774 —(). 34534 —().31429 —() 28504 
Group II: t = 0.10 to 0.04 

stig — 0.45837 —0.44452 —0.43009 —(0.41531 


—0.49249 —0.48286 —0.47130 
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The body was assumed to have a cen- 


Number of 0.9. 
tral cylindrical mid-section of 50 per cent of the body 


total length and spheroidal nose and tail (see Fig. 1). 


The angle of attack was Q”.. 


The peak velocities on such a body are located 
approximately at 10 and 90 per cent of the body 
length, as shown. They have the following values: 
at M = 0, 


(Au/U) mar. = 0.0710 


at M = 0.9, 
(Au/U) mar. = 0.1242 
as calculated by the method of reference 2 and Goeth- 
ert’s rule. 
Applying Eq. (1) and calculating the two terms with 
the group of c, values corresponding to the interval 
t, = 0.10 to fg = 0.20, we find, at M = 0.9, 


(Au/t) mar. = 0.1193 (for nh = 0.8) 


The error in the peak velocity is therefore 0.0049— 
ie., less than 0.5 per cent of the free-stream velocity. 
This is of the same order as the approximation of the 
theory of reference 2 and can be considered as satis- 


factory. 


(5) COMPARISON OF THE COMPRESSIBILITY 


CORRECTION FORMULA WITH EXPERIMENTAL DaTA 


The formula and the procedure explained above 
have also been applied to one of the bodies tested in 
reference 4. 

The body having the nose fineness ratio of 15 per cent 
was selected because it seemed to have the best experi- 
mental points. 

The perturbation velocities at Mach Numbers M = 
0.35 to 0.7 at a station (No. 4) located approximately 
7.5 per cent of the body length from the nose were used 
to calculate the two > terms from Eqs. (6) and (7) 
(this Station 4 is the nearest to the maximum perturba- 
tion velocity point on the body). 

The perturbation velocity at Station 4 at zero Mach 
Number was then computed from Eq. (2); the value 
found was Au/u = 0.0248. 

Finally, the perturbation velocity at the same sta- 
tion at a Mach Number of 0.9 was calculated using Eq. 
|), and the value of the A,,P,,’c, term was deduced from 


Eqs. (6) and (7). The result of the calculation was 


(Au/U) 4 =0.9 = 0.03875* 
This compares to the measured value, corrected for the 
*The values of Au/u, local perturbation velocity at Mach 


Number 1/, were deduced from the isentropic pressure coefficients 


given in Fig. 18 of reference 4 by the formula 


jy. 2 j ¥ 1 ' , Au \? &y (y-D ) 
Ga. met ae cE (: ail 


where 4 1.4 


wall influence and read from Fig. 18 (Station 4) of refer- 
ence 4 


(Au/u) 4-99 = 0.03960 


The absolute error therefore, was, in this case, 0.00085 
u—t.e., 0.085 per cent of the free-stream velocity and 
about 2 per cent of the value of the maximum per- 
turbation velocity Au/u. In view of the high Mach 
Number and the magnitude of wind-tunnel corrections, 
this can be considered as extremely satisfactory. 

Taking the perturbation velocities on the same body 
at Station 5, on the other side of the velocity maxi- 
mum, for the same Mach Number J = 0.9, the fol- 
lowing values were found: perturbation velocity at 
M = 0, 


Au,/u = 0.02806 
calculated perturbation velocity at 1 = 0.9, 
(Au/u) y=0.9 = 0.0422 
Against the value read on Fig. 18 of reference 4 
(Au/u)y=o9 = 0.0435 
The error in this case is 0.0013u—1.e., less than 0.15 
per cent of the free-stream velocity, which is still well 
below the approximation of the theory of reference 2. 
The formula proposed in Eq. (1) therefore appears 
satisfactory in the case of slender bodies of revolution 
with velocity peaks near the ends, even at very high 
subsonic Mach Numbers. 


BODIES WITH APPRECIABLE CONTOUR 


SLOPES 


(6) CASE OF 


Eq. (1) was derived by assuming the slopes s of the 
body contour small everywhere and the radial com- 
ponents of the perturbation velocity s(u + Au) neg- 
ligible (see reference 2, page 9). 

If this is not the case, it can be easily deduced from 
the Appendix that the total local velocities W,,, W,, at 
the same point of the body at Mach Numbers M and 


zero are related by 


W1. all ty M? | to?M? 
“= + — | — — — 
u u 2(1 — pu?) 2(1 — yp’) 


= (1 — M?)ty? 
[log (1 — wy( APr'a)| — = *] (9) 
1 2(1 — yp?) 


the perturbation velocity x component, Au, being now 


related to W by 


W? = (u + Au,)? + AV,? = (u + Au,)? + 


s?(u + Au,)? = (1 + s*) (u + Au,)*? (10) 


The pressure coefficient C, is given by Eq. (8) of Sec- 
tion (5) (footnote); it is easy, therefore, to improve the 
approximation of Eq. (1) near the body ends, if neces- 
sary, by using Eq. (9) in conjunction with Eqs. (8) and 
(10) to find W,, and C,, from W,,. 
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Appendix 


Derivation of Eq. (1) 


Eq. (2) taken from reference 2 and Goethert’s rule give, as perturbation velocity on the ‘“‘reduced’”’ body in in. 


compressible flow, * 
Au Au 1 — p? ve L 
= 1 — M?) = ; 1+ — > A,*P,’0,® | — 
E ), u , E eer | rs = . | | 


Q, = cn log #? + d,; tp? = to? (1 — M?) 


s 


Assuming 


we obtain 
Qn" = Cc, log f? (1 — M?) + d, = c log (1 — M*) + Q,° 
and since 
A,*® = (1 — M”) A, 


>A P.Oe = (1 wi M?*) [24 Re hy tA, — log (1 = M")As'Px'¢n | 


Ate _ ur?) -( Seal 4 ] x 
Uu ; i l1— pe + ty?(1 = M?) 


fo2(1 — M?) 
! += — (1 = MP) [EX Ay?P,/Qu? + log (1 — MY) AntPa'cy i i 


A i—s /* fo? 
a | ae x - — LA,"Py'Qu? — log (1 — MY) A,*P,'cyh . 


Uu L— pw? + (1 = 2 
dl ‘| 1 — g° we _ i 
1— M*tL1 — pw? + 421 — M) 
Developing the uv factor in powers of fo? assumed small 


| l— 2 y x [ fy2(1 — M?) ] oe : fo?%(1 — M?) 
1 — w+ &2(1 — M?) 1 — pw? + &? (1 — M?) 2(1 — w? + 4? (1 — M?)] 


Hence, 


| if a cd - be | 2 it 7 <i - — bh? 
1— MALI — p? + &? (1 — M?) (ot — pst +49 (1 — 1 1 — p?) 


since f9?(1 — M*) is much smaller than | — yu? in the range of validity of the theory. 
Also, 


1 yw? Mp _— ” ” 
; = A, P,’ ),° = l = A Bad oy ),°-1 = 
ee —e +H (1 - ey z ~ ‘ | F — p+ (1 — zit | si 2 ey ‘ 


1 — p* + 2 - — bt 
l — )A,° P,'O, 
=u wi] aa a] [+ 5 - Eas Pv] - 
P — p? + hh? ( - it (1 — M? iz [1 + wea i + a] J 
l — p? + hh? (1 = “3 2(1 — p’) Uu 
[ _ to? l = | = [ 4 ty? M? | rs ty? 
21 — -p?) - 2(1 — p?) u 2(1 — p’) 


Au;/u being the perturbation velocity on the original body in incompressible flow. 
Hence, finally, Eq. 1 


Au, to? Au, lo? 
= E mr =| 2 en [log (1 — Tey pe A,Py’c]f1 — (1 — M’) | Q.E.D. 
u 211 — w*)J u , 2(1 — yu?) ‘ 





* Superscript or subscript R, the corresponding superscripts on the original body in incompressible flow being 0 


(List of References for This Article Can Be Found on Page 206) 
J g 
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A Preliminary Study of Reynolds Number 
Effects on Base Pressure at .J/ = 2.95’ 


SEYMOUR M. BOGDONOFFYt 


Princeton Unwersity 


SUMMARY 


A preliminary study was carried out to determine the effect of 
Reynolds Number on the base pressure of a simple cone-cylinder 
Tests were made at WJ = 2.95 over a range 


body of revolution. 
The body 


of Reynolds Numbers from 0.6 X 108 to 18 & 10®. 
was tested in the smooth condition and with transition fixed by 


several different roughness bands. For the results obtained see 


the section on “Conclusions.” 


INTRODUCTION 


i ben PROBLEM OF PREDICTING the pressure on the 
blunt base of a body of revolution traveling at 
supersonic speeds has not been susceptible to theoretical 
treatment. Several investigators have made experi- 
mental studies of the problem in supersonic wind tun- 
nels and, more recently, by free-flight tests of missiles. 
The results obtained have, however, severe limitations 
that greatly restrict the application to more general 
problems or other configurations. The wind-tunnel 
tests are limited in many cases to Reynolds Numbers 
below about 4 X 10°. Some higher Reynolds Number 
tests (over 10 X 10°) have been carried out at the 
lower supersonic speeds, J = 1.5 and 2.0. Few of the 
supersonic tunnels have any appreciable Reynolds 
Number range. In some cases, an attempt to obtain 
the effect of variable Reynolds Number has been made 
by varying the length of the test body. This increase 
in body length introduces two variables that are not 
directly part of the Reynolds Number effect: (1) the 
pressure distribution over the body changes and (2) the 
boundary-layer thickness increases as the effective 
Reynolds Number increases. For a given body, in- 
creasing the Reynolds Number by increasing the den- 
sity or the model*size decreases the thickness. Both of 
these effects may have a considerable influence on the 
base pressure. Actual missiles or airplanes operate 
over a wide range of Reynolds Numbers, 1-80 X 10° 
and Mach Numbers from 1 to 6 or 7. It is perhaps 
possible to cover this range by free-flight tests, but the 
results are confined almost entirely to bodies that are 


Received August 13, 1951. 

*The present study is part of a program of theoretical and 
experimental research on viscous effects in supersonic flow being 
conducted by the Department of Aeronautical Engineering, 
Princeton University. This research is sponsored jointly by the 
Office of Naval Research (U.S.N.), Mechanics Branch, Mathe- 
matical Sciences Division, and by the Office of Air Research 
(U.S.A.F.) 
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fin-stabilized. The fins and the details of the fin-body 
function and position may have a considerable effect on 


the base pressure. 


There are no wind tunnels presently available which 
can cover the entire flight Reynolds Number range at 
any given Mach Number. The standard technique is 
to roughen the model to get high effective Reynolds 
Number for low test Reynolds Number. The question 
thus arises as to what is the value of the effective Rey- 
nolds Number—i.e., what is the free-flight Reynolds 
Number that corresponds to the roughened model test? 
There are several general methods of applying the 
roughness or transition strip, but there is no way of pre- 
dicting the corresponding flight Reynolds Number. 
If the characteristic being studied reaches a constant 
value at high Reynolds Numbers, then the problem is 
simplified, provided the method of producing the transi- 
tion or roughness band has no effect on the final value 
observed. 


The only available results over a wide Reynolds 
Number range are those shown by Chapman, of the 
N.A.C.A., in reference 1. Using two wind tunnels, 
Chapman was able to get tests from a Reynolds Num- 
ber of about 0.4 & 10° to 10 X 10° at a Mach Number of 
2.0. The results show a sharp decrease in base pres- 
sures as the Reynolds Number is increased from the 
lowest values, a leveling off, and then a sudden increase 
in base pressure at about Re = 5 X 10°, with approxi- 
mately uniform pressures up to 10 X 10°. The schlie- 
ren photographs verify that the flow over the rear of the 
body is laminar for the initial phase and turbulent for 
the final phase. The base pressure increases slightly 
with increasing Reynolds Number once the turbulent 
region is established. The turbulent results at J = 
1.5 show a greater increase than those shown at J = 
2.0. 


In an attempt to establish the transition region and 
to study more completely the results at high Reynolds 
Numbers, a series of tests has been carried out in the 
Princeton 4- by 8-In. Variable Density Blow-Down 
Supersonic Tunnel. A simple cone-cylinder body was 
studied at a Mach Number of approximately 3 over a 
Reynolds Number range from about 600,000 to 15,- 
000,000, based on body length. The effect of several 
transition bands was also studied in an attempt to 
evaluate the effective Reynolds Number of a body with 


fixed transition. 
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Fic. 1. Mounting and moving mechanism for the base pressure 
model tests. 
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Fic. 2. Geometry of cone-cylinder test body and support. 


NOTATION 


D = body diameter, in. 

L = body length, in. 

Re = Reynolds Number based on body length 

P, = base pressure, Ibs. per sq.ft. 

P, = static pressure in the free flow ahead of the model, Ibs 
per sq.ft. 

M = Mach Number (free stream) 


EXPERIMENTAL APPARATUS 


(a) Supersonic Tunnel 


A complete description of the design and operation 
of the 4- by 8-in. tunnel used in these experiments will 
be available in a report to be issued in the near future. 
For completeness, however, a brief description follows: 
The tunnel is of the blow-down type using a tank sys- 
tem of 1,000 cu.ft. at 250 atmospheres as the driving 
power. A regulator system between the tunnel and 
the tanks allows the stagnation pressure in the tunnel 
to be varied at will between approximately 50 Ibs. per 
sq.in. absolute (minimum operating pressure for the 
M = 8 nozzle used) and 500 Ibs. per sq.in. absolute. 
The settling chamber is provided with three 40-mesh 
screens, and there is a contraction ratio of approxi- 
mately 80 to the nozzle throat. Stagnation tempera- 
ture and pressure are recorded at the same time as data 
on the model are taken. The test section is + in. wide 
by S in. high and is provided with 9-in. diameter win- 
dows 2'/2 in. thick. A conventional off-axis schlieren 
system using a |-microsec. Charters-Lebess spark at 
15,000 volts was used to take the schlieren photographs 
within a few seconds of the other data recording. 

The tests consisted of runs of several minutes, during 
which time the stagnation pressure was set; two sets of 


data were recorded. This procedure was repeated yp. 
til the stagnation pressure could no longer be held cop. 
stant because of low pressure in the storage tanks 
The check of the two consecutive sets of data at the 
same stagnation pressure showed whether the tunne| 
and instrumentation had reached equilibrium cond. 
tions. If the tests did not check, reruns were made 
The tunnel is not provided with a heater so that there js 
no control of the stagnation temperature. During , 
typical run the stagnation temperature varied approxi. 
mately 10°-15°F. per min. 

Calibration of the tunnel showed a variation of about 
| per cent in Mach Number along the tunnel centerline 
over the length of approximately 8 in. in which the 
model was tested. Complete surveys will be given jn 
the wind-tunnel report noted previously. The tunnel 
also showed a shift in Mach Number due to stagnation 
pressure changes, varying from approximately 2.94 at 
low pressure to 2.96 at high pressures. The mean valye 
over the 8 in. was used for each pressure in the calcula. 
tions, but, since the variation was small, all results are 
labeled as M = 2.95. 


(b) Models 


The models were supported on a sting carried in a 
central body that spanned the tunnel, Fig. 1. During 
starting, the larger models were retracted into a sup- 
porting strut downstream of the test section. Once the 
flow was established, the model was moved forward 
into the test section free from any effects of the down- 
stream supports. All models and sting systems were 
constructed to the geometry shown in Fig. 2. Three 
models were used to cover the range of Reynolds Num- 
bers with considerable overlap. The models were of 
0.25-, 0.50-, and 1.0-in. diameter and were all con- 
structed of paper laminate Micarta. This was done to 
eliminate any heat-transfer effects to or from the model 
due to the drift in stagnation temperature. The tests, 
therefore, correspond to equilibrium conditions over the 
body. The models were finished in a high-speed pre- 
cision lathe using a sharp tool, and no attempt was 
made to polish the surface after this operation. 
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Fic. 3. The variation of base pressure coefficient with Reynolds 
Number for a smooth model 
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REYNOLDS NUMBER 


The base pressure was obtained by drilling four holes 
in the base of the two larger models half way between 
the sting surface and the outside surface of the model. 
The holes were connected inside the model and a single 
sassage led out through the hollow sting support. For 
the 0.25-in. diameter model, only two holes were used 
yecause of the size limitation. On the 1.00-in. diam- 
eter model sting, a small collar that pressed against the 
model base was used to aid in model-sting alignment, 
but tests without the collar showed no effect on the base 


pressure. 
EXPERIMENTAL RESULTS AND DISCUSSION 


(a) Smooth Model Tests 


The results of tests of the three models, all geometri- 
cally similar and in a smooth condition, are shown in 
Fig. 3. The base pressure is presented in coefficient 
form as the ratio of the measured base pressure to the 
free-stream static pressure ahead of the model. The 
general shape of the curve below a Reynolds Number of 
approximately 6 X 10° is similar to that obtained in the 


only other comparable study—-Chapman’s results at 





|M = 2.0. Representative schlieren photographs over 
the Reynolds Number range are shown in Figs. 4a and 
| 4b. 


In the Reynolds Number range from 0.6 to 1.0 X 10°, 
the base pressure dropped sharply in a manner similar 
toChapman's results. The schlieren photographs show 
a typical ‘‘laminar’’ wake region downstream of the 
base characterized by the sharp free boundary layer 
and the absence of any trailing shock for a considerable 
distance from the wake pattern. There must be a com- 
pression fan arising from the turning of the wake flow 
along the sting, but these waves do not coalesce to form 
ashock for some distance. As the Reynolds Number is 
increased, the base pressure decreases rapidly to a mini- 
mum near 1.0 X 10°, while the schlieren photographs 
show a shortening of the wake region and a moving in 
f the trailing shock, although the free boundary layer 
stillseems to be laminar. A tentative explanation, diffi- 
cult to substantiate from the recorded data or schlieren 
photographs, is that transition from a laminar to tur- 
bulent free jet is occurring at the wake-sting junction. 
The turbulent free boundary layer can probably sustain 
ahigher pressure ratio than the laminar layer and so 
gives a lower base pressure. At Reynolds Numbers 
above 1.0 K 10°, the flow has developed the character- 
stic “‘turbulent’’ formation with the free jet showing 
definite ‘‘fuzziness’’ of the free boundary layer and a 
strong trailing shock approaching close to the edge of 
thewake pattern. Further increases in Reynolds Num- 
ber make no significant changes in the wake pattern 
even though the base pressure increases sharply at about 
3 X 10° and then decreases slightly at high Reynolds 
Numbers. The region from 1-3 X 10° is probably char- 
acterized by transition of the free laminar layer to tur- 
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R, = 4-18 x 10° 


Schlieren photographs of smooth base pressure model, 


Fic. 4a. 
2.95, for a range of Reynolds Numbers 


D = 0.25in., M = 


bulent in the wake region before the occurrence of the 
trailing shock. The region above 3 X 10° is tentatively 
associated with transition to turbulent boundary layer 
on the body moving forward as the Reynolds Number 


increases. Theoretical studies of this phenomena in 
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No Flow 


Fic. 4b. Schlieren photographs of smooth base pressure model, 
D = 1.00 in., 7 = 2.95, for a range of Reynolds Numbers 
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Fic. 5. The variation of base pressure coefficient with Reynolds 


Number for several transition strip thicknesses 
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two, as well as three, dimensions are now under way at 
the Supersonics Laboratory. 

Chapman observed an almost uniform base pressure 
region from Re = 2-5 X 10° before the sharp increase. 
but nothing comparable was obtained from these tests, 
Considerable difficulty was experienced in getting con. 
sistent data in the transition region from Re = 19 
3.5 X 108, although the model seemed to be in good cop. 
dition. No attempt was made to polish the model be. 
fore each test. Most of the tests fell between the dotted 
curves shown in Fig. 3, with the few that fell outside 
this region shown on the figure. A large percentage of 
the tests clustered around the solid line shown on the 
figure, but only representative tests are shown to indi- 
cate the degree of the scattering. This transition re 
gion would naturally be influenced by body finish and 
air-stream turbulence level, and, therefore, direct com 
parison of the results with those of Chapman is not 
completely valid. 

In the high Reynolds Number region, where presum- 
ably the boundary-layer flow over the body is primarily 
turbulent, a small but steady decrease in base pressure is 
observed beyond a Reynolds Number of 6 X = 10¢ 
This is again in agreement with the trends obtained 
from Chapman's tests at 1/ = 1.5 and 2.0 for turbulent 


At = 


creased considerably with increasing Reynolds Num- 


boundary layers. 1.5 the base pressure in- 
ber; at 1/ = 2.0, the curves are almost flat; and at 
2.95, from the tests reported herein, the base pressure 
decreases with increasing Reynolds Number. Further 
increases in Mach Number would be expected to show 
greater decreases with increasing Reynolds Number if 
the trend shown in these tests continues. 


(b) Tests of Model with Roughness 


The experimental equipment permitted the attain- 
ment of fairly high Reynolds Numbers where the bound- 
ary layers naturally were turbulent. A few tests were 
made to evaluate the commonly used expedient of fix- 
ing transition to simulate high Reynolds Numbers and 
to examine whether or not the roughness (or band) or its 
The 0.25- and 1.00- 
Transi- 


relative height is the true variable. 
in. diameter models were used for these tests. 
tion was fixed after the cone-cylinder juncture by a 
band of masking tape 0.10 in. wide placed on the body, 
as shown in Fig. 2. Thickness was varied by using one 
or more layers of the tape, which was about 0.007 in. 
thick. 
carried out are shown in Fig. 5, and some typical schlie- 


The results of the four series of tests which were 


ren photographs are shown in Figs. 6a and 6b. The 
high base pressures at low Reynolds Numbers assoct- 
ated with laminar flow are, as expected, not obtained. 
The results obtained are, however, not related in any 
simple way to the data obtained at higher Reynolds 
Numbers for natural transition. Depending on the 
degree of roughness, below about 5 X 10® the roughness 
tests may fall above or below the smooth body data. 
The actual data depend critically on the ratio of transt- 
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Ry = 200% x 10° 


6 


Fic. 6a. Schlieren photographs of base pressure model with 
transition strip, D = 0.25 in., tf = 0.007, m = 2.95, for a range 
of Reynolds Numbers. 


Below 1 X 
10°, the addition of roughness gave base pressures con- 
siderably below those obtained in the smooth con- 
Between 1-3 X 10° (the ‘‘transition’’ region), 
the rough tests gave base pressures above the smooth 


tion band height to body diameter, ¢/D. 


dition. 


values but approached the smooth values as the rough- 





Fig. 6b. Schlieren photographs of base pressure model with 
transition strip, D = 1.00 in., ¢ = 0.023, M = 2.95, for a range 
of Reynolds Numbers 
ness was increased. Above 3 X 10°, the addition of 
roughness gave lower base pressures than the smooth 
The decrease in pressure 
The tests of 


models for all cases tested. 
depended on the thickness of the band. 
the 0.25- and 1.00-in. diameter bodies with bands of 
approximately the same ¢/D ratio gave comparable re- 
sults in the range where the tests overlapped. The 
tests of the 1.00-in. model with (/D = 0.023 and 0.045 
came together above about 15 X 10° and appear to ap- 
proach a constant value of about 0.37 for p,/ pi. Fur 
ther tests are needed to verify this, and tests at consider- 
ably higher Reynolds Numbers for a smooth body are 
needed to show whether such is also the case for na 


tural transition. 
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Logically, there seems to be little real reason to hope 
that the addition of roughness will completely simulate 
high Reynolds Number flows. The base pressure prob- 
lem is, as has been stated by several investigators, essen- 
tially a mixing problem. The dead-air region behind 
the blunt base must mix with the supersonic free stream 
through a layer that is the detached boundary layer 
from the body. The process must be a function of the 
thickness, as well as the character, of this layer. The 
attainment of high Reynolds Numbers by natural 
means results in a turbulent boundary layer that is 
thin. 

The introduction of a transition band will undoubtedly 
result in turbulent boundary layers, but considerably 
thicker than the natural turbulent layer. At the 
same time, the roughness tests show that thickening 
the layer gives decreasing base pressures, the same 
effect as was caused by increasing the Reynolds Num- 
ber. The results are not explainable by any of the 
present theories. 


(c) Effect of Sting Diameter and Length 

Chapman made a fairly complete study of the effect 
of support sting diameter and length at JJ = 2.9 in 
reference 1. For laminar flow, ratios of support diam- 
eter to body diameter below 0.6 and ratio of support 
length to body diameter over 2.8 gave no effect on base 
pressure. For turbulent flows, the ratio of support 
length to body diameter of over 2.8 also gave no effect. 
However, for ratios of support diameter to body diam- 
eter from 0.6 down to 0.3 (the smallest tested), the base 
pressure is still rising slightly. For the configuration 
used in these tests (a ratio of support diameter to body 
diameter of 0.25 and a length to body diameter ratio of 
3.9), it seems reasonable to assume that the results ob- 
tained are, at most, only slightly affected by the sup- 


port system. 
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CONCLUSIONS 
Preliminary tests were made at M/ = 2.95 of a simple 


cone-cylinder body of revolution over a range of Rey. 
nolds Numbers from 0.6 X 10° to 18 X 10°. The body 
was tested in the smooth condition and with transition 
fixed by several different roughness bands. On the 
basis of these tests, the following conclusions were 
reached: 

(1) The variation of base pressure with Reynolds 
number is characterized by three distinct regions 
(a) a ‘“‘laminar’’ region where the free wake appears to be 
laminar and where increases in Reynolds Number giye 
extremely large decreases in base pressure; (b) a “‘trap. 
sitional”’ region where transition from laminar to turby. 
lent flow seems to be occurring in the free wake and 
where increasing Reynolds Number causes a leveling 
off from the laminar region and then a sharp increase 
in base pressure; (c) a ‘‘turbulent’’ region where the 
boundary layer on the rear of the body is turbulent and 
where increasing Reynolds Number results in a slow but 
steady decrease in base pressure. 

(2) Extrapolation of data in regions (a) or (b) of (1 
will not give valid results at higher Reynolds Numbers. 
For the body tested, data above 6 X 10° is needed to 
predict base pressures accurately at high Reynolds 
Numbers. 

{3) The use of transition bands of the type tested on 
the body at low Reynolds Number gave results that 
not be correlated with data obtained on the 
The base 


could 
smooth body at high Reynolds Number. 
pressure depends critically on the ratio of band thick- 
ness to body diameter. 
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Note on Beams and Plates 


George Gerard 

Guggenheim School of Aeronautics, College of Engineering, New 
York University 

November 26, 1951 


INTRODUCTION 


T# CLASSIFICATION AS A BEAM OR AS A PLATE of a long rec- 

tangular structural element loaded by pure bending along 
the two shorter edges and free along the unloaded edges depends 
When anti 
clastic curvature is completely unrestrained such that Ry = 
-R,/v, then the flexural rigidity is EJ and the element is con- 


When the anticlastic curvature is completely 
1 


upon the nature of the cross-sectional distortions. 


sidered a beam. 
restrained such that Ry = 0, the flexural rigidity is E/(1 
and the element is considered to be a plate. 

The transition between a beam and a plate under bending has 
From a similar standpoint, 


py?) 


recently been examined by Ashwell.! 
Houbolt and Stowell? investigated the transition between a col- 


umn and a compressed plate. It is the purpose of this note to 
show that the classification of a structural element 
or a plate does not depend solely upon the geometric parameters 


Extension of the results given in references 1 


as a beam 


of the element. 
and 2 indicates that an element of a given size behaves initially 
asa beam and that the flexural rigidity may grow as the bending 
strain is increased. True plate behavior in which anticlastic 
curvature is completely restrained is approached only as a limiting 


case 
SYMBOLS 
b = width 
D = bending rigidity of a plate per unit width, D = 
Et?/12(1 — v?) 

E = modulus of elasticity 

| = moment of inertia 

R = defined by Eq. (11) 

V = moment per unit length 

W = total moment at cross section 
R = radius of curvature 

t = thickness 

u = deflection 

Y, ¥, 2 = rectangular coordinates 

€ = strain 

€ = outer fiber strain, « t/2R 
} = Poisson's ratio 

a = WR/EI 

¢ 3/4(1 y? V edit 

Ww 2ey/b 
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BENDING THEORY 


The cross-section distortional behavior of the structural ele 
ment can be examined by bending the entire element to a con 
stant radius Ry = R, and finding the total moment carried at 


each cross section 


The pertinent distortion relationships are 


Mx = D(A/Rx + v/Ry) (1) 
My = D(1/Ry + v/Rx (2) 
By assuming that the deflections are small 
1/Ry = —0O*w/Ox? (3) 
1/Ry = —O*w/dy? (4) 
The total moment carried at each section is 
~ Pb/2 
M=2 f Mx dy (5) 
0 


It is now possible to define a flexural rigidity factor, a, as fol 
lows: 


(6) 


a = MR/EI 
From Eqs. (1) and (5), the factor 
» *b/2 ( 
a = l T v 
(1 — v*)b 0 \ Ry 
-R/vy > Ry > 


Furthermore, since 7/EI = 


dy ( 


Note that between the limits, © ; the effective 
ness factor, 1 Z a Z (1 — v?)7! 
é/z and 1/R = e/sz, then, from Eq. (6), the flexural rigidity factor 
a = é/e. This is the effective strain due to bending relative to 
the strain which would result if unrestricted anticlastic curvature 


occurred 


Ashwell' has obtained the cross-sectional distortion as given 


by the following: 


sin w sinh w + (¢ cos w cosh w 


1 , . 
w/t=|3 l (¢ ) 
y 
(S 


where 
4 sinh ¢ cos ¢ = cosh ¥¢ sin ¢ 

C1 = 
9. 


sinh 2 + sin 2¢ 
+ ¥ 


For a given value of Poisson's ratio, the parameter ¢ is solely a 


function of «(b/t)?. The significance of this parameter is ap 
parent when a structural element of a given 6/f is subject to in 
strains the cross-sectional 
As 
the strain is increased, the central portion of the cross section 
tends to suffer less distortion finally retaining its original planar 
Further straining merely results in greater 
Thus, all struc 
initially 


creasing bending strains. At small 


shape is characteristic of unrestricted anticlastic curvature 


ity at large strains 
concentration of curvature at 
tural elements loaded in the manner described behave 


the free edge 


as beams and then tend to “grow” into plates 








208 JOURNAL OF THE 











Fic. 1. 









Dota of Ref 2 


Eg (ul) 








400 


ao or / 10 100 











€.(at)? 


Effectiveness factor for »v = 0.3. 
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Fic. 2. Flexural rigidity = 


FLEXURAL RIGIDITY 


The increase in flexural rigidity associated with this growth is 
given by the flexural rigidity factor, Eq. (7). By substituting 
the value of R given by Eq. (4) into Eq. (7), integrating directly, 
and noting that from symmetry dw/dy vanishes at the center of 


1 2vR ~) 
a= 1 —- (9) 
1 — »? b \dy/b2 


From Eq. (9), it is evident that the flexural rigidity factor of an 


the element, 


element of given dimensions bent to a radius R depends only upon 
the slope at the free edge of the element. 
By use of Eq. (8), the slope at the free edge was found to be 


dw 7 2y 3) 
dy/oj2 —- [3(1 — v?)]'/*\R 


Substitution of Eq. (10) into Eq. (9) results in the following 
simplified form for the flexural rigidity factor: 


*sin? yg + sinh? ¢ 


: : (10) 
sin 2g + sinh 2¢ 


— ky?)/(1 — v?) (11) 


a = (1 
where 
k *) sin? g + sinh? ¢ 
¢g/ sin 2¢ + sinh 2¢ 
Values of k for a range of the parameter ¢ are given in Fig. | 
The use of asymptotic expressions for the functions appearing 
in k reveals that, when ¢ is small, k = 1 and that, when ¢ is large, 
k ~1/g. The values of a which correspond to these limits are 
a = 1 when gis small, but note that a > (1 py)? 
This behavior of Eq. (11) indicates that the full flexural rigidity 
associated with a plate is approached only as a limiting case. A 


as g@—> &, 


similar occurrence was noted by Houbolt and Stowell? in an in- 
vestigation of buckling of columns and plates. 

By assigning a value to Poisson’s ratio in Eq. (11), a plot of the 
effectiveness factor as a function of the parameter e(b/t)? can be 
constructed by use of Fig. 1. This has been done for a value of 
vy = 0.3, and the results are shown in Fig. 2. 
plot the results given in reference 2 in terms of the same param- 


It is possible to 
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eters used in Fig. 2 
These data are shown in Fig. 2 to indicate that the trend for a'j 


similar in different types of loading. 


by a transformation of the coordinates 
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Critical Stress of Plate ( umns 


On Nonsteady Supersonic Flow About Pointed 
Bodies of Revolution 


John W. Miles* 

Fulbright Lecturer, Auckland University College, Auckland, New 
Zealand 

November 14, 1951 


N”’ rOO LONG AGO, Lighthill described ‘‘how certain difficulties 
in the application of the linearized theory to flow past bodies 
of revolution, which are continually being rediscovered, may by 
surmounted.’’! Intending no claim of discovery (or rediscovery), 
we feel it pertinent to reiterate one of these points in connection 
with Dorrance’s analysis? of the unsteady flow problem, which js 
subject to much of the same criticism originally advanced by 
Lighthill* in connection with Tsien’s analysis‘ of the steady-flow 
problem 

Notwithstanding Dorrance’s statement that his ‘‘solutions are 
all subject to the usual limitations of the linearized theory,” he 
retains terms in the velocity potential that are inconsistent with 
problem implies the 


the linearization. Linearization of the 


restrictions 


k=(Qi1): é6< l and (16)? << | la) 


ké < l and (kMé)? << 1 lb 
where k, / and 6 denote the reduced frequency, Mach Number 
and slenderness ratio, respectively (cf. Lin, Reissner, and Tsien’s 
discussion of the two-dimensional problem® and the heuristic 
treatment of reference 6. These conditions have since been es 
tablished rigorously as necessary, but not sufficient, for non 
planar bodies.’) Consequently, the linearized (with referencé 
to an expansion in 6) equation for the velocity potential reduces 
to Laplaces equation in every plane transverse to the undisturbed 
flow, and the relevant part of the potential for transverse motio: 
is then given by® 

-R*(x)} Op 
"(s, 7, @ §) é 


d(x, 7r,0,t) = 
r or r=R 


Dorrance’s results, going beyond this, are unjustifiably compli- 
cated [We do not claim that they would compare less favorably 
with experiment than their simpler counterparts but only that 
they are not rationally consistent with the (implied) initial 
hypotheses. | 

Finally, we remark that the “low-frequency” results given by 
Dorrance in Eqs. (31), (32), (36), (37), (41), and (42) are the linear 
ized results for all frequencies satisfying Eq. (1b) above but ar 


not new (cf. reference 6) 
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Generalization of the Theodorsen Function to 
Stable Oscillations* 


A.|. Van de Vooren 
National Aeronautical Research Institute, Amsterdam, Holland 


November 23, 1951 


SUMMARY 


It is shown that the Theodorsen function C(k), to be used in the expres 


sions for the unsteady aerodynamic forces, becomes infinitely large for stable 


motion when it is referred to an infinite length of the wake. A formula is 


denoting the generalization of C(k) to finite wake length 


given for C(k, S) 
This formula is 


s—which is valid for both stable and unstable motions 
required for the calculation of the aerodynamic forces, but the motion of an 
airplane can be obtained correctly when C(k) is only known for harmonic 
and unstable motion, even when the motion of the airplane turns out to be 


stable 


INTRODUCTION 
the 


I HE GENERALIZATION TO UNSTABLE OSCILLATIONS of 
C(k) function, originally introduced by Theodorsen' for 
harmonic motions, has first been considered by Jones?, who con 


cluded that the relation 


C(k) = H,®(k)/(H,@(R) 4 (1) 
retains its validity for complex values of k, provided Jmk < 0, 
which corresponds to unstable motions. For stable oscillations, 
however, Jones concluded that a certain contradiction, connected 
with the multivaluedness of the Hankel functions, would arise if 
Eq. (1 


were also retained in this case. 


Not agreeing with this result, Luke and Dengler recently pub 
lished a paper’ in which they take the point of view that Eq. (1) 
is valid independently of the sign of Jmk, while they make the 
> 0—i.e., that the reduced fre 


additional assumption that Rek 
In the present paper, the C(k) 


quency must be positive or zero 
function will first be generalized to finite values of the wake length, 
both for stable and unstable oscillations, and next it will be shown, 
that the C(k, s) function obtained in this way, is oscillatory di 


vergent for s — © in the case of stable oscillations 


GENERALIZATION OF THE C(k) FUNCTION TO FINITE WAKE 
LENGTHS AND ARBITRARY VALUES OF /mk 
Let the function C(k, s) be defined by 
. ' ‘ 
C(k = U dx ++ 
r 2 l 
i : x ] 
Udx (2) 
1 V x* l 


tuthor wishes to thank Martin Goland for correspondence on the 
the the author's 


* The 


Subject nee this has much advanced crystallization of 


ideas 
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FORUM 


which is seen to be the extension of Theodorsen’s definition to the 


case of a finite wake of length sb (6 = semichord). In Eq. (2) 
U denotes the vortex strength in the wake, which is a function of 
s — x, since the vortices are supposed to be at rest with regard to 


the fluid elements. For unstable and harmonic motions it may 
be assumed that the vortex strength approaches to a distribution 
of the type exp 7k(s x)ifs —x—~ oo, With this assumption, 
which will be verified later in the paper, it was shown in refer- 
ences 2 and 3 that Eq. (1) yields the value of C(k, s) if s ~ «= 
andIimk Z 0. 

However, the same reasoning may not be applied for stable 
motions, since in this case the initial disturbance influences the 
wake essentially. The function U is, in general, determined by 
the Kutta condition, which yields the following integral equa- 
tion for lL’ (see reference 1): 

stl ox +1 
U(s +1 —-x)dx = e*s (3) 
1 Vx? -— 1 


where e*** denotes the vertical velocity of the airfoil at the three- 
quarter chord point (in case of an airfoil with control surface an 
additionai complex factor A should be added to the right-hand 
side, but this does not modify the analysis) 

It will now be shown that, using Laplace transformations, the 
need not be solved explicitly for the calculation of 
Replacing x by y + 1, Eq. (3) changes into 


function L’ 
C(k, s) from Eq. (2). 


s 
y+2,, * 
: U(s y) dy = e* 
J 0 \ b 


to which a Laplace transformation may be applied with the result 


that 


L{U(s)} = Lihet*} 


——_ 


J. |s 
“ 


Making use of the well-known relations‘ 


K,(z) = f, e~ 2°9Sh 4 nosh nu du, Rez>0 
v0 


it follows, after substitution of s = cosh u 1, that 
f js+2i a 5 +2 
L< = e~ Ps ds = 
> ; 7 3 
e?} Kol p) + Ay p)t 
if Re p > O, and, hence, 
I l 
L}U(s){t = _— - (4) 
e(p ik) Kop) + Ki(p) 


This result suffices for the calculation of the integral in the num 
erator of the right-hand side of Eq. (2). Since this integral is 


equal to 


y) dy 
and since 
L}¥(s + 1)/VY s? 4 2s)! = e?K\(p 


it is found by aid of Eq. (4), that 


. y+] 1 Ki(p 
I U(s y) dy = : 

Jo Vy + 2y f p tk Kop) + Ki(p) 
Further, it is known from reference 2 |Eqs. (47) and (119)]| that 
for Re p > 0 

Ki(p l 
Set ns - Eafe) (5) 
pik p) + Kilp)} 2 
if k,( ~) = 2, where &,(s) is the Wagner function. Hence, 
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C which is the generalization of the C(k) function to finite y k p(s) 
: ake 
lengths and arbitrary values of /m k. p(s) 
By making use of Eq. (5), it can easily be verified, that for un The 
stable motions (Jmk < 0) and s ~ ~, C(k, s) approaches to th 
g : ‘ : the 
ry value given by Eq. (1). It will be shown in the following section 
that, for stable motions (Jmk > 0), C(k, s) will be oscillatory qj 
vergent in s. ; i.e 
Returning to the velocity in the wake, it follows immediately at 
; oe 1c 
from Eq. (4) that se 
Su 
D a E | co + ja eb(s—l dp tion: 
. E\B Us) = 
G ~~ JS 2at J jo (p — tk) {Ko p) + Ki(p)! (k 
C(k, 
where ¢ must be larger than the real part of any of the singular 33 
. ° . . {} de 
ties of the integrand. It has been shown by Schwarz? that th, : 
) function Kop) + Ai(p) has no zeros. Hence, the integrand has o 
C A only one pole—namely, at the point p = 7k. It follows that for 
Imk < 0, when the pole lies in the right half plane, the function 
U(s) will behave asymptotically like exp ks, which is in accord C(R, 
ance with the usual assumption. The same conclusion holds 
for/mk = 0 
A For Jmk > 0, the path of integration in Eq. (7) will be ABC 
(Fig. 1). The integral taken along the closed contour 
Fic. 1. Path of ie in ng complex p-plane, when ABCDEBFGA, considered as function of s, will have the char 
m a rs . Ngee 
- acter e'*5, Hence, for s > ©, the integral along ABC will be It w 
. equal to the integral along CDEBFGA, increased by (O(e" In F 
I J th — dy ra I Lik(s)} 4 It will now be shown that the integral along CDEBFGA is of for | 
‘ / . a 7 — “ l . > . . 
i. 0 Vy? + 2y { y higher order than e'®S and, hence, will determine the asymptotic (dat 
character of U(s). only 
ikL {et®S! Li ki(s)} - ; ’ . a 0.90 
2 Since for |p| ~ o, the function Ay(p) + A,i(p) is of order T 
a e ?’/¥ p, it is seen that the integrals along CD and GA vanish the 
a For small values of p there exists ae 
. y+ a 
} U(s — y) dy = ip tion 
0 Vy? + 2y Kop) + Ki(p) = }Jiip) + iJi(ip)} (c a ) 4 witl 
=, i : 
= 1 ki(s) + tk ki(s —o eke dg ( l 
° 6 > R(p Atl 
Substitution of this result and Eq. (3) in the expression for C(R, s) for. 
yields where R(p) is regular in p near p = 0. After substitution into 
' i the integral and multiplication of numerator and denominator 
C(k, s) = J ky (s)e ikS 4 ip k(x )e — tks dx § (6) by p, it will be clear that the integral along EBF also vanishes A 
2 | 0 f The integrals along DE and FG yield the 
an } 
tior 
be 
0 ebls 1) dp flap 
+ 
I 
i . I ‘ . . . ip 
(p — ik) ~ {Jhip) + iJip)}( C + In ,) + RP) 
Pp 2 
» . x bls 1) dp 
0 l , eee . ip 
(p ik) iJo(tp) + iJi(tp)} | C + In — 271i} + R(p) 
p 2 C 
P= eP'S~1) Dei p2{ Sip) + iJiip)} dp 0 
= 0 
0 . f “ayy . ip { “2\t . ip 
(p — ik)| 1 — p {Jolip) + Wi(ip)s | C + In © ) + pR(p) 1 — piJolip) + Si(tp)s | C + In > — 2er] 4 pR(p : 
. > Xu 
C 
C 
C 
Since the behavior of this integral for large s is determined by Hence, the conclusion is reached that for a stable motion the 
the values of the integrand for small |p|, the integrand is ex- vorticity diminishes as s * and not as e**S (Im k > 0) ' 
panded in terms of powers of p. This leads to the result that the 
last integral may be replaced asymptotically by CALCULATION OF C(k, s) IN THE CASE OF A 
: STABLE MOTION 
—_ ePs p? dp : 
° For the purpose of calculating C(k, s) from Eq. (6), the fune- Fr 
which is of order 1/s°. tion k,(s) will be approximated as follows 
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p(s) = 2 — 0.38¢ 0.04555 _ () §7¢ — 0.3005 if0<s<s, (8 
i\ = ; = —_ 3) 
p(s) = 2 — (2/s) ifs; <s 
The value s is determined by the relation 

9 

= »« 0.3008 


0.045558) 4 ().67e 


= ().33e 


§ 
1 


s; = 43.131. It is seen from Table 1C of reference 2 that 


Bis 
the second approximation, which has the right asymp 


fors > § 
totic character, is better than the first. 


Substitution of Eq. (8) into Eq. (6) leads by elementary reduc- 
tions to the result 

0.671Rk 
ik + 0.300 


j 0.33ik 
l ») - 


Ce, s) = 2 as ik + 0.0455 


(.33(0.0459 ) ik 0.0455) s 0.67(0.300) (ik + apn | 
¢ 
ik + 0.0455 ik + 0.300 { 
ifs <8 
‘ Pe kx 
l 2 ‘ ‘e (9) 
Clk, s) = 2 ets _ Dik dx 
2 Ss - x 
0.330k f a (ik + 0.0455) 5s; t 
=e ‘ 
ik + 0.0455 
0.67 ik i _ ,— (th + 0.800)s11 ee 
‘ ra 1 
ik + 0.300 


It will be clear that C(k, s) will be oscillatory divergent if Jm k > 0. 
In Fig. 2 the real and imaginary parts of C(k, s) have been plotted 
for k = 0.05 (harmonic oscillation) and k = 0.05 + 0.0133977 
(damped oscillation). The approximation . (8)] introduces 
only a slight error, since the exact value C(0.05, 
(0.90901 — 0.130647, while Eq. (9) yields 0.90828 — 0.133607. 
The values of C(k, s) following from Eq. (6) could be used for 
the calculation of the aerodynamic forces by the usual formulas 
For example, the circula 





») is equal to 


as originally given by Theodorsen.! 
tional part of the lift force due to a damped translational motion 
ARS ts given by 


with a vertical velocity w(s) = ¢ 


Lis) = —2mrpvbC(k, s)w(s) 


Although C(k, s) tends to infinity, the lift force vanishes of course 


for s > 


CALCULATION OF THE MOTION OF AN AIRPLANE 


Although the C(k, s) function is required for the calculation of 
the aerodynamic forces, it will now be shown that the motion of 
an airplane can be calculated exactly if the orthodox C(k) func- 
tion is known for 7m k < 0, even if finally the motion turns out to 
be stable. To this purpose a wing model that is restricted to 
flapping motion will be considered 
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be the dimensionless equation of motion, where m denotes the 
mass, w» is the natural frequency, and a dot denotes a differenti- 
At the moment s = 0, a sudden disturbance 


Application of a La- 


ation towards s. 
acts on the model [2(s) = 0 for s < 0]. 
place transformation to this equation leads to 


m2(0) 


(10) 
m(p? + wn?) + rp?L}\ki(s)} 


L}\2(s)} = 


Hence, 2(s) will eventually vanish if all poles of the right-hand 
side of Eq. (10) have negative real parts, but it will become in- 
finitely large if there are poles with positive real parts. Thus, 
the stability question is reduced to an investigation of the 
zero’s of the denominator of Eq. (10) in the right half plane ( Re p > 


0). Substitution of p = jk makes this denominator equal to 


m( —k? + wy?) + 2xfkC(k), Imk <0 (11) 


for a stability analysis the function C(k) needs only to be 
For systems with 


Hence, 
known for unstable motions (see reference 6) 
more degrees of freedom, Eq. (11) will be replaced by the sta- 
bility determinant. 

For stable motions, 2(s) follows from Eq. (10) as 


ic 
f L\2(s)} e?S dp 
ja 


where the integration curve is the imaginary axis, together with a 
half circle at the side of positive real parts to avoid the singular 
In this way the resulting stable motion can be 
and un- 


2s) = 


point p = 0. 
calculated when C(k) is known for Jmk < 0 (harmonic 


stable motions). 
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SUMMARY 


It is shown that Theodorsen’s circulation function cannot be applied to 
stable damped motions by simply replacing the imaginary parameter by 
It is pointed out that the lift deficiency function of von Kar 


applied to 


complex one 
man and Sears is the fundamental relation that may be directly 
all nonuniform flow problems. The reciprocal relations between Wagner's 
step function response and Theodorsen’s steady-state sinusoidal response are 


correctly derived directly from the lift deficiency function 


INTRODUCTION 


VooREN! AND OTHERS have objected to the applica 
function C(iw) to forced 


‘ J AN DE 
tion of 
motion of the 


variable p defining the motion is negative, 
stable motion, and C is used by simply replacing tw by p. 


Theodorsen’s circulation 
form e”', wherein the real part of the complex 

corresponding to a 
Vari- 
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ous other papers have been written?~4 to show the equivalence 
or reciprocity of the step function response of Wagner to Theo- 
dorsen’s circulation function for steady-state sinusoidal response 

The present note uses the work of Sears‘ to show the limitations 
involved in these relations. It is pointed out that the lift defi 
ciency function of von Karman and Sears® is the fundamental 
function that satisfies all the reciprocity relations and is applic 
able by means of the Duhamel integral to all motions, including 
even gust loadings, as long as the deviations from steady state 
may be superimposed linearly. 

In some simple cases the generalized Theodorsen circulation 
function is exactly equivalent to the lift deficiency function but is 
not directly applicable, for example, to the damped or stable mo- 
tions or to the gust-type vertical velocities. 

From reference 5 we write, for the lift on a flat plate of infinite 
span in incompressible flow at a constant forward speed, l’ = 
s/t, 


L(s) = Los) + Lil(s) + Los) (1) 
where 
0/2 
Los) = oU { yx, S) dx = pUT (s) (2) 


is the quasi-steady lift computed for the instantaneous velocities 
by assuming they are steady state and thereby neglecting all the 
effects induced by the wake, 


ad ftc/2 
Li(s) = pl f x yo(x, 5) dx (3) 
ds c/2 


is the so-called apparent mass contribution that depends only 


upon the instantaneous acceleration, and 
es 
LAs) = —L)o(0)¢(s) — | o(s a)ly'(a) da (4) 
0 


is the entire effect induced by the two-dimensional wake, as 
sumed a flat rigid sheet shed by the plate. In Eq. (4) we obtain 
L, by means of the Duhamel integral and ¢, the indicial response 
or lift deficiency, giving the transient response of 2 to a step 
increment of Zo. From reference 4 we have the modified Laplace 
(or Carson ) transform of ¢ as 
ocx , F Kop) 
or) =» f e "§ &s)ds = — : (5) 
J0 Kop) + K\(p) 


APPLICATIONS 
For a sudden change in Ly we have Wagner's step response R, 
defined by Eqs. (1) through (4) as 


Ri(s) = L(s)/Los) = [1 o(s)| (6) 


For the steady-state sinusoidal motion, the lift due to circu 
lation is given by Theodorsen’s circulation function 


Cli Ki (iw) Lo(s) + Los) 1 

(iw) = ae eee = = | (7q@) | (7) 
Ko(iw) + Ki(iw) Li(s) = 

in agreement with Eqs. (2), (4), and (5). The real part of p 


may be taken as zero in Eq. (5), since all the singularities of 
¢(p) lie on the negative real axis, as shown by Sears.‘ 

For an unstable divergence, where uw (the real part of p = 
u + tw) is positive, we write Eqs. (2) and (4) for the nonoscillatory 
real divergence 


Lo(s) = Lo(O)e*" (8) 


LAs) = LA(O)d(s) — Lo(O)u f os a )e#" do 
u 


8 
= LO )b(s) Lol su f, os ole #9) d(s a (9) 
e 0 
or as a limit when s > 
Li @)/Lilo) = —o(p) (10) 


so that 


[Lolo ) + LA )]/Lilo) = [1 — o(u)] = Cig 
This proves Theodorsen’s circulation function may ty used witt 
iw replaced by uw 2 0 only in the limit as s; > “oa 
Obviously, this is not valid for the case of stable. lamped 9 
cillations where w < 0. For this application writ« = 1 
iw, where 8 >0. Then 

Lis) >= Li O)e Bs » 

es 

Ls) = LAO)P(s) + Li(O)B { ds ae" dg 

° 0 

Now, unlike Eq. (9), Eq. (13) cannot be put directly in the for 


of the Laplace transform. However, the Laplace transforn 


lL» may be obtained by the convolution or Borel produet rule 


L.(p B 
= »(p) + o(p Bp 
LO) ’ p / 
Since 
Lil p) ) 
— = opp = / - 
LAO p B _ 
therefore 
/ / Ap c—¥ 
wip) + wp) = AO = 
Kip) + KAi(p) p+ B 


This proves that Theodorsen’s circulation function cannot give 
the lift directly for stable damped motions where the real part 
of p is negative, since the inverse transfer of Eq. (16 depends 
upon L)(p) However, the inverse Laplace transform of Eq 
(16) does define a new function that fulfills the same purpose for 
each given motion. This is the type of function that is tabulated 
in reference 6 and could conceivably be called a generalized 
Theodorsen circular function for this given type of stable damped 
motion Exactly the same type of generalized function would 
be required for 8 < 0, corresponding to Eq. (11), the divergent 
types of motion, when s is finite and p = uw + iw 

It is important to note, however, that the lift deficiency function 
directly gives L, for a step response and, therefore, for any arbi- 
trary motion by means of Eq. (4). Then, as long as Lo’( ~) = 0 
or s is finite, we may combine Eqs. (4) and (6) for any uniform 


motion following a sudden step as 


(es 
Lis) + Los) = LA0)Ri(s) 4 [ Ri(s a)Ly (oe) da 17 

J0 
For a gust-type vertical velocity, Eqs. (1) through (5) are still 
applicable, although Eq. (17) is not. For example, for a sharp- 
edged gust, Sears! has already shown how simple it is to obtain 
the lift from Eqs. (1) through (5) by using the convolution prod- 
uct rule as in Eq. (14). In this way the Kussner gust response R; 
may be easily obtained,?:* and this response to a sharp-edged 
gust may be used in place of R; in Eq. (17) to determine the lift 

upon entering any smoothly graded gust 


THE RECIPROCAL RELATIONS 


The lift deficiency function directly shows the reciprocity 
relations between the step response and the steady-state sinu- 
soidal response as may be seen by writing Eqs. (6) and (7) in 


conjunction with the inverse Laplace transform as 


1 ¢ sac _ dw 
|] Ri(s)|] = o(s) = / jl C(1w ) |e’? 18) 


«Tle Ww 


This avoids many mathematical difficulties, since ¢(~) has all its 
singularities on the negative real axis.‘ For example, directly 
from Eqs. (5) and (18) we obtain the following reciprocity rela- 


tions, which replace the incorrect results given in reference 2 


F(w) = 1 w { |] Ri(s)| sin ws d 19) 


. 


where 
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(e «x 
G(w) = - of l Ri(s)| cos ws ds (20) 
0 
| i y : c : dw 
Ris) == + |F(w) sin ws + G(w) cos ws] (21) 
2 rl a 
where 
C(iw) = F(w) + iG(w) 
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N A RECENT PAPER on the generalization of Theodorsen’s lift 
function C to growing and damped oscillations, Luke and 
Dengler! criticize some of the analysis given in an earlier report 
of mine on the same subject.2 It should perhaps be pointed out 
that their comments refer only to the concluding paragraph of 
Appendix I, where the problem of damped oscillations is tenta 

tively discussed, and do in no way reflect on the validity of the 
main part of the paper which deals mainly with growing oscilla 

In this particular report it is shown that for unstable mo 

> 0 and S = distance traveled 


tions 
tions of the type e‘*" w)*' where pu 
in half chords, the corresponding lift function may be defined in 
terms of Hankel Functions by the relation 

C(w ip) = Ay (—2)/[Ai(—2) 1H) ( —s)] (1) 
Since the Hankel Functions are multiple 
By interpreting 


where z w lp 
valued, difficulties arise in their interpretation 
sas suggested by McLachlan* and by using his formulas for 


Hy (—s) and H,\” (—=), Eq. (1) yields 


C(z) = AHy!?(2)/ [Ay ?(s) + 1H9"?(s) | (2) 


when » > 0. This gives the usual expression for C when p = 0 is 


assumed—namely, 


C(w) = Hy?)(w)/(Hi?'(w) + 1Hy?(w) (3) 
As a matter of mathematical interest the limiting form of Eq 
1) when y is negative McLachlan’s 
formulas were again used, but a different formula for C(w) was 
obtained for this case. Luke and Dengler maintain that the prin 


ciple of analytic continuation should be satisfied and show that, 


wi 


was also considered 


by interpreting —= differently (—z = ze™' instead of —s = se 
as used in references 2 and 3), (3) can be derived when p» — 0 for 
damped oscillations. They also maintain that Eq. (2) is appli 
cable generally for both damped and growing oscillations. The 
writer feels, however, that consideration of C(z) has no real signifi- 


cance in relation to the problem of damped oscillations, since 
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the analysis on the basis of which Eq. (1) was derived breaks down 


when wp < 0 and S—~ o. For this case the lift function must 
depend on the initial disturbance and the distance traveled, and 
the only satisfactory way of treating the problem is to make use 
of the Wagner &,(s) function and to apply the superposition 
method. This method was suggested and used in reference 2 

In a recent note van de Vooren‘ has considered the problem of 
stable oscillations in detail, and his conclusions clearly indicate 
that the tables of C(w iu) for negative values of « given by 


Luke and Dengler are of no practical value 


REFERENCES 


Tables of the Theodorsen Circulation 


and Dengler, Max A al 
Aeronautical Sciences, Vol 


' Luke, Yudell 
Function for Generalized Motion, 
18, No. 7, pp. 478-483, July, 1951 

2 Jones, W. P., 
R. & M. 2117 

McLachlan, N. W 

* van de Vooren, A 

Oscillations, N.1..L.. Amsterdam 


Journal of the 


Aerodynamic Forces on Wings in Non-lUni m Motion 


Bessel Functions for Engineers 


I., Generalization of the Theodorsen Function te 


Notes on the Calculation of the Response of 
Stable Aerodynamic Systems 


M. A. Dengler, M. Goland, and Y. L. Luke 
Research Engineer, Associate Director for Engineering, and Re- 
search Mathematician, respectively, Midwest Research Institute, 
Kansas City, Mo. 
November 28, 1951 
— THREE PRECEDING READERS’ FORUM NOTES are indeed 
interesting and informative and appear to throw suspicion 
on certain of our past techniques used to calculate the stability 
and response of stable aerodynamic systems and, in particular, 
on our use of the exponential form of the approximation for the 
Wagner growth-of-lift function and on our definition of the gen 
eralized Theodorsen circulation function in the stable half plane 
The present remarks, unfortunately, offer no definitive proof of 
our approach but do show that it is adequate for engineering 
purposes and suggest that many facets of the general theoretical 


question remain as yet unanswered 


(1) In our past work, primary interest is attached to calcu 
lation of the stability and response of aerodynamic systems, par 
ticularly stable systems. It is significant that in a variety of cases, 
where the Wagner growth-of-lift function was used in its exact 
tabulated form (the equations of motion then being integrated 
numerically and these results being compared with those from 
an operational solution based on an exponential approximation 
for the Wagner function), excellent agreement was obtained in 


both stability and response studies 


To fix this idea more precisely, consider a two-dimensional 
wing permitted to flap only, with deflection Z. The wing is 
elastically mounted in a uniform wind stream and, with regard to 
transverse motion, is initially at rest. A step-function lateral 
force is then applied, and the resulting flapping motion is studied 
According to a nondimensional scheme, the equation of motion 
for the wing can be written exactly as 

2: T 2: 
Zt+a f, K(r — u)Z(u) du + wn?Z = wp? 
where 7 is time, A(7r) is the Wagner function, w, is the natural fre 
quency of the suspension in the usual reduced form, and the ap- 
plied force is arbitrarily set equal to w,?. The constant a depends 
on the wing lift-slope-coefficient wing chord, air density, and wing 


mass 
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For the purpose of this note, solutions for Z were obtained by 
two methods for the two cases w, = 0.5 anda = 0.05 and 0.20. 
The choice of w, was governed by the desire for a frequency that 
is typical of aerodynamic systems where unsteady effects are 
significant; the value a = 0.05 affords a lightly damped system, 
and a = ().20 results in fairly heavy damping. 

The solutions for Z were obtained (a) by an extremely accurate 
small-interval numerical integration of the equation of motion, 
using the Lagrange coefficients and employing the exact tabulated 
values for K(7r); and (b) by an operational solution utilizing the 
three-term exponential approximation for A(r) suggested by 
R. T. Jones. 

In both cases, the Z-time histories determined by the two methods 
were indistinguishable from each other within accuracies beyond 
those of engineering interest. The detailed results are not pre- 
sented here merely in the interests of space brevity. 

These and other similar studies suggest that the exponential 
approximation for A(7) is adequate for calculations pertaining to 
stable, as well as unstable, systems. 

As a second point of interest, consider the definition of 
the 


(2) 
the generalized Theodorsen circulation function C(k) = n 
stable half plane. It is first to be noted that the definition 


C(k) = kil jk)/ [Rol jk) + ki(jk)| 


affords results that are consistent with the use of the exponential 


K(r) approximation. Stability roots for stable systems calcu- 


lated on either basis are in essential agreement with each other. 


Next, consider the definition 


t+1 he . 
f xe —I®* dx /(x? — 1) 
Pe 


T+1 
f (x + 1)e—I** dx/(x? — 1)' 


. 


C(k, r) = 


lim C(k, r). 
cl 

numerator and denominator of this expression are oscillatory 
Now, write both the numerator and 
Let 
each term of the series represent the area under one-half of a wave 
along the r-axis; the series is made up by orienting each term so 
that the areas follow in the order of their occurrence along the 
The two 


and let C(k) = For convergent motions, both the 


divergent as 7 increases. 
denominator as infinite series, each constructed as follows: 


r-axis. (Thus, adjacent terms have opposite signs. ) 
series so formed are alternating, divergent in character. 

Now, sum these two series by means of Euler’s transformation 
method, and then take the ratio of the two sums. For C(k) de 
fined in this manner, the result will be equivalent to the definition 


C(k) = Ki/(Ko + K;) 


Thus, consider jk = —0.0868 + 0.49247. For this value, the 
Euler sum of the numerator series is 0.732 + 2.458); for the 
denominator series it is 0.650 — 1.6997. Hence, C(k) = 0.594 — 
0.1647. According to C(k) = Ki/(Ko + K,), the value should 
be C(k) = 0.585 — 0.171). 
C(k) values are due to the use of the trapezoidal rule in the nu- 


The differences between the two 


merical integrations involved in forming the Euler summable 
series. 

With regard to the note by van de Vooren, it is of interest that 
a treatment of his Eq. (6), based on a similar infinite series con 
struction and Euler summation, leads to a finite definition for 
C(k), with a value in agreement with that given by the K,/(Ky 4+ 
K,) fraction. 

These results suggest that, mathematically at least, a finite 
value of C(k) can be constructed for stable motions. The physi- 
cal problem of interpretation remains, and these remarks only 
serve to indicate that as yet all the answers to this elusive ques- 
tion are not at hand. 
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Comments on ‘‘Shear Center for Thin-Walle 
Open Sections Beyond the Elastic Limit” 


Herbert R. Lissner 
Professor and Head, Engineering Mechanics Department, Way, 
University, Detroit ‘ 


November 27, 1951 


ly regard to the paper by Handelman' I should like to point out 

that for the problem described no shear center can exist 
which is of significance to the designer or stress analyst Th 
physical significance of shear center is that, if the applied Joga; 
pass through this point, bending will occur without twisting 
that is, no torsional stresses will be present in the beam. Frop 
the analysis presented, it is immediately apparent that, for th 
portion of the beam which remains entirely elastic, torsion wij 
In fact, for 
the beam partially stressed beyond the proportional limit gy 


be present and torsional stresses will be produced. 
loaded through the described ‘‘shear center,’’ all cross sections of 
the beam are subjected to twisting. 


Therefore, the point determined in reference 1 should definitely 
not be referred to as the shear center, since it does not represent 
the physical properties of that point. 
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Shear Center for Thin-Walled Open Sections Beyond 
Vol. 18, No. 1] 


'‘Handelman, G. H., 
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Pressure Distribution on an Airfoil in 
Nonuniform Motion* 


S. Neumark 
Royal Aircraft Establishment, Farnborough, Hants., England 
November 15, 1951 


I A PAPER PUBLISHED IN 1938 by von Karman and Sears,' the 
following general formulas for lift and moments on a thin 
two-dimensional airfoil in an unsteady vertical motion have been 


given: 
L = pUT» — p(0/dt if 1 Yor dx + uf, ydt/V # 1 
l 1 1 l 
M= pul 1 Tox dx — 5p(0 ot if v0 x? — : dx (| 





l ; a 
pu, ydt/V e — 1 


yo(x, t) being quasi-steady vorticity on the airfoil, Ip(t) being 
The 


first terms in Eqs. (1) represent quasi-steady parts, the second 


quasi-steady circulation, and y(é, ¢) being wake vorticity 


ones are contributions of apparent mass, and the last ones are 
due to the vortex wake. The formulas were derived by con- 


sidering the momentum of the entire vortex system 


It could be expected, by examining these formulas, that an 
analogous expression should exist for the resultant pressure dis 
tribution TI(x, ¢) along the airfoil, because 1 and M are simply 
integrals, from x = —1 to +1, of and xl, respectively. Such 
a formula, however, cannot be derived simply by von K4rman 
and Sears’s method, nor is it easy to obtain it from existing com- 
plicated expressions of Kiissner? and L. Schwarz. 


Chief British Ministry of 


Supply for permission to publish this article 


* Acknowledgment is made to the Scientist 
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I have succeeded in deriving the long anticipated formula by 
directly applying the Bernoulli theorem for unsteady flow (and 
avoiding all doubts as to convergence of integrals along the vor- 
tex wake, especially for the case of steady-state oscillations), in 


the following form: 


* 


x 


i . | —- e dé 
Yon dx + pl : : - = (2) 
1 ts \ 1+4J1 Ve 


Eq. (1) and 


pUyo + p(0/OF) 


nm 


re the roles of the three terms are the same as in 
on(x, t) denotes the noncirculatory part of the airfoil vor- 
The latter, as discovered by Betz‘ and expounded by 


whe 
where 7 
ticity ‘Yo. ; : 
Schwarz, can be determined by a remarkable “solving kernel in- 
version” of the integral equation of vortex sheet 
1 
l yo(s, t) ds 


w(x,t) = (3) 


2r -—1 £=-—S 


where w represents the downwash distribution on the airfoil, 
being a starting point of the theory and known in each particular 


case of unsteady motion. The Betz solution may be written 


T w(s, OWI 2d 
WS, S* as (4) 


| ad, 


rV 1 —x? rV1—x?/-!1 
with the circulation I'\(¢), satisfying Joukovsky condition at the 


trailing edge, given by 
al ha | 
f ' ve dx=2f wis, t)V¥(1 +s) /1 —s)ds (5) 


Splitting up yo into “circulatory” part yoc and ‘‘noncirculatory”’ 
one yon has proved essential for establishing the formula for pres- 
sure distribution |Eq. (2) ]. 


. 

Eqs. (1) are now easily obtained as simple consequences of Eq. 
The respec- 
As to the 


(2), by integrating 1 and xII from x = —1 to +1. 
tive first and third terms are obtained immediately. 


second terms, we have 


1 x 1 x 
f pdx (0 ot f 1 Yon dx) = p(o ar) f dx ff , Yon = dx 


(integrating by parts) = 


1 1] 
—p(O a) ff 1 Yon x dx = —p(0/dt) f 1 vox dx (6) 
and 
1 (Px 1 x 
f px dx (0 ar f , Yon dx) = (2/ar) ff i* ax f Yon = ax 
: | 1 
(integrating by parts) = — = p(o a) ff 1 Yon x?dx = 


l 1 1 l 
— 5 A(0 o(f yox? dx - f 1 Ye x? ax)= — 5 p(0/dt) 
.2 " l . I . .2 l 1x - 
, Yor? dx - 9h = — 30 or) oat x? —3 dx (7) 


so that the connection has been proved 

Eq. (2) assumes a particular, more definite, form in the case of 
steady-state oscillations. Denoting angular frequency by n, and 
reduced frequency by k = /U, we have the well-known relation- 


ships 
Oyo /Ot = inyo, 7 = — Toe ~ #8 (Ko + K;) (8) 
and, hence, Eq. (2) becomes 
x , | 
Io Ko = FF 
Il = pl (> + ik Yon dx — ( < -) . - (9) 
xr\Ko + Ky Vv +x 
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As an example, let us consider vertical translatory oscillations 
The downwash velocity w, whiie oscillating in time, does not vary 


along the airfoil; therefore, from Eqs. (5) and (4), 


Ip = 2xrw, Yon = —2wx/V 1 x?/ 
4 (10) 
yo = 2wvV (1 — x)/(1 + x) \ 
and the pressure distribution becomes 
a Ky i-£ 
Il = 2pUw + ik~WV/1 — x? (11) 


Ko+KiNi+x 


This consists of twoterms. The first one is the well-known solu- 
tion for steady vertical translation, multiplied by the complex 
‘*Theodorsen’s® function’? K,/( Ky + K,) of ik—i.e., decreased in 
magnitude and lagging in phase behind the airfoil oscillation 
according to the vector diagrams in Fig. 6 of reference 1 (lower 
curves): the second one is elliptic, equal to that for the quasi- 
steady effect of angular velocity kw about mid-chord axis multi- 
plied by i—i.e., leading the airfoil oscillation by 90°. This is a 
rather striking picture of the oscillating pressure distribution. 
The pictures of lift and moment follow immediately. 

It is easy to find solutions for other types of oscillations—e.g., 


for the rotary ones. 
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Single Degree of Freedom Flutter of an Aileron 


H.L. Runyan, H. J. Cunningham, and C. E. Watkins 
Langley Aeronautical Laboratory, N.A.C.A. 
January 11, 1952 


IT A RECENT TREATMENT of the problem of aileron flutter in a 
single degree of freedom by K. P. Abichandani and R. M. 
Rosenberg, it was concluded that the flutter speed of an aileron 
hinged ahead of (approximately) the 34 per cent aileron chord 
is imaginary. This conclusion was based, however, on the failure 
of these authors to realize that the aileron divergence speed, 
rather than the flutter speed, is imaginary for these hinge con 
ditions. In the notation of Abichandani and Rosenberg' the 
divergence speed is given by Eq. (1) 


| 
va = \ Qnpb? (ce ‘ ee) 7 (7; - os) . ft | 
3? 2n? 2x? 2? 


where Cg is the coefficient of torsional restraint, / is the distance 
between aileron leading edge and hinge axis, the 7”s are defined 
by Theodorsen and Garrick? and the ¢'s are defined by Kiissner.* 
Both quantities, 


( ¢1¢31/7?) + (935/22?) 


and 
(7; — TsT\0)/2r?|] + (TT 12/22?) 


in Eq. (1) are positive. Thus, if / is sufficiently small, the de- 
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nominator of Eq. (1) is negative and vq is imaginary and phys- 
ically does not exist. 

With this interpretation, the conclusions of Abichandani and 
Rosenberg that the ratio of flutter frequency to aileron natural 
frequency is less than unity should be reversed. That is, the 
flutter frequency should be greater than the aileron natural 
frequency for sufficiently small /. This is in agreement with 
the results given by Runyan, Cunningham, and Watkins,‘ and 


Runyan. 
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A Note on Conical Flows in Airfoil Theory 


John W. Miles* 

Fulbright Lecturer, Auckland University College, Auckland, New 
Zealand 

November 14, 1951 


"eam IN A RECENT NOTE in these columns and discussing 
Lampert’s? use of conical flow methods in subsonic airfoil 


theory, states that (in this connection) ‘‘ . . . conical flows . . . are 


merely a mathematical tool for solving nonconical flow problems, 


AERONAUTICAL 


-MARCH, 1952 


SCIENCES 


and it must not be construed that they have any significance in 
themselves.’’ While we agree with Goodman as to the necessity of 
satisfying the appropriate boundary conditions at infinity, we wish 
to point out that the use of homogeneous functions in attacking 
Laplace's equation was proposed originally by Green? wel] over a 
century ago and that the work of Donkin‘ on the same subject jg 
(Donk. 
in’s work has been noted previously by R. T. Jones, who suggested 


large measure anticipated modern “conical flow”’ theory. 


the use of homogeneous functions in the work of reference 2). 
Moreover, while the concept of conical flows in supersonic gir. 
foil theory is properly attributed to Busemann,° the general ideg 
of introducing homogeneous functions in attacking the waye 
equation is due to Bateman.® 

“conical flow,” 


the term so des 


criptive in supersonic applications, is perhaps inappropriate 


Finally, we remark that 


in subsonic work, but the antecedent requirement to the success- 
ful introduction of homogeneous solutions in any physical prob. 
lem is invariance, at least for individual parts of the field, toa 
scale transformation (cf., Hayes’s’ treatment of conical flows and 
Lamb’s discussion of diffraction by a half plane*). 
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